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Infrared  spectroscopy  is  used  to  investigate  the  effect  of  transport  currents  on 
the  ab-plane  superconducting  state  properties  of  thin  c-axis  oriented  YBa2Cu306+x 
films  on  sapphire  and  magnesium  oxide  substrates. 

The  films  were  first  characterized  in  both  the  normal  and  superconducting 
states,  through  optical  reflectivity  and  transmittance  measurements  over  a  wide  fire- 
quency  range  (40-350  cm~^)  in  the  far-infrared.  We  observed  fairly  high  reflectivity 
and  low  transmittance  in  the  superconducting  state  at  firequencies  below  150  cm~^ 
This  is  in  contrast  to  the  lower  reflectivities  and  higher  transmittances  which  are 
observed  in  the  normal  state  over  the  same  frequency  range. 

It  is  weU  known  from  nonlinear  transport  studies  on  thin  superconducting  films 
below  their  transition  temperature,  Tc,  that  large  currents  exceeding  the  critical 
current  Ic  destroy  superconductivity,  triggering  a  transition  to  the  ohmic  lossy 
regime  typical  of  the  normal  state.  We  studied  the  effect  of  transport  currents  on 
the  superconducting  state  transmittance  in  the  far  infrared  region  below  350  cm~^ 
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The  first  set  of  spectra  was  taken  at  75  K—  about  10  K  below     for  both  films— 
and  a  set  of  spectra  was  also  taken  at  low  temperature  (40  K)  for  the  film  on 
sapphire  substrate. 

The  films  on  magnesium  oxide  showed  no  electro-optic  response  up  to  very 
high  currents,  while  at  high  currents,  the  film  on  sapphire  showed  an  increase  in 
transmittance  at  low  frequency  accompanied  by  a  suppression  of  the  transmittance 
at  higher  frequencies.  At  lower  temperature  a  similar  effect  is  observed  in  the  film 
on  sapphire,  even  at  currents  much  lower  than  the  critical  value.  We  believe  the 
changes  observed  are  spurious,  and  may  be  attributed  to  ohmic  heating.  This 
has  been  further  verified  by  simulations  of  the  dielectric  response  using  a  simple 
two-fluid  model  of  the  superconducting  state  with  a  current  dependent  superfluid 
oscillator  strength. 
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CHAPTER  1 
INTRODUCTION 

High  temperature  superconductivity  (HTSC)  in  the  layered  cuprates  has 
been  a  central  issue  of  contemporary  condensed  matter  physics  for  more  than  17 
years  [1].  These  materials  are  interesting  not  only  because  of  their  novel  electronic 
properties,  but  also  because  of  the  enormous  potential  which  exists  for  fabricating 
many  practical  devices  incorporating  them.  Infrared,  microwave,  terahertz  and 
optical  spectroscopies  have  been  used  to  study  the  nature  of  charge  transfer,  the 
dynamics  of  free  carriers,  phonons,  the  Josephson  plasma  resonance  and  many 
other  collective  phenomena  which  are  observed  in  this  class  of  compounds  [2,  3]. 
Despite  considerable  effort  by  theorists  and  experimentalists  alike,  the  pairing 
mechanism  responsible  for  HTSC  in  these  materials  still  remains  a  mystery. 

High  temperature  superconductors  appear  to  be  attractive  candidates  for  the 
fabrication  of  optical  modulators  which  are  used  in  fiber  optic  communications. 
At  present,  fiber  optic  transmission  rates  are  limited  by  the  rate  at  which  optical 
data  can  be  modulated.  The  current  modulation  techniques  which  include  direct 
modulation,  electro-absorption  modulation,  and  Mach  Zender  modulation,  hmit 
the  transmission  rates  to  about  40  Obit  per  second.  Direct  and  electro-absorption 
modulation  are  hmited  by  the  carrier  mobility  of  gallium  arsenide  (GaAs).  The 
modulation  rate  of  Mach-Zender  modulators  is  limited  by  the  difficulties  in 
producing  rapid  high  voltage  pulses  required  to  switch  the  electro-optic  lithium 
niobiate  (LiNbOa)  material  used  in  these  modulators.  HTSC  can  easily  be  switched 
from  a  superconducting  state  to  a  resistive  state  by  passing  a  critical  amount  of 
current.  The  resistive  state  is  characterized  by  high  transmission  at  far  infrared 
frequencies,  as  compared  to  the  opaque  superconducting  state.  It  is  believed  that 

1 


2 

the  rate  of  switching  may  be  as  high  as  a  terabit  per  second.  This  high  rate  of 
switching  would  vastly  improve  the  transmission  rate  of  current  modulators  in  fiber 
optic  networks. 

Another  figure  of  merit  of  any  optical  modulator  is  its  extinction  ratio.  The 
extinction  ratio  for  a  superconducting  modulator  is  the  ratio  of  transmitted  light  in 
the  resistive  state  to  the  transmitted  light  in  the  superconducting  state.  It  will  be 
shown  that  a  superconducting  modulator  based  on  HTSC  films  offers  no  advantage 
over  the  current  modulators  which  are  available  commercially. 

This  dissertation  describes  a  semi-systematic  study  of  the  effects  of  trans- 
port current  on  the  far  infrared  transmittance  spectra  of  thin  superconducting 
YBa2Cu306+x  films.  I  have  confined  the  investigations  to  transmittance  measure- 
ments below  350  cm~^  There  are  two  main  reasons  for  this.  The  first  is  because 
the  effect  of  the  dephasing  and  subsequent  weakening  of  superconducting  order  is 
expected  to  be  seen  in  the  far  infrared  region.  Secondly,  the  changes  are  expected 
to  be  quite  small  and  optical  transmittance  is  far  more  sensitive  than  reflectivity  to 
such  small  changes. 

The  next  section  gives  a  brief  description  of  the  organization  of  the  remaining 
chapters  of  this  treatise. 

This  dissertation  consists  of  eight  chapters  including  the  introduction.  I  have 
tried  to  organize  the  chapters  in  a  logical  sequence,  so  that  all  the  concepts  needed 
are  easily  accessible. 

Chapter  2  describes  the  basic  optical  theory  of  the  propagation  of  light 
through  optical  media.  I  shall  discuss  the  various  techniques  which  are  used 
to  extract  the  optical  constants  of  a  material  from  measurements  of  optical 
reflectivity  and  or  transmittance.  Chapter  3  gives  a  discussion  of  the  basic  theory 
of  interferometric  spectroscopy  and  a  description  of  the  interferometer  used  in 
this  work.  A  review  of  our  current  understanding  of  the  electrodynamics  of  high 
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Tc  superconductors,  both  in  the  normal  and  superconducting  states,  is  given  in 
chapter  4,  where  various  models  used  to  analyze  optical  data  are  introduced. 

A  brief  description  of  the  film  preparation  techniques  and  their  characteri- 
zation by  a  standard  4-probe  method  is  presented  in  chapter  5.  Chapter  6  deals 
with  the  results  of  the  optical  reflectivity  and  transmittance  measurements  and  the 
extracted  optical  constants  for  the  thin  YBagCuaOe+x  films  used  in  this  work.  In 
chapter  7  I  present  the  results  of  electrooptic  measurements  in  the  superconducting 
state  of  the  same  films  used  in  chapter  6.  The  difficulties  encountered  in  interpret- 
ing the  data  are  discussed.  Finally,  chapter  8  concludes  this  dissertation  with  a 
review  of  the  results. 


CHAPTER  2 
OPTICAL  THEORY 

2.1  Introduction 

The  response  of  a  medium  to  an  external  oscillatory  electromagnetic  field 
is  quantified  by  its  complex  frequency-dependent  response  functions.  The  most 
commonly  used  response  functions  in  infrared  spectroscopy  are  the  complex 
conductivity  a{uj)  —  (Ji{io)  +  ia2{uj),  the  complex  dielectric  function  e{uj)  = 
ei{uj)  +  ie2{uj)  and  the  complex  index  of  refraction  N{uj)  —  n{u>)  +  iK{uj). 
Any  of  these  complex  functions  gives  a  complete  description  of  the  absorptive 
and  dispersive  nature  of  the  interaction  of  light  and  matter.  They  may  easily 
be  extracted  from  optical  reflectivity  and/or  transmittance  of  the  material  as 
a  function  of  the  energy  of  incident  light.  The  next  few  sections  of  this  chapter 
describe  how  these  functions  arise  in  a  natural  way  from  the  formal  theory  of 
the  interaction  of  light  with  matter.  I  shall  present  the  theory  as  applied  to  both 
isotropic  and  anisotropic  birefringent  media  where  the  functions  take  the  form  of 
tensors. 

2.2    Linear  Isotropic  Homogeneous  Conducting  Media 
This  section  describes  the  propagation  of  light  through  a  linear  isotropic 
homogenous  conducting  medium  [4,  5].  I  will  also  introduce  the  relevant  equations 
for  the  reflectivity  from  and  transmittance  across  a  boundary  formed  between  two 
such  semi-infinite  media. 
2.2.1    Wave  Propagation 

The  propagation  of  light  through  a  medium  is  based  on  two  fundamental 
but  equally  important  sets  of  equations.  The  first  more  fundamental  set  is  the 
macroscopic  Maxwell  equations.  The  vahdity  of  these  equations  is  based  on  the 
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assumption  that  the  wavelength  of  hght  far  exceeds  the  dimensions  of  the  unit 
cell  of  the  material  involved.  An  assumption  of  this  type  is  usually  satisfied  in 
most  cases  over  the  entire  range  of  interest  in  the  infrared  region  20  to  45,000 
wavenumbers  (20-45,000  cm~^).  In  Gaussian  units  the  Maxwell  equations  in 
differential  form  are  given  by 

VD  =  4npf  (2.1) 
V-B   =   0  (2.2) 

VXE  =  (2.3) 

c  ot 

VxH   =   ^^J;  +  i^  (2.4) 

c  c  ot 

where  E  and  D  are  the  electric  field  and  electric  displacement,  B  and  H  are  the 
magnetic  induction  and  magnetic  field,  c  is  the  speed  of  light  in  free  space  and  p; 
and  J f  are  the  free  charge  and  free  current  densitiy.  The  second  set  of  equations 
are  the  linear  constitutive  relations  between  the  field  vectors  given  by 

D  =  eiE       B  =  liH       Jf  =  GiE  (2.5) 

where  Ci  and  ai  are  the  frequency  dependent  dielectric  constant  and  conductivity  of 
the  material  respectively  while    is  the  complex  permeability  of  the  material  which 
describes  its  magnetic  properties.  These  relations  assume  that  the  electrodynamics 
of  the  medium  are  local.  By  using  the  linear  constitutive  relations  in  Maxwell's 
equations  and  ehminating  D  and  B  Maxwell's  equations  take  the  form: 

  47r 

Vf;   =   —pf  (2.6) 

^1 

V.H   =   0  (2.7) 
^     n  udH 

V  xE  =  2.8 

c  ot  ^  ^ 

^  47r    ^     ei  dE 


In  the  absence  of  free  charge  (2.6),  (2.8),  and  (2.9)  after  the  ehmination  of  H  give 
the  wave  equation  for  E, 

=  qi^^  +  l^f .  (2.10) 

d^t  at 

Consider  a  monochromatic  plane  wave  of  angular  frequency  uj  and  wavevector 
q  given  by 

^   \  =  \   ^°   UxpWo  ■  r  -  {2.11) 

where  Eo  and  Ho  are  complex  vector  amplitudes.  Upon  substitution  of  these 
expressions  into  Maxwell's  equations  we  obtain  the  following  results,  in  the  absence 
of  free  space  charges: 

q-E  =  0  (2.12) 
qH=0  (2.13) 

qxE   =    '^H  (2.14) 

c 


UJ  ^  4:TT 

iq  X  H   =   -i-E  H  

c  c 


UJ 


w  Atx 
a^+i—{l-ei)E=—[Jd  +  Js  +  Jp]  (2.15) 


47r 

where  Eqns.  (2.12),  (2.13)  and  (2.14)  indicate  that  the  vectors  E,  H  and  q 
form  an  orthogonal  triplet,  while  Jd,  Jf  and  Jp  represent  the  displacement,  free 
and  polarization  currents  respectively  in  Eqn.  (2.15).  It  follows  from  the  latter 
equation  that  the  properties  of  the  medium  can  easily  be  described  by  a  complex 
conductivity  a{(jo)  —  ai{uj)  +  ia2{oj)  where  0-2(0;)  =  a;(l  —  ei(a;))/47r  or  complex 
dielectric  function  e(a;)  =  ti{uj)  +  ^£2(0;)  where  e2{u})  —  4irai{ui)/uj.  It  follows  that 
Eqn.  (2.15)  may  be  rewritten  in  terms  of  these  complex  functions  as 

LU  UJ  Att 

iq  X  H  = -i-eE  = -i-E  +  — aE  (2.16) 

c  c  c 

After  combining  (2.16)  with  (2.14)  assuming  nonmagnetic  media  for  which  //  -  1 
we  obtain  the  relation     =  (^)^e.  It  is  easy  enough  to  define  a  complex  refractive 


index  N  which  leads  to  the  dispersion  relation: 

g  =  ^iV  =  -(n  +  m)  =  -x/i.  (2.17) 

c        c  c 

The  field  amplitudes  of  E  and  H  are  then  related  by 

H=ViE  =  NE.  (2.18) 

Whenever  a  plane  wave  is  incident  normally  on  a  conducting  medium  both  E  and 
H  are  attenuated  as     c'^'-e*^^"'-'^*).  This  is  a  damped  wave  with  a  skin  depth 
S  —  c/uK.  Since  the  power  dissipated  per  unit  volume  is  proportional  to  the 
square  of  the  field  amplitude,  the  corresponding  length  scale  for  power  loss  is  the 
absorption  coefficient  a  =  2ojk/c  =  2/ 5. 

Finally  the  power  flux  or  average  rate  of  flow  of  energy  per  unit  area  given  by 
the  Poynting  vector,  .     ,  , 

(S)  =  ^^{E  X  H*)  (2.19) 

OTT 

which  is  parallel  to  the  wave  vector  q  of  the  wave  for  an  isotropic  medium.  The 
physical  significance  of  the  real  part  of  the  optical  conductivity  oi  can  easily  be 
seen  from  the  relation  for  the  power  dissipated  per  unit  volume.  Using  Ohm's  law 
we  can  see  that  the  power  dissipated  per  unit  volume  is  given  by 

^  ^  \U{J  .  E*)  =  ^^{aE  .  E*}  =  \ay\E\\  (2.20) 

The  real  part  of  the  conductivity  is  thus  a  measure  of  the  absorptive  properties  of 
the  medium. 

2.2.2    Extraction  of  Optical  Constants 

In  many  cases  the  first  step  towards  the  extraction  of  the  optical  constants  of 
a  material  involves  measuring  the  optical  reflectivity  TZ{u))  of  a  thick  single  crystal. 
Here,  the  material  is  thick  enough  that  all  light  incident  upon  it  is  either  totally 
reflected  or  absorbed  with  little  or  no  transmission  (i.e.,  the  crystal  is  thicker  than 


the  penetration  depth  S).  The  difficulty  encountered  is  that  the  reflectivity  is 
desired  over  the  entire  frequency  from  a;  =  0  to     =  oo,  so  that  Kramers- Kronig 
analysis  can  be  carried  out.  This  section  concentrates  on  the  extraction  of  the 
frequency  dependent  optical  constants  from  the  optical  reflectivity  measured  over 
a  wide  frequency  range  on  thick  single  crystals  at  normal  incidence.  For  the  sake 
of  brevity  we  shall  not  consider  the  corresponding  formalism  which  deals  with 
light  incident  on  a  material  at  oblique  incidence,  since  normal  incidence  is  always 
assumed  in  all  cases  in  this  work. 

Consider  a  plane  wave  in  medium  1  with  a  complex  index  Ni  which  is  incident 
normally  on  the  boundary  between  itself  and  a  second  linear  isotropic  conducting 
medium  of  complex  index  iV2.  Let  {Ei,  Hi),{Er,  Hr)  and  {Et,Ht)  represent  the 
incident,  reflected  and  transmitted  electric  and  magnetic  field  amplitudes  at  the 
boundary  respectively —  which  are  all  parallel  to  the  surface.  The  boundary 
conditions  mandate  that  the  resultant  electric  and  magnetic  fields  at  the  boundary 
be  continuous  in  the  absence  of  surface  currents.  We  have  the  following  equations: 

Ej  +  Er      =  Ef 

(2.21) 

Hi  —  Hr    =  Ht- 

Prom  the  relation  between  E  and  H  given  in  equation  (2.18)  it  is  easy  to  write  the 
second  equation  in  (2.21)  in  terms  of  E  and  the  refractive  indices  using 

H^    =  N^Ei 

Hr  =  NiEr  (2.22) 
Ht   =  N^Ef 

This  allows  us  to  calculate  the  complex  amplitude  reflection  coefficients  for  the 
reflected  and  transmitted  flelds,  and  hence  the  fraction  of  the  incident  power  which 


is  reflected  and  transmitted.  The  amplitude  reflection  coefficients  are  given  by 

(2.23) 


4   =  =  r{u)e"»^^^) 


Ei         N1  +  N2' 

The  first  equation  in  (2.23)  is  written  in  exponential  form  in  terms  of  a  frequency 
dependent  amplitude  and  phase — the  significance  of  this  will  be  explained  later.  In 
the  laboratory  the  actual  measured  quantity  is  the  refiectivity  TZ{uj)  or  the  fraction 
of  the  incident  power  reflected  where  medium  1  is  free  space  with  iVi  =  1  and 
medium  2  has  index  N2  —  n  +  in.  The  observed  laboratory  reflectivity  is  given  by 

niu)  =  ff*  =  r(LV?  =  (1  ~"')  +  ^^  (2.24) 

It  is  possible  now  to  rewrite  the  first  equation  in  (2.23)  in  terms  of  the  reflectivity 
7?.(u;),  the  phase  change  upon  reflection  (j){uj),  and  the  complex  refractive  index  on 
reflection. Then  (2.23)  becomes 

V^e-*  =  f  =  (2.25) 

(l  +  n)  +  z«;  ^  ' 

where 

'-*=-r3^-  (2-26) 

The  reflectivity  'R{u))  and  phase         are  related  through  causality  by  the 
Kramers- Kronig  relations  [6,  7] 

00 

'^('^)  =  -  /   -i  da;'  2.27 

0 

or 

00 

,      1  Lu'-uj  d\nTZ{u;')  ,  ,  ,  , 

0 
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The  final  result  is  that  once  7^(u;)  and  (J){u)  are  known  its  is  easy  enough  to  invert 
eqn.  (2.25)  to  obtain  n{u)  and  k.{uj)  which  gives 

n(u;)  =  LZ^^^L   (2.29) 

1 +  n{u)- 2 y/R(u)  cos  (t){uj) 


,(u)  =  2v/^sin0(a;)  

We  may  also  define  the  transmittance  across  a  boundary  between  two  semi- 
infinite  media  of  complex  indices  Ni  and  N2-  The  fraction  of  the  incident  power 
transmitted  across  the  boundary,  for  light  incident  from  a  medium  1  with  complex 
index  Ni  and  passing  into  a  medium  N2  is  the  transmittance 

r(.).S|P|.4».  (2.31) 

Eqns.  (2.17)  provide  useful  relationships  between  the  different  optical  con- 
stants 

2  2 
ci    =   n'^  —  K'^, 

62   =  2nK, 

o"i    =  u£2/4:n, 

'  (2.32) 

(72   =  u{\  -  ei)/47r, 
a   =  2u!K,/c. 

We  have  discussed  the  use  of  Kramers-Kronig  analysis  to  extract  the  optical 
constants  of  a  thick  single  crystal  provided  the  optical  reflectivity  is  known 
over  a  broad  firequency  range.  The  reflectivity  as  measured  in  the  laboratory  is 
usually  only  known  over  a  finite  frequency  range,  and  various  extrapolations  of 
the  reflectivity  are  usually  employed  at  the  high  and  low  frequency  ends  of  the 
reflectivity  spectra.  These  extrapolations  usually  introduce  errors  in  the  optical 
constants  extracted  using  this  technique.  A  further  difficulty  involves  the  fact  that 
thick  single  crystals  are  not  always  available,  and  the  optical  reflectivity  under  such 
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circumstances  cannot  be  used  along  with  the  Kramers-Kronig  analysis  to  extract 
the  needed  optical  constants.  The  next  section  briefly  deals  with  the  theory  of 
stratified  media  which  is  more  relevant  to  the  current  work  which  deals  with  thin 
films  on  transparent  substrates. 

2.3    Stratified  Media 

A  stratified  medium  is  any  medium  whose  optical  properties  are  constant 
throughout  each  plane  perpendicular  to  a  particular  direction  [8] .  A  common 
example  of  stratified  media  is  multilayers  which  are  formed  from  a  succession  of 
thin  plane  parallel  films.  The  analysis  of  the  optical  characteristics  of  stratified 
media  is  very  different  from  that  of  thick  single  crystals  which  was  discussed 
previously.  In  this  section  the  basic  theory  of  the  optical  response  of  stratified 
media  will  be  presented. 

Consider  a  single  layer  of  thickness  d  and  complex  refiractive  index  N2  nested 
between  two  semi-infinite  media  of  indices  A'^i  and  N^,  (Figure  2-1).  Light  is 
incident  normally  from  the  left.  The  ampUtudes  of  the  incident,  reflected  and 
transmitted  electric  fields  are  Ei,  E[  and  E3  respectively.  The  corresponding 
magnetic  fields  (not  shown  in  the  diagram)  point  in  such  a  way  that  the  Poynting 
vectors  of  the  various  waves  point  in  the  directions  of  their  wavevectors  q  in  the 
diagram.  The  fields  E2  and  E'2  are  the  resultant  electric  fields  of  the  forward  and 
backward  travelling  waves  in  the  medium. 

The  electric  and  magnetic  vectors  are  parallel  to  both  interfaces  of  the  layer. 
The  boundary  conditions  require  the  resultant  tangential  electric  and  magnetic 
fields  at  each  interface  to  be  continuous.  At  the  left  interface  the  boundary 
conditions  give 

Electric  :  E^  +  E\    =   Eo  +  EL 

(2.33) 

Magnetic  :   N^Ei  -  N^E[    =   N2E2  -  iVa^;^. 
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Figure  2-1:  The  resultant  electric  fields  and  their  associated  wavevectors  of  a  plane 
wave  incident  normally  on  a  single  layer  of  thickness  d. 


(2.34) 


Similarly  at  the  right  interface  the  boundary  conditions  give 

Electric:  E2e''i^'^  +  E'.^e-''^^'^   =  E^ 

Magnetic:   N2E2e''i^'^  -  N2E'^e-''i-"^   =  N3E3. 

Here  the  complex  wave  vector  q2  =  {2ir/Xo)N2]  where  Aq  is  the  free  space  wave- 
length. In  order  to  simphfy  the  notation  I  shall  denote  the  complex  phase  shift 
by      After  eliminating  the  amplitudes  E2  and  E2  we  can  rewrite  the  remaining 
two  equations  in  matrix  form  as 

z 


1 

1 

+ 

-iVi 

El 

cos  62   =—  Bill  02 

—iN2  sin  82  cos  ^2 

Eqn.  (2.34)  maybe  rewritten  in  a  more  compact  form  as 


1 

N3 


Er, 


1 

1 

1 

+ 

r  =  M2 

-Ni 

t. 


(2.36) 


The  reflection  coefficient  is 


r  = 


E[ 
El 


(2.37) 
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and  the  transmission  coefficient  is  given  by 


t  = 


El 
El 


(2.38) 


The  matrix  M2  is  known  as  the  transfer  matrix 


-:r-  sin  ^2 

N2 


(2.39) 


cos52 

—iN2  sin  82        cos  82 
The  elements  of  this  matrix  are  functions  of  the  complex  refractive  index  of  the 
layer  N2  and  the  layer  thickness  d. 

It  is  easy  enough  to  generalize  this  formulation  to  a  case  where  there  are  N 
layers  numbered  2,3,4,...A''  +  1  with  complex  indices  N2,  N3,  N4, . . .  N^+i  and 
thicknesses  022,  c?3,  c^4,  •  ■  •  ^^at+i-  This  forms  a  multilayer  whose  transfer  matrix  is  the 
product  of  the  transfer  matrices  of  the  individual  layers.  By  applying  the  boundary 
conditions  at  each  interface  it  can  be  shown  that  the  multilayer  satisfies  a  similar 
matrix  equation  which  becomes 


t;  (2.40) 


where  the  MjS  are  the  various  transfer  matrices  of  the  layers  and  M  is  their 
product.  The  matrix  M  completely  characterizes  the  optical  response  of  the 
multilayer  and  gives  an  expression  for  the  amplitude  reflection  and  transmission 
coefficients  of  a  multilayer.  If  we  denote  the  elements  of  the  transfer  matrix  M  of 
the  multilayer  by  A,  B,  C,  and  D  respectively  such  that 


1 

1 

1 

1 

+ 

r  =  M2M3M4  ■  •  •  Mjv+i 

t=  M 

Ni 

Nn+2 

Nn+2 

M2M3M4  •  •  •  Mjv+1  =  M  = 


A  B 
C  D 


(2.41) 
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then  the  reflection  and  transmission  coefficients  become 

^  NiA  +  NN+2N1B  -  C  -  NN+2D 
NiA  +  NN+2N1B  +  C  +  NN+2D 

 2Ni  

~  NiA  +  NN+2N1B  +  C  +  NN+2D' 

The  reflectivity  K{u)  and  transmittance  T{uj)  then  become 


(2.42) 
(2.43) 


and 


rr*  = 

|r|2 

Nn+2 

tt*  = 

Nn+2 

If 


(2.44) 


(2.45) 


These  equations  are  quite  comphcated  even  for  a  single  layer.  If  the  optical 
constants  of  A''  —  1  layers  are  accurately  known  for  an  N  layer  system,  it  is  not 
possible  to  analytically  invert  these  equations  if  we  know  the  reflectivity  and 
transmittance  over  the  same  frequency  range.  It  is  possible  however  to  invert  the 
reflectivity  and  transmittance  of  a  single  layer  using  the  technique  outlined  in  [9]. 

2.4    Anisotropic  Media 
The  first  part  of  this  chapter  dealt  with  hnear  isotropic  homogeneous  media; 
where,  Maxwell's  equations  combined  with  the  linear  constitutive  relations  allowed 
wave  solutions  in  the  media.  The  media  were  assumed  to  be  isotropic  so  that  the 
optical  excitations  were  independent  of  the  direction  of  the  electric  field  in  the 
crystal  and  the  vectors  D  and  E  were  parallel.  We  now  consider  the  case  where 
the  medium  is  electrically  anisotropic  so  that  D  and  E  are  not  necessarily  parallel 
[10]. 

We  begin  with  Maxwell's  equations  for  a  conducting  medium: 


^      „     47r  ^  IdD 


V  X  E  =  , 

c  dt  ' 


(2.46) 
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and  consider  the  propagation  of  a  monochromatic  damped  plane  wave  of  frequency 
u  and  wave  vector  q  =  uj{N/c)s.  Here  N  represents  the  complex  refractive  index 
of  the  medium,  and  s  represents  the  unit  vector  in  the  wave  normal  dhection.  The 
complex  field  vectors  E,  D,  H,  B,  and  J  vary  as, 


exp 


lui  I  — r  •  s  —  t 


Substituting  this  result  into  Eqn.  (2.45)  gives 


Air 


Nsx  H  =  -D  +  —Jf,       NsxE  =  B 


(2.47) 


We  can  easily  ehminate  H  from  Eqns.  (2.46)  using  the  constitutive  equation 
B  =  fiH,  so  that  we  get 


47rz  , 
D  +  — Jf 


=  N^[E  -  s{s  ■  E)] 


(2.48) 


For  an  arbitrary  choice  of  orthogonal  cartesian  axes  in  an  anisotropic  crystal, 
the  components  of  the  electric  displacement  vector  D  and  the  electric  field  vector 
E  are  connected  by  a  dielectric  tensor     ;  while  the  components  of  the  free  current 
density  J/  are  related  to  the  E  field  components  by  a  conductivity  tensor  Uij.  In 
the  linear  regime  the  constitutive  relations  become 


Dy  6yxEx    "h   EyyEy    "I"  Ey^E^i 

Dz     =     ^zxEx  +  ^zyEy  +  ^zzEz, 


(2.49) 


and  -.-        •  ;^  '•  -  -  <  • 

-  Jx       ~      (^XxEx  +  (JxyEy  +  (TxzEz, 

^  \       Jy   =   CTyxEx  +  (TyyEy  +  ayzEz,  /  (2.50) 

Jz       =      <^ZxEx  +  CTzyEy  +  (JzzEz- 

The  nine  quantities  Exx,  Cxy,  •  •  •  are  constants  of  the  medium  and  constitute  the 
dielectric  tensor  while  the  constants  axx,  (^xy,  ■■■  form  the  conductivity  tensor. 
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These  relations  take  a  simple  form  if  we  use  principal  axes  where  the  tensors 
are  diagonal.  The  principal  axes  of  both  tensors  coincide  for  crystals  of  higher 
symmetry  classes  (of  at  least  orthorhombic  symmetry).  The  constitutive  equations 
then  take  the  simple  form 


Dx  =  €xEx,    Dy  —  ^yEy,    Dz  —  ezEz, 


(2.51) 


After  substituting  eqn.  (2.50)  into  eqn.  (2.47)  and  introducing  a  set  of  complex 
dielectric  constants 


(2.47)  gives 


47rz 

^3  =     — ^3      0  =      ^)  (2-52) 


N!^A^_  (2.53) 


Equation  (2.52)  represents  a  set  of  Unear  homogeneous  equations  in  Ej  which  have 
a  non  trivial  solution  if  the  complex  refractive  index  iV,  the  components  of  the 
complex  dielectric  constant  e^,  and  the  components  of  the  vector  s  satisfy  Presnel's 
equation, 

o2  o2  „2 


_^_^  +  3-^  +  3-^  =  0.  (2.54) 
N"^     M^x      iV2  A^e^ 
Let  us  now  consider  the  case  of  a  uniaxial  crystal.  Such  a  crystal  has  an 

optical  axis  which  lies  along  a  direction  of  high  symmetry.  The  optical  axis  is  one 

of  the  principal  axes  of  the  crystal.  Since  the  crystal  is  uniaxial  all  other  orthogonal 

axes  which  lie  in  a  plane  normal  to  it  are  optically  equivalent.  A  transverse  wave 

whose  electric  field  is  polarized  parallel  to  this  axis  travels  as  an  e-wave.  Similarly 

a  transverse  wave  whose  electric  field  is  polarized  perpendicular  to  the  optic  axis 

travels  as  an  o-wave.  Let  us  choose  the  z-axis  of  our  cartesian  coordinate  system  to 
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lie  along  the  optic  axis,  then  for  a  uniaxial  crystal  we  have 

Cz   —   Ce   =  (e-wave) 

^         ^  (2.55) 

e^  =  ey   =   Eo   =   No  (o-wave) 

where  eo  and  No  are  the  dielectric  function  and  refractive  index  corresponding  to 
the  o-wave  and  ie  and  Ng  are  the  corresponding  functions  for  the  e-wave. 

The  theory  of  the  reflectivity  from  and  transmission  through  a  uniaxial  layer 
has  already  been  developed  by  [11,  12].  Consider  light  of  frequency  u  which  is 
incident  from  a  medium  1  with  refractive  index  A^i  on  a  thin  unia^cial  layer  of 
thickness  Az  and  indices  No  —  \fTo  and  N^.  =  ^/e^  at  an  angle  of  incidence  0i  (see 
Fig.  2-2).  It  finally  passes  into  a  third  medium  of  complex  index  N2  at  an  angle  62. 
The  xy-plane  is  the  reflecting  surface  and  (a,  p,  7)  are  the  direction  cosines  of  the 
optic  axis  which  satisfy  the  relations:      -I-     -t-  7^  =  1  and  —1  <  a,  /9, 7  <  1.  The 
incident,  reflected  and  transmitted  electric  fields  have  a  z  dependence  given  by 

s  —  polarization 

incident  :  {0,1, 0)e''''^ 

reflected  :  (r^p  cos  ^i,  7'ss,  r^p  sin ^i)e~''i^ 

transmitted  :  (tgp  cos  ^2,  ^ssi —^sp  sin  ^2)6^^^'^"^^', 


(2.56) 


(2.57) 


p  —  polarization 

incident   :    (cos^i,  0,  -  sin^i)e''>^, 
reflected   :    {rpp  cos  9i,rps,rpp  sin  9i)e~^'^^^ , 
transmitted   :    (tpp  cos ^2,  ^ps, -^pp  sin ^2)6^''^^^"^^^ 
The  wave  vectors  used  in  these  equations  are  gi,2  =  {u>/c)Ni^2  cos^i  2- 

The  form  of  fields  in  the  layer  can  be  found  in  [11].  Let  us  consider  the  special 
case  of  normal  incidence  where     =  0  =  62  and  where  the  surrounding  media  are 
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Figure  2-2:  A  slab  of  anisotropic  crystal  with  parallel  faces.  The  laboratory  axes 
are  labelled,  x,  y,  z,  with  the  upper  and  lower  faces  parallel  to  the  xy-plane.  Light 
is  incident  on  the  upper  plane,  2  =  0  and  exits  at  z  =  Az.  Principal  ax;es  of  the 
crystal  are  denoted  by  a,  b,  and  c;  only  the  c-axis  is  shown  with  its  orientation  in 
the  lab  frame  given  by  its  direction  cosines.  The  plane  of  incidence  is  the  xz-plane. 


the  same  iVi  =  iV2.  In  this  case  the  reflection  and  transmission  coefficients  become 


fpp  — 


a2  +  /?2  ' 


(re  -  To), 


(2.58) 


^  aHe  +  PHo 

aP 


:{te  -  to). 


(2.59) 


The  reflection  and  transmission  coefficients  Tq,  r^,  tg  and  te  are  those  for  a  slab  of 
thickness  Az  and  are  given  by 

Soil  -  e^*^"^^) 


1  -  (s„)2e2i>^oAz' 
~  1  -  (sg)2e2i;CeAz  ' 

_  (l  +  g,)(l-g,)e'-^°^^ 
1  -  (soye^'^oAz  ' 

^  a  +  5e)(l_-Se)e^ 
1  -  (Se)2e'^=^^  ' 


(2.60) 


(2.61) 
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where  Ko  =  kNo,  Ke  =  kNgNe/Nj  where  the  single  bounce  reflectivities  are 

N,-No 
So  =  —  

iVi  +  iVo 

with 


(2.62) 


=     =  to  +  i\le  -  eo).  (2.63) 

These  equations  allow  us  to  calculate  the  reflectances  7^s,p)  ^-iid  transmittances  T^^p 
of  s  and  p-polarized  light  where 


Ir  1^4- Ir  P 

\'  ss\     1   \'  ap\  1 

1 ^pp 1      f    1 ^ps 1  ) 

l^pp 

1    "^^  l^psl 

(2.64) 

In  the  case  where  we  have  unpolarized  light  at  normal  incidence  the  reflectance  and 
transmittances  become  averages, 

7^=i(7^,  +  7^p)  =  ^(|r„p  +  |^e|2), 

T  =i(T,  +  Tp)  =  ^(|t„p  +  |i,|2).  (2.65) 

The  generalization  to  a  multilayered  sample  consisting  of  a  single  anisotropic 
layer  and  several  isotropic  layers  follows  directly.  It  can  be  seen  from  Eqns. 
(2.64)  that  the  total  reflectivity  and  transmittance  is  the  average  of  the  reflectiv- 
ity/transmittance  of  two  slabs  with  reflection  and  transmission  coefficients  (ro,  rg) 
and  {to,te).  In  a  multilayerd  system  the  matrix  theory  would  be  modified,  there 
would  be  two  sets  of  matrix  equations,  one  for  a  material  of  index      with  a  num- 
ber of  isotropic  layers  and  the  other  set  for  a  material  of  index  NoNe/N^  and  the 
same  set  of  isotropic  layers.  The  resulting  reflectivity  and  transmittance  would  then 
be  the  average  of  the  reflectivities  and  transmittances  of  both  sets  of  multilayers.  . 


20 

2.5    Drude-Lorentz  Model 

The  extraction  of  the  optical  constants  of  the  stratified  media  composing  a 
multilayer  usually  involves  assuming  a  particular  model  of  the  dielectric  response, 
and  fitting  the  reflectivity  or  transmittance  of  the  multilayer  by  varying  the 
parameters  of  the  model.  The  most  common  model  of  the  dielectric  response  used 
in  the  infrared  region  is  the  Drude-Lorentz  model.  This  model  although  ad-hoc  has 
proven  to  be  useful  in  describing  the  optical  characteristics  of  many  materials. 

The  basic  model  assumes  that  the  dielectric  function  can  be  written  as  a  sum 
of  three  terms  as 

e  =  eD  +  h  +  Coo-  (2.66) 

The  Lorentz  term  ii  represents  the  response  of  bound  atomic  electrons  which 
are  bound  to  their  cores  by  harmonic  forces  in  the  presence  of  viscous  damping 
which  represents  the  dissipative  losses  of  the  medium.  This  term  is  used  to  describe 
the  bound-carrier  interband  transitions,  phonons,  magnons,  and  many  other 
infrared  active  processes.  The  Lorentz  term  is  given  by 

h  =  J2   2     f   ■      >  (2-67) 
^  LUj  —  uj^  —  I'yjU 

where  uj,  ujpj,  and  Tj  are  the  center  frequency,  plasma  frequency  and  damping 
constant  respectively  of  the      Lorentzian.  The  plasma  frequency  Upj  is  related 
to  the  number  density  rij,  and  the  effective  band  mass  of  bound  carriers  m*,  by 

=  A-KUje^/m*.  The  plasma  frequency  is  sometimes  written  as  Upj  =  y/S'jUj 
where  Sj  is  an  oscillator  strength  which  represents  the  contribution  of  the 
oscillator  to  the  static  dielectric  constant. 

The  first  term  in  (2.65)  is  the  Drude  term.  It  describes  the  response  of  carriers 
which  are  essentially  free.  Typically  this  is  the  optical  response  of  a  metal.  The 
Drude  term  describes  intraband  processes  where  the  free  carriers  scatter  from 
impurities,  phonons,  or  scattering  processes  between  the  carriers  themselves.  The 
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Drude  term  has  the  form 

=  -  VtVt'  (2.68) 
where        is  the  Drude  plasma  frequency  defined  by  uJ^j)  =  4nne^/m*  where  n 
is  the  density  of  free  carriers,  m*  is  the  effective  mass  of  the  free  carriers  and  r  is 
the  relaxation  time  of  the  carriers  associated  with  scattering.  The  real  paxt  of  the 
Drude  conductivity  is 

47r  1  +  U^T^        1  +  U^T^ 

in  the  d.c  limit  the  conductivity  reduces  to 

=         =  (2.70) 

An        m*  ^  ^ 

2.6    Sum  Rules 

One  of  the  most  important  tools  used  in  the  analysis  of  optical  data  are  the 
sum  rules.  The  sum  rule  states  that  the  total  area  under  the  real  paxt  of  the 
conductivity  is  conserved.  This  is  simply  a  statement  of  charge  conservation  which 
becomes 

°° 

0 

In  many  cases  it  is  not  feasible  to  integrate  over  the  entire  frequency  range  and  it  is 
often  carried  out  up  to  some  finite  upper  frequency  hmit  Uc-  This  integral  gives  the 
effective  number  of  carriers  N^ff  per  atom  which  participate  in  the  optical  process 
up  to  that  frequency  and  is  usually  referred  to  as  a  partial  sum  rule.  This  sum  rule 
is 

/7r  Tl  c 

0 

where  m*  is  the  effective  mass  and  Ua  is  the  density  of  formula  units. 
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Sum  rules  are  very  important  in  understanding  the  superconducting  state 
electrodynamics  of  the  high  temperature  superconductors  (HTSC).  The  supercon- 
ducting state  is  marked  by  a  depression  of  the  low  frequency  conductivity  below 
the  normal  state  conductivity  above  Tc.  The  missing  spectral  weight  appears  as  a 
delta  function  at  the  origin  which  represents  the  response  of  the  condensate.  This 
effect  is  expressed  by  the  Ferrel-Glower-Tinkham  (FGT)  sum  rule 


where  ps  represents  the  superfluid  density.  The  sum  rules  have  also  been  applied  to 
the  study  of  the  redistribution  of  spectral  weight  with  doping  in  the  cuprates.  This 
will  be  discussed  in  greater  detail  in  chapter  4. 


oo 


(2.73) 
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CHAPTER  3 
INTERFEROMETRIC  SPECTROSCOPY 

3.1  Introduction 

One  of  the  most  versatile  tools  available  to  any  spectroscopist  who  is  interested 
in  the  optical  characteristics  of  materials  in  the  solid  state  is  the  Fourier-transform 
interferometric  spectrometer.  The  purpose  of  any  spectroscopic  instrument  is 
to  extract  the  frequency  spectrum  of  the  light  field  [13-16].  This  process  can  be 
accomplished  in  different  ways  depending  on  the  technique  used.  In  a  grating 
instrument  the  spectrum  is  measured  by  scattering  light  from  a  periodic  grating. 
Each  frequency  component  is  scattered  in  a  different  direction  and  the  separation 
is  done  by  using  a  sht  to  isolate  the  desired  component.  The  spectrum  can  then  be 
determined  by  passing  each  frequency  component  in  turn  onto  a  suitable  detector. 
This  technique  is  quite  time  consuming  and  inefficient  because  of  the  necessity  of 
sequential  sampling.  A  far  more  efficient  technique  involves  the  use  of  a  two-beam 
Michelson-type  interferometer.  The  spectral  components  in  the  interferometer  are 
mapped  onto  a  set  of  sine  functions  of  the  retardation,  and  can  easily  be  extracted 
by  computing  the  Fourier  transform  of  the  detector  signal.  The  purpose  of  this 
chapter  is  to  present  the  basic  theory  of  Foiurier  spectroscopy,  and  also  to  give  a 
brief  description  of  the  interferometric  spectrometer  used  in  this  work —  the  Bruker 
113v. 

3.2  Basic  Theory 

The  concepts  which  are  relevant  to  the  understanding  of  the  modern  inter- 
ferometric spectrometer  will  be  presented  in  this  section.  I  shall  discuss  the  ideal 
spectrometer  as  well  as  some  of  the  problems  which  are  encountered  in  practice. 
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3.2.1    Interferogram  and  Spectrum  ■ 

Figure  3-1  shows  a  simple  Michelson  interferometer.  Light  emitted  by  a  source 
is  incident  on  a  beamsplitter  where  it  is  partly  transmitted  into  one  arm  and  partly 
reflected  into  the  other.  The  beams  travel  to  the  end  mirrors  and  are  reflected,  and 
return  to  the  beamsplitter  where  they  recombine.  The  combined  beams  travel  to  a 
detector.  The  light  from  a  broadband  source  can  be  represented  by  a  TEM  wave 
having  Fourier  components  E{u)e^'^"^,  here  |fc|  =  2'ku,  where  u  is  the  frequency  of 
light  in  units  of  cm~^  which  is  related  to  the  wavelength  A  in  cm  by     =  1/A.  If 

Moving  mirror  2  ^ 

\  2;/2 
T 


O  

Source 

Beam  Splitter  Fixed  mirror 


O  Detector 

Figure  3-1:  A  simple  Michelson  interferometer.  The  vertical  arm  labelled  2  is  2;/2 
longer  than  the  horizontal  one  labelled  1;  light  in  the  vertical  arm  has  an  extra 
path  difference  of  x. 

V  is  in  cm~^  then  \v\  is  in  cm.  As  the  wave  passes  through  the  interferometer  its 
phase  advances  and  is  modified  by  the  beamsplitter  which  has  amplitude  reflection 
and  transmission  coefficients  rfic"'^'-  and  iftc'''^' ;  and  the  end  mirrors  with  reflection 
coefficients  rie"*i  and  r2e"^2  The  resultant  field  at  the  detector  is  a  superposition 
of  the  fields  from  both  arms 

E  =  Ei+E2  (3.1) 
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where 

El  =  Es{u)he''t"rie''>"ne''^^e^'''''^  (3.2) 

and 

E2  =  Es{u)ne''t'^r2e''f'Hbe"^'e^''"'^e^'''''''  (3.3) 

where  R  is  the  common  optical  path  length  of  both  rays  and  x  is  the  additional 
path  difference  of  ray  2.  An  ideal  interferometer  has  identical  arms  with  (rie"^'  = 
rgc^'^^  =  —  1)  and  an  ideal  beamsplitter  for  which  Vb  =  tf,  =  l/V^.  The  signal  at  the 
detector  is 

00  CXD 

S'(x)  =  27r  j\E\^du  =  -K  J p{u)[l  +  cos{2ttux)]  di^  (3.4) 
0  0 

where  p{u)  =  \Es{i')\'^  is  the  spectral  power  density  of  the  field.  The  first  term 
in  the  integral  gives  the  total  power  in  the  field.  It  can  be  seen  that  in  the  ideal 
case  the  detector  signal  at  zero  path  difference  5(0)  is  twice  that  at  infinite 
path  difference  5(00)  since  the  second  term  in  the  integral  vanishes  due  to  rapid 
oscillations  of  the  cosine  function.  Since  we  are  not  interested  in  the  first  term, 
subtracting  out  the  first  term  and  rewriting  the  integral  in  complex  notation,  and 
extending  to  negative  frequencies  we  have 

00 

I{x)  =  S{x)  -  S{oo)  =  ^  y  p(i/)e2'^'''^di/.  (3.5) 

—00 

I{x)  is  called  the  interferogram  which  is  the  Fourier  transform  of  the  spectrum 
p{iy).  The  spectral  power  density  may  easily  be  recovered  by  taking  the  inverse 
transform, 

00 

|p(i/)  =  J  /(x)e-2--dx,  (3.6) 
—00 

so  that  I(x)  and  p(i/)  form  a  Fourier  transform  pair 

H^)  ^  \p{^)-  (3.7) 
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3.2.2    Truncation  and  Apodization 

In  reality  it  is  impossible  to  measure  an  interferogram  out  to  infinite  path 
differences.  The  interferogram  must  be  truncated  at  some  point  in  space.  This  is 
achieved  by  multiplying  the  interferogram  by  a  function  A{x)  which  serves  to  cut 
off  the  interferogram  at  some  point. 

The  function  A{x)  is  usually  chosen  to  be  a  nonzero  symmetric  function  on 
the  interval  —^L<x<  |L  or  a  one-sided  function  on  0  <  x  <  L.  The  truncated 
interferogram  becomes 

I'{x)  =  I{x)  ■  A{x)  ^  p{u)  ^  d{u);  (3.8) 

so  that  the  calculated  spectrum  is  the  true  spectrum  convolved  with  the  transform 
of  the  truncation  function.  The  convolution  integral  on  the  right  is  defined  as 


oo 

p{iy)*d{u)  =  j  p{v')h{u-u')dv'.  (3.9) 

— oo 

In  general  both  p{u)  and  d^u)  may  be  complex  since  I{x)  and  or  A{x)  may  not  be 
symmetric  functions  and  so  possess  complex  transforms.  If  the  width  of  A{x)  is  L 
then  the  with  of  its  transform  d{u)  is  of  the  order  1/L,  and  thus  the  resolution  of 
the  interferometer  is  limited  by  the  maximum  path  difference  of  the  interferogram. 
Consider  a  simple  case  where  the  truncation  function  is  a  window  function 

\x\  >  \L 

^  (3.10) 

\A  <  \L 

W{x,  L)  ^  w{u,  1/L)  =  Lsm{-KuL)/{-KuL). 

The  window  function  is  not  a  convenient  truncation  function  because  its  transform 
is  the  sine  function  which  has  very  large  side  lobes  on  either  side  of  the  principal 
maximum.  These  large  side  lobes  may  cause  intensity  errors  in  the  vicinity  of  very 
sharp  spectral  features.  In  order  to  combat  this  undesirable  situation,  the  window 


function  is  usually  multiplied  by  an  apodization  function,  whose  transform  when 
convolved  with  the  sine  function  serves  to  suppress  the  side  lobes  at  the  expense  of 
resolution. 

There  are  many  apodization  functions  which  are  used  in  practice.  The 
choice  of  apodization  function  intimately  depends  on  the  desired  resolution  and 
accuracy.  Table  3-1  lists  a  few  such  apodization  functions  and  their  corresponding 
transforms. 

Table  3-1:  Apodization  functions  and  their  corresponding  transforms. 

A{x)  a(j7)  Half  width 

Unapodized  1  Lsinc(y)  1.22/ L 

Triangular  1  -  (2|x|/L)  lLsmc^{y/2)  1.80/L 

Cosine  cos{nx/L)  |L[sinc(y  +  |)  +  sinc(y  —  |)]  1.58/L 

Happ-Genzel    a  +  b cos{27rx / L)    L  {sinc{y)){a  -  [by^ /{y"^  -  tt^)]}  1.82/L 

Quartic  [1  -  (2x/L)^]^  4(27r)Uj|(y)/j/i  1.92/L 

All  A{x)  =  0  for  |x|  >  \L\y  =  ttuL 

3.2.3  Sampling 

In  practice  the  interferogram  is  recorded  by  digital  computer  which  performs 
the  actual  transform.  The  interferogram  is  usually  sampled  and  the  spectrum  cal- 
culated at  a  large  number  of  discrete  equally  spaced  points.  The  integral  transform 
is  then  replaced  by  a  sum.  Mathematically,  the  sampling  of  the  interferogram  can 
be  expressed  as  the  product  of  a  Dirac  comb  and  the  actual  interferogram.  The 
comb  and  its  transform  are  given  by 

oo 

()(x,Ax)  =  ^  5{x-jAx)  ^  (1/Aa;)()(i/,1/Aa;).  (3.11) 

j  =  -oo 

The  interferogram  is  then  sampled  at  points  jAx;  the  sampled  interferogram 
becomes 

I{x)  ■  ()(x.  Ax)  ^  {1/Ax)p{u)  *  0{u,  1/Ax);  (3.12) 
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where 

oo 

p(z/)*()(j/,l/Ax)  =  ^  I{jAx){exp[-2mjiAx)u]}Ax.  (3.13) 

j  =  -oo 

The  periodicity  of  the  complex  exponential  gives  rise  to  a  number  of  interesting 
consequences.  It  can  easily  be  seen  from  equation  3.13  that  for  a  fixed  value  of  j  if 
u  is  increased  to    +  1/Ax  the  exponential  remains  unchanged.  The  sum  implies 
that  the  convolution  with  the  comb  generates  infinitely  many  periodic  replicas 
of  p{u)  having  a  period  1/Ax  =  2i>m',  which  is  known  as  aliasing  (folding).  This 
implies  that  if  the  maximum  frequency  in  the  spectrum  Um  is  less  than  half  the 
period  of  the  comb  the  periodic  spectra  or  aliases  won't  overlap  and  there  will 
be  no  error  associated  with  the  sampling;  this  statement  is  usually  referred  to  as 
Nyquist's  theorem. 

It  is  often  the  case  that  one  is  interested  in  a  spectral  range  which  is  narrower 
than  the  range  from  0  to  Vm]  for  instance  if  we  are  interested  in  a  spectral  band- 
width Vi  <  V  <  Vm  where  vi  ^  0.  It  is  then  possible  to  increase  the  sampling 
interval  Ax  of  the  interferogram  beyond  l/{2i/rn)  to  l/2{i/m  —  ^i)-  This  speeds  up 
the  time  required  to  scan  the  interferogram.  The  range  of  the  spectrum  can  easily 
be  reduced  by  optical  or  mathematical  filtering  [17].  An  interferogram  which  is 
oversampled  simply  wastes  computer  time;  while  undersamphng  leads  to  errors  in 
the  calculated  spectra  due  to  spectral  folding. 
3.2.4    Phase  Errors 

In  an  ideal  interferometer  the  measured  interferofram  is  symmetric  about 
the  zero  path  difference  (ZPD)  or  "white-light"  position.  It  therefore  follows  that 
the  transform  of  a  two-sided  interferogram  extending  from  —L/2  to  L/2  gives 
a  purely  real  spectrum.  An  asymmetric  interferogram  may  arise  in  a  number  of 
different  ways.  The  simplest  example  of  a  phase  error  occurs  when  the  arms  of 
the  interferometer  are  not  identical  so  that  0i  ^  02-  Differences  in  the  absolute 
reflectivity  do  not  affect  the  phase.  Other  possible  causes  of  phase  error  include 
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the  electrical  phase  shift  associated  with  the  detector  and  also  sampling  error.  The 
transform  of  an  asymmetric  interferogram  is  a  complex  function  of  frequency.  The 
main  purpose  of  any  phase  correction  procedure  is  to  recover  the  original  spectrum. 

It  can  easily  be  shown  that  in  the  presence  of  a  phase  error  ^(t/)  the  interfero- 
gram becomes 

oo 

I{x)  =  I  y  p(iy)e^2w^'^^+'^('^»  di/  (3.14) 

-00 

so  that 

oo 

p(j.)e''^M=  j  /(x)e-''^''^^  dx.  (3.15) 

-oo 

The  error  may  be  easily  corrected  by  multiplying  the  transform  of  the  inter- 
ferogram by  e"^*^^"^'  which  gives  the  correct  spectrum  p{u)-  This  multiphcative 
procedure  was  first  proposed  by  Mertz  [18].  Another  technique  involves  talcing  the 
magnitude  of  the  calculated  spectrum 

where  piiv)  =  piv)  cos0(i/)  and  P2{i^)  =  pii')  sincplu)  are  the  real  and  imaginary 
parts  of  the  spectrum  calculated  from  a  two-sided  interferogram  extending  from 
— L  to  L.  Other  techniques  of  phase  correction  have  been  proposed  by  other 
authors  [19,  20].  In  practical  terms  the  phase  correction  is  usually  carried  out  by 
measuring  a  short  two-sided  interferogram  which  extends  from  —Li/2  to  Li/2.  The 
complex  transform  of  this  low  resolution  interferogram  allows  the  phase  0(j^)  to  be 
determined. 
3.2.5  Noise 

One  of  the  great  advantages  of  Fourier  transform  spectroscopy  over  conven- 
tional grating  spectroscopy  is  the  advantage  in  the  signal-to-noise  ratio  which 
achieved  due  to  multiplexing.  This  fact  was  first  pointed  out  by  Fellgett  [21].  The 
reduction  of  noise  in  the  calculated  spectrum  is  always  a  welcomed  feature  in  any 
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spectroscopic  technique.  An  advantage  of  this  kind  is  not  always  realized,  as  it 
depends  strongly  on  the  type  of  noise  which  is  present  in  the  spectrum. 

Noise  in  transform  spectroscopy  may  arise  in  a  variety  of  ways.  The  radiation 
field  invariably  contributes  noise  and  so  does  the  detector  and  its  associated  elec- 
tronics. Noise  may  be  characterized  as  additive  or  multiplicative.  The  noise  which 
is  produced  by  a  detector  and  its  electronic  components  tend  to  add  to  the  interfer- 
ometric  signal;  provided  that  the  electronics  have  strictly  hnear  characteristics;  this 
is  called  additive  noise.  The  radiation  field  produces  multiplicative  noise.  There  are 
two  main  categories  of  multiplicative  noise  namely;  photon  noise  and  signal  noise 
[22,  23].  Photon  noise  is  associated  with  the  quantum  nature  of  the  light  field.  The 
number  of  photons  N  arriving  at  a  given  point  will  fluctuate  randomly  by  ±y/N 
and  therefore  gives  rise  to  photon  noise.  Signal  noise  is  produced  by  fluctuations 
in  the  source  output  or  else  by  scintillation  of  the  intervening  medium  between  the 
source  and  the  detector.  Either  type  of  noise  may  have  a  frequency  independent 
power  spectrum — white  noise — or  a  frequency  dependent  power  spectrum.  A  com- 
mon example  of  frequency  dependent  noise  is  amplifier  noise  (also  called  1//  noise) 
which  diverges  as  1//,  where  /  is  the  frequency. 

The  relationship  between  random  noise  on  the  interferogram  and  noise  on  the 
corresponding  spectrum  as  well  as  the  Fellgett  advantage  can  easily  be  derived. 
The  sampled  interferogram  may  be  viewed  as  a  sum  of  a  noiseless  interferogram 
and  a  noisy  interferogram  B{x) 

B{jAx)  =  Ej[Io  +  /(jAx)]"  (3.16) 

Here  it  is  assumed  that  the  noise  scales  as  the  nth  power  of  the  detector  signal, 
where  Iq  is  the  dc  component  which  was  neglected  earlier.  Ej  is  a  statistically 
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random  parameter  whose  ensemble  averages  satisfy  the  following  relations 

\  ^  ^  '  .  (3.17) 

{E,)   =  0 

The  dc  component  Iq  is  related  to  the  average  value  of  the  spectrum  by  Iq  = 
p[uk)NAi',  where  N  is  the  number  or  sampUng  points  and  Av  =  1/L  where  L  is 
the  length  of  the  interferogram.  The  ensemble  average  of  the  noise  on  the  spectrum 
becomes 

=  {\h\')  =  E'iAx)'  J2l^o  +  lijAx)^.  (3.18) 

3 

When     »  1, 

6^  ^  E''{AxfNll''  (3.19) 

so  that 

h  =  E\p{vk)TN'''^'''{Au).  (3.20) 

The  comparison  if  the  signal-to-noise  ratio  of  an  inteferometer  and  that  of  a  grating 
monochromator  of  equal  resolution  Au  can  easily  be  obtained.  The  signal-to-noise 
ratio  of  the  interferometer  is 

(SNR).  =  ^  =  ^M^AT-^  (3.21) 

and  that  for  a  grating  monochromator  is  given  by 

The  ratio  of  these  two  expressions  gives 

N^^"^     n  =  0       detector  noise 


(SNR), 


1  =  I       photon  noise  (3.23) 

N'^/"^   n  =  1       source  noise. 
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This  last  equation  clearly  indicates  that  the  Fellgett  advantage  is  obtained  only  for 
the  n  =  0  case  where  the  noise  is  the  thermal  detector  noise  (Johnson  noise).  One 
interesting  consequence  of  this  observation  is  that  at  higher  frequencies  in  the  near 
infrared  to  ultraviolet  regions  of  the  spectrum  where  phonon  noise  is  the  dominant 
source  of  noise  the  signal-to-noise  ratio  of  an  interferometer  offers  no  advantage 
over  the  corresponding  grating  monochromator.  The  Felleget  advantage  is  only 
obtained  in  the  infrared  region  where  the  dominant  noise  source  is  thermal  detector 
noise.  Another  way  of  seeing  the  Fellgett  advantage  is  in  terms  of  the  measurement 
time.  The  noise  parameter  E  is  directly  proportional  to  where  t  is  the  time 

taken  to  measure  a  single  data  point  on  the  interferometer  or  a  single  point  in  the 
spectrum  on  a  grating  monochromator.  Then  it  is  easy  to  see  that  by  substituting 
E  =  Kt~^/'^  into  equations  (3.23)  with  n  =  0  and  (3.22)  respectively  that  the 
signal-to-noise  ratio  in  terms  of  the  measurement  time  becomes 

(SNR),  =  J^r^f/,  (3.24) 


and 


where  T  —  Nt  \s  the  total  measurement  time  for  the  inteferometer.  This  shows 
that  the  signal-to-noise  ratio  of  the  monochromator  is  only  equal  to  that  of 
the  interferometer  if  the  time  taken  to  measure  a  single  spectral  point  with  the 
monochromator  is  equal  to  the  total  measurement  time  of  the  interferometer. 

3.3    Bruker  113v  Spectrometer 
The  far  infrared  measurements  which  will  be  presented  later,  were  performed 
using  the  Bruker  113v  fast-scan  Fourier  transform  interferometric  spectrometer. 
It  covers  a  frequency  range  from  about  20-5000  cm~^;  where  the  far  infrared 
range  lies  between  20-600  cm'^  Figure  3-2  is  a  schematic  diagram  of  the  Bruker 
113v.  The  spectrometer  consists  of  four  main  chambers:  the  source  chamber,  the 
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interferometer  chamber,  the  sample  chamber  and  the  detector  chamber.  During  the 
measurement  the  entire  spectrometer  is  evacuated  to  prevent  absorption  by  carbon 
dioxide  and  water  vapor. 

The  source  chamber  houses  two  different  lamps;  a  mercury  arc  lamp  which 
is  used  in  the  far  infrared  range  and  a  globar  source  for  the  mid  infrared  range. 
Light  from  these  sources  is  channelled  by  a  series  of  mirrors  through  an  automated 
aperture  to  the  main  interferometer. 

The  interferometer  chamber  consists  of  a  two-sided  scanning  mirror,  optical 
filter,  control  interferometer  and  beamsplitters.  There  are  5  beamsplitters  which 
together  cover  the  entire  spectral  range.  Three  of  them  are  different  thicknesses 
of  mylar  which  are  used  over  different  spectral  ranges  in  the  far  infrared  and  the 
fourth  is  a  composite  germanium-potassium  bromide  (Ge/KBr)  beamsplitter  which 
is  used  in  the  mid  infrared.  The  final  beamsplitter,  also  used  in  the  far  infrared, 
is  made  of  metal- mesh,  but  is  rarely  used  because  of  the  losses  associated  with 
it.  The  control  interferometer  serves  the  important  purpose  of  ensuring  that  the 
interferogram  is  properly  sampled.  It  is  two  interferometers  rolled  into  one.  As 
the  scanning  mirror  moves  the  control  interferometer  generates  two  interference 
patterns.  It  has  a  white  light  source  and  a  helium-neon  (He-Ne)  laser  source  each 
of  which  produces  its  own  interference  pattern.  The  laser  wavelength  is  accmately 
known  and  so  the  zero-crossings  of  its  interference  pattern  is  used  to  mark  when 
the  interferogram  is  sampled.  The  white  light  interference  pattern  is  sharply  peaked 
at  at  zero  path  difference  and  this  position  marks  where  the  sampling  begins  on 
each  scan  of  the  moving  mirror.  This  is  the  technique  which  is  used  to  digitize  the 
interferogram.  The  two-sided  mirror  introduces  twice  the  path  difference  of  the 
simple  Michelson  interferometer  discussed  earher  for  a  given  mirror  displacement. 
After  passing  through  the  main  interferometer  the  light  passes  into  the  sample 
chamber. 
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Transmittance  measurements  are  usually  made  on  samples  mounted  in  the 
front  of  the  sample  chamber  at  the  transmittance  focus  i  in  the  figure,  while 
reflectivity  measurements  are  usually  done  in  the  back  of  the  chamber  where  the 
sample  focus  lies  at  the  point  j.  The  sample  stages  which  are  used  for  each  type 
of  measurement  are  not  shown  in  the  figure.  Finally,  the  light  enters  the  detector 
chamber. 

Signal  detection  is  performed  by  a  liquid  helium-cooled  silicon  bolometer  which 
is  shown  in  fig.  3-3.  This  detector  which  is  used  in  the  far-infrared  spectral  range, 
consists  of  three  main  parts:  1)  the  actual  detector,  2)  a  liquid  helium  can  which 
is  in  contact  with  the  detector  and  which  serves  to  cool  it,  and  3)  a  preamplifier 
which  amphfies  the  signal.  The  detector  is  maintained  at  cold  operating  temper- 
ature by  evacuating  the  detector  chamber  and  filling  the  helium  can  with  liquid 
helium.  Certain  design  characteristics  ensure  that  the  detector  remains  cold  for 
long  periods;  among  these  are  the  outer  nitrogen  can  which  is  filled  with  nitrogen 
and  a  thermal  radiation  shield  which  surrounds  the  detector  chamber.  At  optimum 
performance  the  detector  remains  cold  for  nineteen  to  twenty  hours.  Light  which  is 
focused  on  the  detector  is  amplified  by  the  preamplifier  and  sent  on  to  an  anolog 
to  digital  converter  of  interferometer  for  digitization.  The  most  notable  feature  of 
a  bolometer  of  this  type  is  the  very  high  signal-to-noise  ratio  which  is  obtained  in 
the  far-infrared  region  compared  to  other  detectors  which  operate  over  a  similar 
frequency  range. 
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A  Source  Chamber  ^  Automatic  beamsplitter  changer 

a  Near-,  mid-  orfar-IR  sources  ^  jwo-side  movable  mirror 
b  Automated  Aperature         ^  Control  interferometer    C  Sample  Chamber 
B  interferometer  Chamber      g  Reference  laser  D  Detector  Chamber 

c  Optical  filter  h  Remote  control  alignment  mirror      k  Near,  mid  orfar- 

IR  detectors 


Figure  3-2:  Schematic  diagram  of  Bruker  113-v  Fourier  transform  spectrometer. 


Figure  3-3:  Helium  cooled  bolometer  detector.  All  dimensions  are  in  inches. 


CHAPTER  4 

INFRARED  PROPERTIES  OF  HIGH-Tc  SUPERCONDUCTORS 

4.1  Introduction 

The  issue  of  superconductivity  in  the  cuprates  is  by  no  means  a  simple  one. 
The  interpretation  of  the  results  of  various  experiments  on  the  cuprates  have  often 
been  the  subject  of  intense  debate;  and  have  spawned  several  theories  relating  to 
the  underlying  mechanism  of  superconductivity  in  these  materials.  In  spite  of  an 
unprecedented  effort,  the  mechanism  of  high-Tc  superconductivity  in  the  cuprates 
still  remains  a  mystery. 

Infrared  spectroscopy  still  remains  a  premier  experimental  tool  in  the  in- 
vestigations of  the  electrodynamics  of  the  high-Tc  superconductors.  Immense 
progress  in  the  growth  of  high  quality  single  crystals  and  thin  film  specimens  has 
led  to  a  closer  convergence  amongst  the  spectroscopic  results  of  different  groups. 
The  physics  of  highly  correlated  electronic  systems — of  which  the  cuprates  are 
a  member — are  reflected  in  the  optical  constants  at  fairly  low  frequencies  which 
can  be  probed  by  infrared  (IR),  optical,  terahertz  (THz)  and  microwave  spectro- 
scopies. A  combination  of  model  independent  analysis  of  infrared  data  based  on 
the  sum  rules,  combined  with  comprehensive  studies  of  many  different  families  of 
the  cuprates  by  IR-THz/optical  spectroscopies  have  served  to  highlight  universal 
trends  in  the  optical  response.  The  main  purpose  of  this  chapter  is  to  briefly  review 
the  generic  optical  features  of  the  cuprates  and  to  present  some  of  the  techniques 
commonly  used  to  ajialyze  optical  data. 
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4.2    Generic  Trends 

4.2.1  Normal  State  Phase  Diagram 

Superconductivity  in  the  cuprates  appears  when  the  parent  antiferromagnetic 
insulating  phase  is  doped.  This  introduces  a  moderate  density  of  charge  carriers. 
Doping  is  achieved  either  by  chemical  substitution  (eg.  La2Cu04  La2_xMxCu04 
where  M  =  Ba,  Sr,  Ca);  or  by  doping  with  oxygen  (La2Cu04  —*  La2Cu044.5). 

The  generic  trends  of  the  evolution  of  the  optical  response  of  the  cuprates  in 
the  normal  state  with  temperature  and  doping  can  be  inferred  from  the  normal 
state  phase  diagram  [24].  There  are  a  number  of  common  elements  which  are 
worth  mentioning.  It  is  found  that:  1)  the  undoped  insulating  compounds  all  show 
antiferromagnetism  with  a  Neel  temperature  which  decreases  with  doping;  and 
becomes  conducting  upon  further  doping;  2)  superconductivity  occurs  for  doping 
levels  beyond  a  critical  level;  3)  a  transition  temperature  which  first  increases  with 
doping  (underdoped  region),  reaches  a  maximum  at  optimal  doping  and  finally 
decreases  with  further  doping  (overdoped  region);  4)  in  the  underdoped  regime 
the  superconducting  state  is  preceded  by  a  pseudogap  whose  onset  temperature 
decreases  with  increasing  doping;  5)  a  fairly  narrow  temperature  region  above 
the  transition  temperature  %  which  is  associated  with  the  formation  of  a  bosonic 
resonance  and  the  onset  of  high  in-plane  coherence. 

These  features  are  all  generic  features  of  the  hole-doped  cuprates  except  La- 
214  where  the  pseudogap  is  quite  weak  [25]  and  the  Cu  spin  relaxation  shows  no 
sign  of  a  spin  gap  at  lower  temperature  [26,  27].  The  pseudogap  and  bosonic  mode 
are  still  in  dispute  in  the  electron-doped  cuprates  such  as  Nd2-xCexCu04. 

4.2.2  Crystal  Structure 

The  cuprates  possess  a  layered  structure  which  is  often  classified  as  a  per- 
ovskite  crystal  structure.  This  is  a  label  given  to  a  generic  class  of  crystals  whose 
structure  resembles  that  of  SrTiOa  (such  as  YBa2Cu307)  and  a  derived  one  K2NiF4 


(La2Cu04  structure)  [28].  The  salient  structural  feature  of  the  cuprates  is  the 
presence  of  n  copper  oxide  planes  CUO2  separated  by  by  the  so-called  charge 
reservoirs.  Superconductivity  is  thought  to  be  associated  with  the  CUO2  planes 
while  the  charge  reservoirs  play  a  secondary  role.  The  La2-xSrxCu04  (La-214)  and 
Nd2-xSrxCu04  (Nd-214)  materials  all  have  a  single  CUO2  plane  and  are  usually 
called  the  single-layered  compounds;  they  are  usually  characterized  by  the  lowest 
transition  temperatures  at  optimal  doping.  The  most  frequently  studied  cuprates 
YBazCusOe+x  (Y-123)  and  Bi2Sr2CaCu208+x  (Bi-2212)  have  n  =  2  CUO2  planes 
per  unit  cell  and  higher  TcS  than  the  single-layered  compounds.  Many  different 
families  of  the  cuprates  may  be  synthesized  with  n  =  1, 2, 3, 4  and  the  transition 
temperature  increases  with  the  number  of  planes  n  up  to  n  =  3.  The  transition 
temperature  saturates  or  even  decreases  for  higher  n. 

The  undoped  cuprates  are  usually  described  as  Mott-Hubbard  insulators.  The 
large  coulomb  repulsion  between  Cu  ions  splits  the  Cu  3d  band  into  an  upper  and 
lower  Hubbard  band  with  a  gap  U.  The  energy  U  is  the  cost  of  creating  a  doubly 
occupied  Cu  site  in  the  upper  Hubbard  band.  Another  important  energy  scale  is 
the  hopping  matrix  element  t.  If  the  spins  on  adjcent  sites  are  oppositely  oriented 
they  can  hop  from  one  site  to  the  next  effectively  lowering  the  energy  of  the  system 
by  an  amount  proportional  to  -t'^/U.  This  is  equivalent  to  an  exchange  interaction 
J  =  U'^/U  which  leads  to  antiferromagnetism.  The  oxygen  2p  band  Hes  inside  the 
Hubbard  gap  so  that  the  charge  transfer  gap  Eg  is  less  than  U.  The  cuprates  are 
thus  referred  to  as  charge  transfer  insulators  [29].  The  optical  conductivity  in  the 
Cu02  (ab  plane)  thus  shows  an  absorption  edge  which  lies  between  1-2  eV  owing 
to  excitations  between  the  O  2p  to  the  Cu  3dx2_y2  states. 
4.2.3    Doping  Effects 

The  evolution  of  the  low  frequency  electronic  state  of  the  cuprates  with  doping 
has  been  studied  by  several  groups  in  the  infrared  [30,  31].  The  main  feature  of 
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the  process  is  that  doping  leads  to  the  development  of  free-carrier  absorption  and 
a  very  broad  absorption  band  in  the  mid-infrared  region  below  the  charge  transfer 
gap.  These  two  bands  develop  at  the  expense  of  the  charge  transfer  excitations 
which  weaken  with  doping.  The  free  carrier  band  develops  more  rapidly  than  the 
mid  infrared  band  which  softens  with  doping.  In  many  cases  the  bands  merge  and 
can  hardly  be  distinguished  at  high  doping  levels. 

At  moderate  doping  levels,  the  free-carrier  spectral  weight  which  is  propor- 
tional to  n/m*  increases  linearly  with  doping,  [30,  32]  where  n  is  the  carrier  density 
and  m*  is  their  effective  mass.  Close  to  optimal  doping  further  doping  does  not 
increase  the  electronic  spectral  weight  and  the  ratio  n/m*  saturates.  This  is  not  a 
generic  feature  as  a  different  effect  has  been  reported  for  Bi-2212.  It  has  been  found 
that  the  spectral  weight  of  Bi-2212  continues  to  increase  in  the  overdoped  regime 
[33]. 

4.2.4    Optical  Anisotropy 

The  strongly  layered  crystal  structure  of  the  cuprates  give  rise  to  anisotropy 
in  the  observed  electronic  properties  in  the  infrared.  The  in-plane  response  as 
measured  by  optical  reflectivity  in  the  normal  state  along  the  Cu02  planes  displays 
conducting  behavior  while  the  c-axis  response  perpendicular  to  the  planes  almost 
universally  displays  insulating  behavior.  This  behavior  is  observed  in  all  families 
as  a  fairly  large  ab-plane  reflectivity  at  low  frequencies  (a;  <  100  cm~^)  and  the 
appearance  of  a  plasma  edge  in  the  reflectivity  in  the  near  infrared  to  visible 
frequency  range  {uj  ^  1  -  2  eV);  while  the  c-axis  reflectivities  are  lower  and  show 
no  such  plasma  edge  above  the  spin  gap  temperature.  The  ab-plane  phonon  modes 
observed  in  the  parent  compound  are  also  strongly  suppressed  with  increased 
doping  due  to  screening  by  the  free  carriers.  By  contrast  the  c-axis  reflectivity 
spectra  are  dominated  by  a  large  number  of  infrared  active  phonons  at  all  doping 
levels,  which  are  also  observed  in  the  parent  compound. 
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There  are  many  physical  quantities  of  interest  which  display  anisotropy  in 
the  cuprates,including  dc  resistivity  [34-36],  optical  conductivity,  and  penetration 
depths  in  the  superconducting  state.  The  degree  of  anisotropy  depends  strongly 
on  the  particular  family  being  studied.  The  existence  of  copper  oxide  chains 
in  YBCO  implies  that  the  ab-plane  electronic  and  optical  properties  are  also 
anisotropic  [37,  38];  ab-plane  anisotropy  has  also  been  observed  in  Bi-2212  [39- 
41].  These  features  have  been  studied  by  polarization  measurements  along  the 
a  and  b-axes  by  several  groups  above.  The  observed  anisotropy  is  not  confined 
to  the  normal  state;  the  superfluid  density  is  also  anisotropic  and  varies  from 
family  to  family.  In  most  cuprates  the  condensate  plasma  frequency  extracted 
from  the  ab-plane  optical  conductivity  lies  between  4000  and  10000  cm~^;  while 
the  corresponding  c-axis  plasma  frequencies  range  between  1  and  1000  cm~^  This 
gives  rise  to  strong  anisotropy  in  the  calculated  penetration  depths  extracted  from 
optical  measurements  (note:  the  penetration  depths  are  just  the  reciprocal  of  the 
condensate  plasma  frequency  in  cm~^  so  the  c-axis  penetration  depths  are  always 
larger  than  the  in-plane  values).  For  Y-123  there  is  superfluid  density  associated 
with  the  chains  so  that  the  b-axis  penetration  depths  A{,  are  smaller  than  the 
corresponding  a-axis  values  Aq.  Light  polarized  along  the  a-axis  probes  only  the 
Cu02  planes  while  light  polarized  along  the  b-axis  probe  both  the  chains  and 
planes. 

4.3    Models  of  Optical  Conductivity 
Over  the  last  decade  and  a  half  a  number  of  models  have  been  used  to  analyze 
the  in-plane  reflectivity  and  transmittance  data  of  the  cuprates.  Each  model 
has  its  advantages  and  its  disadvantages.  The  dynamical  conductivity  obtained 
from  Kramers  Kronig  analysis  when  fitted  to  any  of  these  models  allows  for 
the  extraction  of  various  optical  parameters  which  give  important  information 
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about  the  low  energy  excitations  and  the  electronic  state  of  the  cuprates.  A  brief 
description  of  the  commonly  used  models  will  be  given  in  this  section. 
4.3.1    Two-Component  Model 

One  of  the  earliest  models  used  to  fit  the  normal  state  optical  conductivity 
of  the  cuprates  is  the  two-component  model  [2,  40,  42].  The  model  arose  out  of 
the  observation  that  doping  led  to  the  transfer  of  spectral  weight  into  the  charge 
transfer  gap  of  the  parent  compounds.  At  light  to  moderate  doping  levels  sharp 
features  are  observed  in  the  mid-infrared  region  together  with  a  Drude-like  band 
at  low  frequencies  {u>  <  500  cm"^).  Two  resonances  were  observed  in  the  mid 
infrared;  a  very  broad  resonance  at  0.5  eV  which  softened  with  doping  together 
with  a  second  resonance  at  0.25  eV.  The  feature  at  0.25  eV  was  close  to  the 
anti- ferromagnetic  exchange  energy  J  [43].  The  mid  infrared  absorption  has  also 
been  attributed  to  polaron  modes  [44-46].  The  two  component  model  assumes  a 
dielectric  function  of  the  form 

2  2 
U)  r,  U  ■ 

'"^^^  =  +      iu;]  -  J-  ^^7,)        =  + W-.(u;)  +  6^    (T  >  T,). 

(4.1) 

The  first  term  in  eq.  (4.1)  represents  the  Drude  free-carrier  response  with  plasma 
frequency  Upr,  and  relaxation  rate  l/r.  The  sum  represents  additional  contributions 
to  the  dielectric  response  in  the  mid  infrared  which  are  modeled  with  Lorentzian 
oscillators  with  center  frequency  luj,  plasma  frequency  Upj  and  line  width  -jj. 
Finally  e^o  includes  the  high  frequency  contribution  from  processes  above  the  ■ 
measured  range.  A  two-component  model  is  always  questionable  at  higher  doping 
levels  of  the  cuprates,  where  the  Drude  and  midinfrared  bands  are  sometimes  not 
clearly  resolved. 

In  the  superconducting  state,  a  two-fiuid  model  is  used  where  a  fraction  fs  of 
the  free  carriers  condense  to  form  superfluid,  and  the  remaining  fraction  {1  -  f^) 
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represents  the  response  of  thermally  excited  quasiparticles  which  still  exhibit  a 
Drude  response.  The  two  fluid  conductivity  is  written  as 


2  2 


+  ^mid-ir{(^)-  (4.2) 


47r(l  —  iut)  in 

The  mid-infrared  conductivity  (Jmid-ir  usually  shows  little  or  no  temperature 
dependence.  The  superfluid  response  in  the  dc  limit  which  gives  zero  resistivity 
in  the  superconducting  state  is  represented  by  the  delta  function  6{u)  in  the 
conductivity. 

4.3.2    One-Component  Models 

The  extended  Drude  model.  In  many  cases  the  infrared  conductivity  of  the 
cuprates  below  the  charge  transfer  region  is  fitted  with  an  extended  Drude  model 
[47-49].  The  model  assumes  that  the  low  frequency  conductivity  is  due  to  a  single 
type  of  charge  carrier  with  a  frequency  dependent  scattering  rate.  The  conductivity 
can  be  extended  beyond  the  simple  Drude  form  by  maJcing  the  damping  term 
complex  and  frequency  dependent  such  that  1/r  =  M{uj)  =  M'{ijj)  +  iM"{u)),  where 
M{u))  is  a  complex  memory  function.  This  model  has  been  extensively  used  to 
study  the  carrier  dynamics  of  the  cuprates.  The  memory  function  maybe  rewritten 
as  M{u,T)  —  1/t{uj,T)  —  iuX{uj,T)  so  that  the  conductivity  becomes 

'     AnM{uj,T)-iu     47rl/T(u;,r)-za;[l  +  A(a;,r)]'         ^  ^ 

The  quantities  1/t{uj,T)  and  A(cj,T)  are  the  frequency  dependent  scattering  rate 
and  effective  mass  enhancement  due  to  interactions  of  the  carries  with  a  cloud  of 
excitations.  They  are  related  to  a{uj)  through  the  relations 

(  1 


i.e.,  the  scattering  rate  is  proportional  to  real  part  of  the  reciprocal  conductivity; 
while  the  mass  enhancement  factor        varies  as  the  imaginary  part  of  the  same 
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response  function 

l  +  X{u;)  =  -^-^(^).  (4.5) 

Marginal  Fermi  liquid.  Another  very  commonly  used  model  is  the  marginal 
Fermi  liquid  model  due  to  Varma  et  al  [50,  51].  The  model  assumes  that  the 
charge  carriers  interact  with  a  fairly  flat  spectrum  of  excitations  over  the  interval 
T  <  u>  <  Uc  where  ujc  is  a  high  frequency  cutoff.  The  dielectric  response  for  this 
model  can  be  written  as 

where  E  is  the  quasiparticle  self-energy,  given  by 

'iTT  —  iuj\     ,  2- 


S  -  2\uj  log    -  m'XT-  (4.7) 


here  T  is  the  temperature  A  is  a  coupling  constant  and  uJc  is  the  cutoff  frequency. 
The  limiting  forms  of  this  expression  go  as 

I    TT^AT,  UJ<T 

-^S(a;)  ~  <^  (4.8) 
I    ttAo;,    uj  >T. 

Since  the  imaginary  part  of  the  self-energy  is  effectively  the  scattering  rate,  eq. 
(4.8)  predicts  hnear  variation  of  the  dc  resistivity,  which  is  observed  in  most 
transport  studies  of  the  cuprates.  In  a  similar  fashion  the  real  part  of  the  self- 
energy  gives  the  mass  enhancement  of  the  carriers; 


(i)  • 


where  m*  is  the  frequency  dependent  renormalized  mass  and  rub  is  the  band  mass 
which  appears  in  the  plasma  frequency  Up  =  inne'^/mi,. 

Nested  Fermi  Hquid.  The  anomalous  optical  properties  of  the  normal  state  of 
the  cuprates  have  been  described  in  terms  of  Fermi  surface  nesting  by  Virosztek 
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and  Ruvalds  [52,  53].  This  nested  Fermi  liquid  approach  attributes  the  anoma- 
lous normal  state  properties  to  carriers  which  lie  on  parallel  regions  of  the  two 
dimensional  Fermi  surface  which  are  spanned  by  common  nesting  wavevectors. 
This  nesting  causes  quahtative  changes  in  the  susceptibility  and  self-energy  of  the 
carriers.  It  is  assumed  that  the  dominant  scattering  process  is  electron-electron  or 
hole-hole  scattering,  for  which  the  imaginary  part  of  the  complex  self-energy  (which 
is  the  scattering  rate)  becomes 


here  a  is  a  dimensionless  coupling  constant  and  /?  is  a  constant  of  the  order  of 
unity.  The  above  equation  predicts  a  similar  temperature  and  frequency  depen- 
dence of  the  scattering  rate  as  the  marginal  Fermi  Uquid  model.  Therefore  the 
scattering  rate  is  linear  in  the  higher  energy  scale  u)  or  T.  The  model  also  predicts 
a  strong  mass  renormalization  at  low  frequency, 


Superconductivity  in  the  cuprates  occurs  when  the  temperature  is  lowered 
below  the  transition  temperature  Tc-  The  superconducting  state  is  a  macroscopic 
quantum  state,  where  are  certain  fraction  of  the  carriers  condense  to  form  su- 
perfluid.  This  idea  was  first  proposed  by  London  and  London  in  connection  with 
conventional  superconductors  [54].  The  density  of  the  condensed  carriers  is 
related  to  the  magnetic  penetration  depth  as 


—  $5S(cj)  =  Q;max(/3r, 


(4.10) 


TT      [max(/?T,  a;) 


(4.11) 


where  Uc  is  a  cutoff  frequency. 


4.4    Superconducting  State  Electrodynamics 
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In  a  simple  London  superconductor  the  real  and  imaginary  parts  of  the  complex 
conductivity  are  ai{uj)  =  (n^eV"^*)! 5(0),  and  0-2(0;)  =  {nge^ /m*u)  where  the 
area  under  the  delta  function  is  proportional  to  the  density  of  superconducting 
electrons.  The  simple  London  theory  assumes  local  electrodynamics  where  the 
penetration  depth  \l  is  much  greater  than  the  superconducting  coherence  length  ^ 
[55].  This  condition  is  satisfied  in  the  cuprates  where  for  YBCO  at  optimal  doping 
and  a  Tc  of  the  order  of  93  K  the  ab-plane  penetration  depth  Aa6  ~  1500  A  while 
the  ab-plane  coherence  length  ^ab  ~  15  A  [28]. 

These  simple  expressions  for  the  conductivity  assume  complete  condensation 
of  all  the  Drude  charge  carriers  below  Tg.  It  predicts  a  vanishing  real  conductivity 
a  I  at  all  nonzero  frequencies — a  clearly  untenable  state  of  affairs.  The  two  fluid 
model  eq.  (4.2)  rectifies  this  by  assuming  that  below  Tc  the  total  Drude  oscillator 
strength  is  shared  between  the  condensate  and  the  thermally  excited  quasiparticles. 
Both  the  simple  London  theory  and  the  two  fluid  model  fail  to  describe  the  energy 
gap  2A  in  the  density  of  states  which  is  the  centerpiece  of  the  BCS  theory  of 
conventional  superconductors.  The  gap  in  conventional  superconductors  could 
be  observed  by  infrared  spectroscopy  in  the  dirty  limit  where  the  scattering  rate 
far  exceeds  the  gap,  l/r  ^  2 A.  At  low  temperatures  where  the  quasiparticle 
excitations  are  frozen  out  and  Oi  {uj)  =  0  for  a;  <  2 A  the  gap  appears  as  an 
absorption  edge  in  Ui  {u)  at  the  gap  frequency  uj  =  2 A. 

Many  of  the  early  infrared  studies  of  the  cuprates  in  the  superconducting 
state  assigned  a  gap  frequency  to  certain  features,  such  as  a  shoulder  observed  in 
the  reflectance,  or  to  the  minimum  in  the  extracted  conductivity  at  far  infrared 
frequencies.  These  assignments  have  since  been  discredited  by  the  observation 
that  the  cuprates  are  in  the  clean  hmit  normal  skin  effect  regime,  where  (l/r  <C 
2A;^o       I  <^  ^l)  ie.  the  scattering  rate  is  much  smaller  than  the  gap  and  the 
mean  free  path  /  lies  between  the  coherence  length  and  penetration  depth.  This 
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means  that  when  condensation  occurs  the  majority  of  the  spectral  weight  of  the 
conductivity  will  be  transferred  to  the  delta  function  at  the  origin  [42,  56].  The 
gap  feature  if  present  will  be  very  weak  as  there  is  little  spectral  weight  left  for 
excitations  above  the  gap.  i..^ 

4.5    Survey  of  Experimental  Results 
4.5.1    Doping  Dependence  of  Optical  Constants 

One  of  the  first  systematic  studies  of  the  effect  of  doping  on  the  in-plane  low 
frequency  carrier  response  was  due  to  Orenstein  et  al.  on  microtwinned  crystals  of 
YBa2Cu306+x  [30].  Later,  Uchida  et  al.  [31]  performed  similar  ab-plane  reflectivity 
measurements  at  room  temperature  on  La2_xSrxCu04  single  crystals.  The  results 
of  these  latter  measurements  and  the  extracted  optical  conductivity  a{ijj)  are  shown 
in  fig.  (4-1).  The  next  few  sections  highlight  some  of  the  characteristics  which  are 
common  to  all  doped  cuprates. 

Optical  reflectivity  .  The  reflectivity  shows  a  number  of  characteristic  features 
which  are  worth  mentioning.  Firstly,  the  undoped  phase,  x  =  0,  shows  a  reflectivity 
typical  of  an  insulator;  with  a  phonon  mode  near  0.09  eV  which  precedes  a  very 
broad  flat  reflectivity  at  higher  frequency.  The  structure  which  appears  in  the 
reflectivity,  peaked  at  2  eV  has  been  assigned  to  the  Cu-0  charge  transfer  band. 
Then,  with  doping  the  reflectivity  grows  in  the  midinfrared,  and  a  reflectivity  edge 
develops  at  the  expense  of  the  charge  transfer  band  which  weakens.  The  0.09  eV 
phonon  mode  is  also  masked  with  doping.  Finally  the  plasma  minimum  at  ^ 
1  eV  is  almost  insensitive  to  doping  and  in  fact  displays  a  slight  softening  in  the 
overdoped  regime  x  =  0.34.  The  plasma  minimum  corresponds  to  the  zero  crossing 
of  the  real  part  of  the  dielectric  function  ei(a;).  The  plasma  minimum  lies  between 
the  midinfrared  and  charge  transfer  bands;  the  midinfrared  band  makes  a  negative 
contribution  to  ei(a;)  while  the  charge  transfer  band  makes  a  positive  contribution. 
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The  softening  of  the  plasma  minimum  is  an  indication  of  a  decrease  of  the  spectral 
weight  associated  with  these  bands  in  the  strongly  overdoped  regime. 

Optical  conductivity.  The  corresponding  optical  conductivity  extracted  from 
the  reflectivity  is  shown  in  fig  (4-1).  The  insulating  character  of  the  undoped 
compound  is  clearly  seen  in  the  absorption  edge  which  begins  at  about  1  eV  and 
is  peaked  at  ^  2  eV  and  is  attributed  to  charge-transfer  processes.  At  light  doping 
the  charge-transfer  band  is  weakened  and  a  peak  centered  at  about  0.5  eV  in  the 
midinfrared  develops.  As  the  doping  increases  a  smaller  feature  at     1.5  eV  is  seen. 

The  peak  in  the  midinfrared  band  softens  and  strengthens  with  increased 
doping,  as  does  the  Drude,  peak  and  becomes  just  a  shoulder  at  a;  =  0.15  which 
corresponds  to  optimal  doping.  The  peak  merges  with  the  the  Drude  conductivity 
centered  at  a;  =  0  in  the  overdoped  regime  and  a  clear  separation  between  the  two 
bands  can  no  longer  be  made. 

Oscillator  strength  sum  rule.  The  observed  changes  in  the  conductivity  may  be 
further  quantified  by  considering  the  effective  electron  number  Neff{uj)  per  copper 
atom  contributing  to  the  optical  processes  up  to  a  frequency  lu.  This  is  given  by 
the  paxtial  sum  rule 
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here  nie  is  the  ft-ee  electron  mass  and  Vceii  is  the  unit  cell  volume.  Fig.  (4.2)  shows 
the  results  of  the  calculation  for  La2_xSrxCu04  from  Uchida  et  al.  [31]. 

In  undoped  samples  most  of  the  spectral  weight  is  accumulated  at  energies 
exceeding  the  charge-transfer  gap,  while  below  u  =  1.5      Neff{uj)  is  flat  and 
almost  zero.  As  doping  increases  the  spectral  weight  at  low  frequency  increases 
rapidly  as  states  appear  in  the  gap  region.  For  x  <  0.1,  Neff{uj  <  1.5  eV)  increases 
faster  than  x.  This  is  greater  than  what  is  expected  if  one  assumes  each  Sr  ion 
produces  a  single  hole  per  unit  ceU  with  effective  mass  equal  to  the  free  electron 


(4.13) 


mass.  This  has  been  observed  in  other  cuprates  [57].  It  is  evident  from  the  figure 
that  doping  affects  primarily  the  states  below  3  eV  and  that  spectral  weight  is 
transferred  from  the  charge-transfer  band  to  the  midinfrared  and  Drude  bands. 
At  doping  levels  up  to  0.2  eV  the  spectral  weight  converges  at  3  eV,  but  at  higher 
doping  there  is  an  overall  decrease.  This  is  in  accord  with  the  observation  that  the 
spectral  weight  in  the  high-Tc  materials  saturates  or  even  decreases  above  optimal 
doping  [33] .  The  reason  for  this  strange  behavior  is  still  a  mystery. 
4.5.2    Midinfrared  Band 

The  midinfrared  conductivity  has  been  extracted  from  optical  reflectivity  by 
several  authors.  Kamaras  et  al.  [42]  extracted  the  midinfrared  conductivity  of  Y- 
123  by  employing  a  two-component  fit  to  the  refiectivity  at  temperatures  above  Tc 
and  subtracting  the  Drude  component  from  the  total  conductivity,  and  compared  it 
with  the  total  conductivity  below  Tc.  A  similar  procedure  was  employed  by  Reedyk 
et  al.  to  Bi-2212  [58].  Tanner  et  al.  [59]  used  an  iterative  procedure  to  extract 
the  free  carrier  and  midinfrared  components  of  Bi-2212,  and  Gao  [9]  extracted  the 
midinfrared  band  of  La-214  in  a  similar  way. 

The  midinfrared  band  displays  a  number  of  common  features  which  are  worth 
mentioning:  1)  the  band  displays  an  edge  which  lies  in  the  far  infrared  between 
80-150  cm"^  below  which  the  conductivity  vanishes  2)  a  notch-like  structure 
usually  between  400-500  cm"^  whose  strength  varies  with  temperature  3)  a  broad 
maximum  somewhere  between  750  cm"^  to  1000  cm"^  which  shows  a  fairly  large 
enhancement  below  Tc  as  compared  to  the  weak  temperature  dependence  of  the 
band  for  u  >  1500  cm^^  The  notch-like  structure  has  been  attributed  to  a  phonon 
mode  resulting  fi-om  the  coupling  of  bound  carriers  in  the  midinfrared  to  phonons— 
the  charged  phonon  model.  It  has  been  found  that  the  temperature  dependence  of 
the  midinfrared  band  may  be  fitted  with  a  set  of  Lorentzians  by  fixing  their  center 
frequencies  and  oscillator  strengths  while  varying  the  fine  widths  hnearly  with 
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temperature  as 

7(r)  =  7(0)  +  ksT.  (4.14) 

This  behavior  has  been  investigated  by  Thomas  et  al.  [60]  from  optical  mea- 
surements on  the  cuprate  semiconductor  Nd2Cu04_x  where  a  similar  behavior  is 
observed.  This  linear  temperature  dependence  was  modelled  by  considering  the 
interactions  of  bound  electrons  with  bosons  in  the  antiferromagnetic  lattice.  A 
similar  temperature  variation  is  observed  in  the  Drude  relaxation  rate  below. 
4.5.3    Relaxation  Rate 

One  of  the  attractive  features  of  infrared  spectroscopy  is  the  ability  to  probe 
quasiparticle  dynamics  in  the  superconducting  state.  A  two-component  fit  in  the 
normal  state  combined  with  a  two-fluid  fit  in  the  superconducting  state  allows 
one  to  study  the  variation  of  the  carrier  relaxation  rate  with  temperature.  Such 
an  analysis  of  the  ab-plane  response  has  been  done  for  Y-123  [42,  30,  61],  La-214 
[9,  62],  and  Bi-2212  [41,  63].  In  all  cases  the  width  of  the  Drude  peak  narrows  with 
temperature  in  a  linear  fashion  as  expected  from  transport  measurements.  There 
is  httle  or  no  temperature  dependence  of  the  Drude  plasma  frequency,  so  that 
changes  in  the  observed  resistivity  can  be  attributed  entirely  to  a  depression  if  the 
relaxation  rate  with  decreasing  temperature. 

Liu  [41]  has  carried  out  systematic  studies  of  the  carrier  relaxation  rate  as  a 
function  of  doping.  The  results  of  such  studies  are  that  the  carrier  scattering  rate 
remains  linear  in  the  normal  state  for  underdoped  cuprates  even  at  temperatures 
below  the  observed  pseudogap  temperature  T*  where  characteristic  deviations  are 
observed  in  the  ab-plane  dc  resistivity. 

The  superconducting  state  is  characterized  by  a  suppression  of  the  quasipar- 
ticle scattering  rate  below  Tc.  This  is  a  general  trend.  It  is  widely  believed  that 
a  gap  is  formed  in  the  excitation  spectrum  responsible  for  scattering  the  charge 
carriers  which  depresses  the  scattering  rate.  One  possible  candidate  is  a  gap  seen 
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in  the  spin  excitation  spectrum  (spin  gap)  by  a  variety  of  techinques  including:  in- 
elastic neutron  scattering  INS,  angle  resolved  photo  emission  spectroscopy  ARPES, 
and  dc  transport,  already  in  the  normal  state.  It  is  widely  believed  that  this  close 
correlation  between  the  charge  and  spin  response  in  the  cuprates  is  intimately 
related  to  superconductivity. 
4.5.4    One-Component  Analysis 

The  one-component  analysis — the  extended  Drude  model —  has  frequently 
been  used  to  analyze  the  conductivity  of  the  cuprates  [49,  64].  The  frequency 
dependent  scattering  rate  1/t{u),  at  room  temperature,  is  found  to  be  almost 
linear  over  a  very  wide  frequency  range. 

A  fundamental  feature  of  the  one-component  analysis  of  the  conductivity 
is  the  observation  that  as  the  temperature  is  lowered  below  room  temperature 
a  gap-like  depression  is  observed  in  the  frequency  dependent  scattering  rate  of 
underdoped  cuprates.  Such  a  gap  is  usually  associated  with  a  gap  which  opens 
up  in  the  underlying  bosonic  spectrum.  A  detailed  study  of  the  temperature  and 
doping  dependence  of  the  ab-plane  scattering  rate  and  c-axis  conductivity  has 
provided  valuable  insights  into  the  normal  state  phase  diagram. 

There  are  two  separate  temperature  scales  associated  with  gap  formation  in 
the  normal  state.  The  upper  onset  temperature  T*  corresponds  to  the  onset  of 
a  depression  of  the  c-axis  conductivity  and  is  also  seen  in  as  a  gap  in  the  Knight 
shift  [65]  in  NMR  measurements.  A  normal  state  gap  of  42_j^2  spmmetry  is  also 
observed  at  the  same  temperature  by  ARPES,  which  is  similar  to  that  observed 
in  the  superconducting  state  [66,  67].  This  gap  in  the  density  of  states  is  a  true 
gap.  It  also  appears  as  a  weaker  depression  of  the  low  frequency  ab-plane  scttering 
rate  extracted  from  the  extended  Drude  model.  The  lower  temperature  scale 
Ts  is  usually  a  few  degrees  above  Tc  in  underdoped  samples  and  is  marked  by  a 
reduction  in  the  NMR  relaxation  rate  l/TxT.  The  temperature      is  associated 
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with  the  observation  of  a  41  meV  neutron  resonance  in  Y-123  at  optimal  doping. 
This  represents  a  true  spin  gap  which  corresponds  to  an  ordered  spin  state.  The 
scattering  rate  shows  a  much  stronger  depression  below  Tg  than  that  observed  just 
below  T*;  charge  transport  then  becomes  coherent.  It  is  found  that  the  pseudogap 
temperatures  merge  with  Tc  in  the  optimal  to  overdoped  regions  of  the  phase 
diagram,  where  only  one  gap  structure  below  Tc  is  then  observed.  In  contrast, 
Bi-2212  appears  to  exhibit  a  pseudogap  even  in  the  overdoped  state.  The  scattering 
rate  at  high  frequency  also  becomes  superlinear.  The  linear  scattering  rate  at  high 
frequency  is  independent  of  temperature  in  the  underdoped  case,  develops  a  weak 
temperature  dependence  at  optimal  doping,  and  becomes  strongly  temperature 
dependent  in  the  overdoped  regime  where  it  becomes  superlinear.  A  review  of 
the  current  state  of  affairs  with  respect  to  the  normal  state  phase  diagram  can  be 
found  in  the  most  recent  review  by  Basov  et  al. 

Finally,  the  appearance  of  superconductivity  does  not  lead  to  any  marked 
changes  in  the  gap  feature  in  the  scattering  rate.  The  feature  increases  in  depth 
without  any  noticeable  shift  in  frequency.  This  highlights  the  marked  similarities 
between  the  normal  and  superconducting  states,  and  the  unconventional  nature  of 
the  gap  feature  seen  in  these  materials. 
4.5.5    c-Axis  Conductivity 

The  evolution  of  the  c-axis  conductivity  with  doping  has  been  investigated 
by  Uchida  et  al.  [68]  for  La-214,  and  Cooper  et  al.  [32],  Homes  et  al.  [69,  70]  and 
Bernhard  et  al.  [71]  for  Y-123.  The  c-axis  conductivity  has  a  large  number  of 
phonon  modes  which  remain  relatively  strong  up  to  very  high  doping  levels.  There 
are  a  number  of  features  worth  mentioning;  which  are  distinctly  different  from 
the  corresponding  evolution  of  the  in-plane  conductivity  with  doping:  1)  there  are 
no  features  in  the  conductivity  of  the  insulating  parent  compounds  which  can  be 
identified  with  a  charge  transfer  gap,  2)  for  La-214  doping  leads  to  a  transfer  of 
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spectral  weight  from  the  5-6  eV  range  to  lower  frequencies  which  forms  a  broad 
peak  near  2  eV,  and  a  region  of  lower  conductivity  below  0.25  eV  which  does  not 
show  a  Drude  upturn  until  it  is  overdoped  3)  the  2  eV  peak  position  does  not 
soften  with  doping  4)  the  overall  magnitude  of  the  electronic  conductivity  is  much 
smaller  than  the  corresponding  ab-plane  conductivity,  so  phonons  remain  quite 
strong  even  into  the  overdoped  regime. 

The  variation  of  the  electronic  contribution  to  the  c-axis  conductivity  with 
temperature  has  been  studied  by  the  above  authors.  It  has  been  found  that  in 
the  underdoped  cuprates  the  electronic  c-axis  conductivity  is  fairly  fiat  at  room 
temperature,  but  develops  a  gap  at  lower  frequencies  below  the  pseudogap  temper- 
atiure  T*.  The  conductivity  is  semiconducting,  and  the  depth  of  the  gap  structure 
increases  with  decreasing  temperature.  A  transition  into  the  superconducting 
state  leads  to  a  further  depression  of  the  gap.  It  is  observed  that  the  conductiv- 
ity remains  finite  down  to  the  lowest  temperatures.  This  pseudogap  structure  is 
much  harder  to  discern  in  La-214  than  in  Y-123  owing  to  the  smaller  magnitude  of 
conductivities  in  La-214. 

At  higher  doping  levels  corresponding  to  optimally  and  overdoped  regimes 
the  conductivity  become  more  metallic  in  character.  A  Drude  upturn  can  now  be 
discerned  at  low  frequency  and  the  conductivity  first  increases  with  decreasing 
temperature  above  %  at  low  frequency  and  then  decreases  in  the  superconducting 
state  as  the  condensate  forms.  This  is  similar  to  what  is  observed  for  the  ab-place 
conductivity.  A  fit  of  the  conductivity  to  an  extended  Drude  model  shows  that  the 
scattering  rate  along  the  c-axis  is  so  large  as  to  cause  the  transport  along  the  c-axis 
to  be  incoherent.  The  mean  free  paths  for  c-axis  transport  would  be  comparable  to 
the  unit  cell  dimensions  in  the  cuprates.  The  gap  feature  which  is  observed  in  the 
underdoped  case  above  Tc  now  appears  below  Tc- 
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4.6  Conclusion 

There  appear  to  be  a  number  of  features  which  are  common  to  the  carrier 
dynamics  of  the  cuprates  at  low  frequencies.  The  growth  of  spectral  weight  in 
the  mid  and  far  infrared  spectral  regions  with  doping  and  the  appearance  of 
superconductivity  are  the  most  salient  features.  The  analysis  of  optical  data  within 
the  framework  of  a  variety  of  models,  over  a  wide  doping  and  temperature  range 
have  allowed  a  fairly  precise  characterization  of  the  normal  state  phase  diagram. 

At  this  point  there  is  no  comprehensive  theory  of  the  superconducting  state. 
Two-fluid  studies  of  the  ab-plane  quasiparticle  dynamics  below  Tc  have  revealed 
a  strong  suppression  of  the  quasiparticle  scattering  rate.  This  can  be  compared 
with  the  corresponding  depression  of  the  frequency  dependent  scattering  rate 
which  begins  even  in  the  normal  state  in  underdoped  samples.  There  is  a  sizeable 
quasiparticle  spectral  weight  in  the  conductivity  at  low  temperatures  which  has 
sparked  debate  about  the  symmetry  of  the  superconducting  gap  in  these  materials. 
In  spite  of  considerable  effort  no  feature  in  the  measured  conductivity  can  be 
identified  with  a  gap  as  was  done  for  conventional  superconductors.  The  issue  of 
high-Tc  superconductivity  is  a  work  in  progress  and  infrared  spectroscopy  continues 
to  be  of  paramount  importance.  ^ 
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Figure  4-1:  Room  temperature  ab-plane  reflectivity  and  optical  conductivity  cr(a;) 
at  various  doping  levels  for  single  crystals  of  La2-xSrxCu04.  Both  scales  of  the 
reflectivity  figure  are  logarithmic  (see  ref.  [31]). 


Figure  4-2:  Eff'ective  number  of  carriers  per  unit  cell  for  single  La2_xSrxCu04  crys- 
tals at  various  doping  levels.  The  data  for  the  heavily  doped  samples  x  =  0.27  and 
X  =  0.34  are  shown  as  dot-dashed  and  dashed  curves  (see  ref.  [31]). 


CHAPTER  5 

SAMPLE  PREPARATION  AND  EXPERIMENTAL  TECHNIQUES 

5.1  Introduction 

Two  types  of  YBa2Cu306+x  high- Tc  superconducting  thin  films  were  investi- 
gated in  this  work;  both  films  are  c-axis  oriented,  with  their  c-axis  perpendicular 
to  the  substrate  surface.  The  first  film  is  grown  on  a  sapphire  substrate,  while  the 
second  is  grown  on  a  MgO  substrate.  This  chapter  deals  with  the  details  of  the 
sample  preparation  processes  and  the  characterization  of  the  films  by  dc  transport. 
This  will  be  followed  by  a  description  of  the  optical  experiments  both  with  and 
without  a  transport  current  fiowing  through  the  films. 

5.2    YBa2Cu306+x  Film  on  MgO 

The  Y-123  films  on  MgO  substrates  were  prepared  by  R.  Hadfield  at  NIST. 
The  films  were  deposited  by  pulsed  laser  deposition  from  a  stoichiometric  target  of 
Y-123  [72].  The  process  begins  with  polished  MgO  substrates  of  approximately  15 
mils  thickness,  with  length  and  width  12  mm.  Conducting  silver  paste  was  used  to 
affix  the  substrate  to  a  heater  block.  The  deposition  was  carried  out  under  vacuum 
with  the  substrate  temperature  raised  to  770  °C,  with  a  background  O2  pressure  of 
600  mTorr.  A  120  mJ  excimer  laser  operating  at  a  repetition  rate  of  10  Hz  for  54 
seconds  (two  complete  passes)  was  used  to  ablate  the  stoichiometric  target  onto  the 
heated  substrate  block.  The  chamber  was  then  flooded  with  O2  while  the  sample 
was  slowly  cooled,  to  ensure  optimal  doping.  The  resulting  film  thickness  was 
measured  to  be  500  (±20  A)  using  step  height  profilometry. 

The  chamber  was  then  vented  and  two  gold  electrodes  were  deposited  by 
dc  magnetron  sputtering  through  an  aluminum  mask.  Both  electrodes  were 
approximately  2  mm  wide  and  1200  A  thick  and  were  used  as  current  electrodes 
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(see  fig.  5-3  (b)).  The  samples  were  then  pohshed  to  remove  the  silver  paste 
from  the  substrate.  A  dummy  sample  consisting  of  bare  MgO  substrate  was 
also  produced  under  similar  conditions  in  order  to  account  for  the  effects  of  the 
substrate  in  the  infrared  spectra.  Finally,  two  silver  electrodes,  about  2000  A 
thick  and  2mm  x  2mm  were  sputtered  by  the  current  author  to  facilitate  voltage 
measurements  (see  fig.  5-3  (b)).  The  properties  of  the  samples  are  summarized  in 
table  5-1. 

The  measured  ab-plane  resistivity  of  one  of  the  YBa2Cu306+x  films  is  shown 
in  the  upper  panel  of  fig.  5-1.  The  film  exhibits  a  sharp  superconducting  transition 
near  86  K,  and  displays  a  linear  resistivity  of  the  form  p(T)  =  po  +  aT  above  100  K, 
where  a  =  1.097  //f)  cm/K  with  an  extrapolated  intercept  po  =  38.819  fiQ,  cm. 

The  advantages  of  MgO  as  a  substrate  for  the  growth  of  thin  Y-123  films 
are  twofold.  Firstly,  the  underlying  cubic  crystal  structure  of  MgO  facilitates 
epitaxial  growth  of  Y-123.  The  following  epitaxial  relationship  has  been  found, 
Y-123  [100](001)||MgO  [100](101)  [73].  Secondly,  MgO  is  reasonably  transparent 
in  the  far  infrared  up  to  330  cm~^  at  low  temperatures,  which  makes  it  ideal  for 
optical  transmittance  measurements. 


Table  5-1:  YBa2Cu306+x  thin  film  characteristics. 


Thickness 

Area 

Tc 

AT, 

Buffer 

Substrate 

(A) 

(mm^) 

(K) 

(K) 

layer 

500  A 

12x12 

86 

1.5 

none 

MgO 

600  A 

10x10 

86 

1.5 

Cerium  oxide  (200A) 

Sapphire  (AI2O3) 

5.3    YBa2Cu306+x  Films  on  Sapphire 
The  Y-123  films  on  sapphire  substrates  were  supplied  by  the  research  group 
at  Conduct  us.  The  films  were  deposited  on  the  R-plane  with  a  thin  ^  200-300 
A  buffer  layer  of  cerium  oxide  (Ce02)  by  thermal  evaporation  (TE).  The  film 
thickness  was  about  600 A  deposited  on  a  0.3  mm  thick  sapphire  substrate. 


There  are  two  main  orientations  of  sapphire  which  are  commonly  used  to  grow 
Y-123  films;  the  M-plane  (lOlO)  and  the  R-plane  (1102),  both  referenced  with 
respect  to  hexagonal  axes.  The  hexagonal  (trigonal)  crystal  structure  of  sapphire 
makes  it  anisotropic.  The  optic  axis  lies  in  the  M-plane,  but  at  an  angle  of  57.5667° 
to  the  R-plane  [74].  This  makes  R-cut  sapphire  inconvenient  for  the  growth  of 
Y-123  films  in  cases  where  the  dielectric  anisotropy  is  undesirable.  The  following 
epitaxial  relationships  exist  for  Ce02  buffered  Y-123  films  on  sapphire  substrates: 
CeO2[110](001)  II  Al2O3[1011](1102)  and  Y-123[100](001)  ||  CeO2[110](001) 
[75].  The  main  advantages  of  sapphire  substrates  are  their  hardness,  resistance  to 
corrosive  chemicals  and  high  thermal  conductivity  comparable  to  that  of  copper  at 
cryogenic  temperatures.  They  are  also  very  transparent  at  low  temperature  in  the 
far  infrared  region  below  350  cm~^. 

The  resistivity  of  the  film  as  a  function  of  temperature  was  measured  using 
a  standard  4-probe  ac  technique.  The  lower  panel  of  fig.  5-1  shows  that  the 
superconducting  transition  is  fairly  sharp  with  Tc  =  86  K.  Above  100  K  the 
resistivity  follows  a  liner  law  p{T)  =  po  +  aT  similar  to  that  found  for  the  films  on 
MgO  substrates  but  with  a  higher  normal  state  resistivity.  Table  5-1  summarizes 
the  characteristics  of  the  film. 

5.4  Experiment 

5.4.1    Cryogenic  Apparatus 

In  order  to  perform  infirared  spectroscopic  measurements  at  cryogenic  temper- 
atures, it  is  necessary  to  cool  and  regulate  the  temperature  of  the  sample  under 
investigation.  The  samples  are  usually  glued  to  the  back  surface  of  an  oxygen- 
free  copper  plate  with  a  tapered  hole  (see  fig.  5-2).  Oxygen-firee  copper  has  a 
very  low  brittle  point  and  high  thermal  conductivity  which  makes  it  ideal  for  low 
temperature  applications. 


There  are  four  main  parts  of  the  coohng  system:  a  flow  cryostat,  hehum 
dewar,  transfer  hne  and  temperature  controller.  All  parts  work  in  tandem  to 
maintain  a  stable  sample  temperature.  The  sample  plates  are  mounted  onto  a 
screw-on  attachment  which  screws  into  the  cryostat  cold  finger.  A  small  silicon 
diode  thermometer  is  inserted  into  a  small  hole  drilled  into  the  sample  holder  which 
measures  the  temperature.  To  reduce  the  effects  of  thermal  radiation  heating,  a 
metal  shield  is  screwed  onto  the  cold  stage  of  the  cryostat.  The  whole  arrangement 
is  inserted  into  a  shroud  which  is  part  of  the  reflectivity  or  transmittance  stage. 
The  shroud  is  sealed  with  infrared  transparent  windows  which  allow  light  to  reach 
the  sample.  The  shroud  is  then  evacuated  to  a  pressure  10~^  Torr  or  less.  Helium 
gas  is  then  transferred  from  the  hehum  dewar  to  the  cryostat  through  the  transfer 
line.  This  gas  flow  cools  the  cold  finger  and  hence  the  sample  by  conduction.  The 
temperature  is  maintained  at  a  set  value  by  a  temperature  controller  which  reads 
the  temperature  of  the  diode  thermometer  and  supplies  the  requisite  power  to  a 
heater  on  the  cold  finger.  Depending  on  the  desired  rate  of  coohng  the  helium 
flow  rate  may  be  increased  or  decreased  by  a  flow  meter  which  regulates  the  flow 
through  the  cryostat.  The  model  temperature  controller  used  in  this  work,  the 
Lakeshore  model  330  auto-tuning  temperature  controller,  automatically  adjusts  the 
control  parameters  at  each  temperature  and  gives  fairly  stable  temperatures  within 
0.01  K.  The  diode  thermometer  when  weU  cahbrated  has  an  accuracy  of  ±  1  K. 
5.4.2    Reflectivity  Measurements 

The  reflectivity  measurements  on  the  Y-123  samples  were  performed  on  the 
Bruker  113-V  FTIR  spectrometer  (see  fig.  3-2  p.  35).  The  measiu-ements  were  con- 
fined to  the  far  infrared  range  35-350  cm~\  and  to  temperatures  between  300  and 
40  K.  The  sample  was  mounted  on  a  holder  similar  to  that  shown  schematically 
in  fig  5-2(b).  An  aluminum  mirror,  consisting  of  2000  A  of  aluminum  deposited 
on  a  microscope  slide  served  as  a  reference.  The  reference  was  affixed  to  a  plate 
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of  the  same  dimensions  as  the  sample  holder.  Both  sample  and  reference  holders 
were  mounted  on  the  cryostat  in  a  way  which  allowed  them  to  be  interchanged 
in  the  infrared  beam  (see  fig.  5-2  (a)).  An  arc  lamp  (Hg-arc)  was  used  as  a  light 
source  and  a  helium  cooled  bolometer  as  a  detector.  The  absolute  reflectivity  is 
found  by  dividing  the  measured  spectrum  of  the  sample  in  the  beam  by  that  of  the 
reference  mirror  in  the  beam.  Although  the  absolute  reflectivity  of  aluminum  is 
not  100  %  it  is  greater  than  99  %  over  a  very  wide  frequency  range  extending  into 
the  near  infrared.  The  measured  reflectivity  is  corrected  by  multiplying  it  by  the 
absolute  reflectivity  of  the  reference.  The  reflectivity  of  the  bare  MgO  substrate 
was  also  measured  to  account  for  substrate  effects  in  the  data  for  the  film  on  MgO. 
Similarly,  one  of  the  films  on  sapphire  was  etched  with  nitric  acid,  and  the  reflec- 
tivity of  the  Ce02/sapphire  combination  was  measured.  Chapter  6  deals  with  the 
difficulties  encountered  in  trying  to  analyze  both  the  film/Ce02/sapphire  data  and 
the  Ce02/sapphire  data  for  anisotropic  R-plane  sapphire. 
5.4.3    Transmittance  Measurements 

The  transmittance  of  both  films,  and  their  substrates  were  also  measured  at 
the  same  temperatures  and  over  the  same  frequency  range  as  the  refiectivity.  The 
transmittance  is  found  in  a  similar  manner  to  the  refiectivity.  Fig.  5-3  (a)  shows 
the  schematic  setup  for  transmittance.  The  sample,  and  empty  aperture  which 
now  serves  as  a  reference,  are  mounted  at  right  angles  on  the  tip  of  the  rotating 
cryostat.  The  absolute  transmittance  is  determined  by  dividing  the  measured 
single  beam  spectrum  of  the  sample  by  that  of  the  empty  reference  hole.  The 
transmittance,  unUke  the  reflectivity  is  measured  at  near  normal  incidence.  The 
angle  of  incidence  for  the  refiectivity  measurements  described  in  the  previous 
section  is  about  11°. 

To  investigate  the  effects  of  transport  currents  on  the  superconducting  state 
electrodynamics,  surface  electrodes  for  both  current  and  voltage  were  sputtered 
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onto  both  films.  The  process  used  for  the  deposition,  and  the  arrangement  of 
electrodes  for  the  film  on  MgO  was  previously  described,  and  is  shown  in  fig.  5-3 
(b).  The  film  on  sapphire  was  also  patterned  using  a  small  amount  of  dilute  nitric 
acid  and  a  paper  mask  cut  into  the  form  of  a  half  cross.  The  geometry  of  the 
patterned  film  is  shown  in  fig.  5-3  (b).  The  two  outer  silver  electrodes  used  to 
carry  current  were  2  mm  wide,  5  mm  long  and  2000  A  thick.  Each  of  the  voltage 
electrodes  is  2mm  long,  1mm  wide  and  2000  A  thick.  The  center  to  center  distance 
between  the  voltage  electrodes  is  about  3  mm.  The  center  to  center  distance 
between  the  voltage  electrodes  for  the  film  on  MgO  is  about  4  mm.  Two  current 
sources  were  used  to  supply  current,  together  they  covered  a  range  from  1  nA-1  A. 
The  Kiethly  model  220  programmable  current  source  covers  the  lower  range  1  nA- 
100  mA  and  the  model  228  current- voltage  source  covers  the  range  4  mA-1  A. 
Voltage  measurements  were  made  with  a  Kiethly  model  195  digital  multimeter. 
Contact  was  made  by  soldering  gold  wire  to  the  electrodes  using  indium  solder  and 
a  low  temperature  soldering  iron.  Indium  is  known  to  produce  contacts  of  much 
lower  resistance  than  silver  epoxy. 

The  electrooptic  transmittance  was  performed  in  two  separate  regimes;  near  Tc 
at  75  K  and  at  low  temperature — 40  K.  In  order  to  ensure  that  the  samples  were 
in  the  superconducting  state  near  Tc  the  samples  were  cooled  to  80  K  and  a  small 
current  of  about  1  /xA  was  passed  through  the  films  and  the  voltage  reading  on 
the  multimeter  was  observed.  It  was  found  that  in  both  cases  there  was  a  voltage 
reading,  which  may  be  related  to  the  possibility  that  the  film  may  be  at  a  slightly 
higher  temperature  than  the  sample  holder.  The  temperature  was  further  reduced 
to  75  K  and  the  procedure  was  repeated.  The  voltage  reading  on  the  multimeter 
was  now  zero;  implying  that  the  films  were  in  the  superconducting  state. 

The  transmittance  of  the  film  on  MgO  was  measured  as  a  function  of  current 
at  75  K  at  25,  75  and  200  mA.  The  voltage  was  observed  each  time.  The  current 
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was  also  increased  to  300  mA  at  which  point  the  substrate  cracked.  This  may  be 
as  a  result  of  ohmic  dissipation  at  the  contacts  which  produced  a  temperature 
gradient  between  the  film  side  and  backside  of  the  substrate  which  caused  the 
substrate  to  crack.  The  same  procedure  was  carried  out  for  the  film  on  sapphire  at 
75  K  except  that  the  currents  passed  were  1,  5,  10,  15  and  20  mA,  since  the  critical 
current  was  between  10  and  15  mA.  The  results  of  both  these  measurements  are 
presented  in  chapter  7. 

The  transmittance  of  the  film  on  sapphire  was  also  measured  as  a  function  of 
current  at  40  K.  Four  currents  of  50,  80,  100  and  150  mA  were  passed.  At  150  mA 
the  voltage  began  to  increase  rapidly,  far  beyond  that  expected  for  the  resistive 
state  at  the  same  temperature — which  appeared  to  be  a  sign  of  contact  heating. 
This  will  also  be  discussed  in  greater  detail  in  chapter  7. 
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Figure  5-1:  The  temperature  dependent  resistivity  of  two  YBa2Cu306+x  films. 
Upper  panel:  the  ab-plane  resistivity  as  a  function  of  temperature  for  a  500  A 
YBa2Cu306+x  film  on  a  MgO  substrate.  Lower  panel:  the  ab-plane  resistivity  of 
second  YBa2Cu306+x  film,  600  A  thick,  on  a  sapphire  substrate. 
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two  sample  holders  with 
sample  and  reference  minor 


(a)  light  path  diagram 


fitont  cross  section 

(b)  sample  holder 

Figure  5-2:  Reflectivity  stage  and  sample  holders  used  in  this  experiment,  (a)  A 
schematic  diagram  of  the  reflectivity  stage  and  light  path  showing  light  incident  on 
rotatable  sample  holders,  (b)  A  schematic  diagram  of  the  oxygen  free  copper  holder 
used  to  mount  samples. 
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Two  plates  with  sample, 
and  empty  aperture  as  ref- 
erence on  rotating  cryostat 


(a) 


(b)  (c) 

Figure  5-3:  A  schematic  of  the  transmittance  holder  and  electrode  configuration 
for  transmittance  measurements,  (a)  Schematic  diagram  of  light  path  for  transmit- 
tance measurement,  (b)  Surface  electrode  configuration  used  for  infrared  electroop- 
tic  measurements  for  Y-123  films  on  MgO.  (c)  Patterned  half-cross  and  surface 
electrode  arrangement  used  for  infrared  electrooptic  measurements  of  Y-123  on 
sapphire. 


CHAPTER  6 

FAR-INFRARED  STUDIES  OF  YBasCusOe+x  FILMS 
6.1  Introduction 

This  chapter  deals  with  a  detailed  characterization  of  the  far-infrared  prop- 
erties of  the  two  Y-123  thin  films  which  were  described  in  the  previous  chapter. 
Both  the  reflectivity  Ti{uj)  and  transmittance  T{u]  were  measured  over  a  frequency 
range  of  40-350  cm~^  and  over  a  temperature  range  from  40-300  K.  The  charac- 
teristics of  both  films  have  been  listed  in  table  5-1  on  p.  57;  while  the  dc  transport 
properties  are  shown  in  fig.  5-1  on  p.  63.  The  films  are  highly  ab-plane  oriented, 
and  the  E]-field  of  the  incident  light  lies  within  the  plane.  All  the  optical  properties 
which  are  reported  here  are  therefore  ab-plane  averages. 

The  dielectric  properties  of  the  film  on  magnesium  oxide  in  the  normal  state 
above  Tg  are  analyzed  within  the  framework  of  a  simple  two-component  model.  We 
have  been  able  to  obtain  fairly  good  fits  to  the  transmittance  by  assuming  that 
the  low  frequency  response  is  dominated  by  a  free-carrier  Drude  band,  while  the 
response  at  higher  frequency  is  influenced  primarily  by  a  very  broad  midinfrared 
band  which  is  fitted  with  a  series  of  Lorentzians.  The  results  obtained  from  this 
two-component  approach  agree  very  well  with  the  results  obtained  from  the 
Kramers-Kronig  analysis  of  high  quahty  single  crystals.  Below  Tc  the  dielectric 
features  of  the  superconducting  state  are  studied  with  a  two-fluid  model.  In  this 
model,  part  of  the  spectral  weight  associated  with  the  Drude  band  is  transferred  to 
the  condensate,  which  appears  as  a  pole  in  the  conductivity  at  zero  frequency.  The 
midinfrared  band  is  found  to  be  fairly  temperature  independent. 

The  film  on  sapphire  will  also  be  discussed.  The  chief  drawback  of  the  analysis 
is  the  difficulty  of  extracting  the  optical  constants  of  the  film  independently  of  the 
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substrate.  Due  to  the  anisotropy  of  sapphire,  a  uniaxial  transmittance  model  simi- 
lar to  that  described  in  section  2.4  is  used  to  fit  the  substrate  transmittance  data. 
The  fits  to  this  model  are  in  reasonable  agreement  with  the  measured  data,  over 
a  wide  frequency  range.  The  generic  features  of  both  the  normal  and  supercon- 
ducting states  of  this  film  on  sapphire  are  shown  to  be  quite  similar  to  those  of  the 
film  on  MgO.  The  normal  state  data  of  this  film  was  fitted  using  a  two-component 
model,  and  a  two-fluid  model  was  used  to  analyze  the  superconducting  state  data. 
The  results  of  this  analysis  will  be  discussed  in  later  sections. 

6.2    Results  and  Discussion 

6.2.1    MgO  Substrate 

Magnesium  oxide  possesses  complex  temperature-dependent  lattice  vibration 
spectra  in  the  far-infrared  region.  We  concentrated  on  the  low  frequency  region 
below  350  cm~\  as  MgO  is  opaque  at  all  temperatures  above  this  frequency  due  to 
a  strong  infrared  active  transverse  optic  (TO)  phonon  mode  near  400  cm~^  The 
measurements  were  made  at  2  cm~^  resolution,  and  the  fringes  due  to  the  multiple 
internal  reflections  in  the  substrate  were  carefully  smoothed. 

Fig.  6-1  shows  the  measured  reflectivity  and  transmittance  of  a  0.3  mm  thick 
MgO  substrate  at  several  temperatures.  Below  100  cm~^  the  absorption  is  small 
at  all  temperatures,  as  the  transmittance  approaches  60%  while  the  corresponding 
reflectivity  approaches  40%.  The  transmittance  in  the  region  extending  from  100 
cm"'  to  about  300  cm"'  shows  a  very  marked  temperature  dependence.  The 
spectra  are  characterized  by  structure  near  100,  150  and  290  cm"'  which  have  been 
assigned  to  two  phonon  difference  band  processes  which  lead  to  strong  absorption 
at  room  temperature,  and  which  are  frozen  out  at  lower  temperatures.  In  these 
processes  a  photon  is  absorbed,  and  an  acoustic  phonon  is  annihilated  with  the 
creation  of  an  optic  phonon.  The  possible  assignments  which  correspond  to  regions 
where  the  phonon  densities  of  states  are  large  are  the  W-point  TA— TO  transition 


68 

at  90  cm~\  the  X'-point  LA — LO  transition  at  95  cm"^  and  the  X-point  TA — TO 
transitions  at  150  and  255  cm~^  As  the  temperature  is  lowered  there  is  an  overall 
increase  in  both  the  the  transmittance  and  the  reflectivity  as  the  multiphonon 
processes  are  suppressed.  The  reflectivity  near  350  cm~^  is  fairly  insensitive  to 
temperature  as  is  expected  for  a  strong  TO  mode  which  dominates  the  reflectivity 
above  this  frequency. 

In  order  to  account  for  the  complex  MgO  absorption  processes  which  appear 
in  the  measured  spectra  for  the  film  on  MgO,  we  fitted  the  substrate  transmittance 
with  a  large  number  of  Lorentzians.  The  frequency  dependent  dielectric  function 
e{uj)  was  assumed  to  have  the  form 

S  ijJ^ 

e(a;)  =  6oo  +  V^  y^-  ,  (6.1) 

^  Uj  —u''  —  iwyj 

where      is  the  high  frequency  dielectric  constant,  Sj  is  the  mode  strength,  uj  is 
the  mode  frequency  and  jj  is  the  damping  constant.  The  high  frequency  dielectric 
constant  e^o  was  fixed  at  3.01  at  all  temperatures  as  was  done  in  ref.  [76].  We 
initially  started  with  the  parameters  which  were  suggested  by  ref.  [76]  but  we 
found  that  it  was  impossible  to  obtain  good  fits  in  the  multiphonon  spectral  region 
without  introducing  additional  osciUators.  This  prescription  of  simulating  the 
multiphonon  difference  band  processes  by  a  large  number  of  Lorentzians  allowed  us 
to  account  for  the  substrate  effects  in  the  film  data.  Table  6-1  shows  the  oscillator 
parameters  which  were  used  in  the  fits  to  the  substrate  transmittance  data  at  all 
temperatures.  The  mode  strengths  in  eqn.  (6.1)  are  written  in  terms  of  a  plasma 
frequency  Upj  =  uj        in  the  table.  The  lower  panel  of  fig.  6-1  shows  the  fits  to 
the  substrate  transmittance  data  using  this  procedure. 
6.2.2    Sapphire  Substrate 

In  contrast  to  MgO,  sapphire  belongs  to  a  trigonal  (hexagonal)  crystal  system. 
This  imphes  that  any  macroscopic  quantity  such  as  the  dielectric  constant  is 


described  by  a  2nd-rank-tensor  which  has  two  independent  parts.  A  complete 
characterization  of  the  dielectric  response  requires  two  independent  measurements 
of  e{u;);  one  with  the  electric  field  parallel  to  the  crystalline  c-axis  and  the  other 
with  the  electric  vector  anywhere  in  the  plane  perpendicular  to  the  c-axis,  since 
all  such  directions  are  optically  equivalent.  The  particular  cut  used  here — R-cut 
sapphire —  makes  it  impossible  to  perform  such  a  characterization  by  polarized 
optical  measurements.  The  c-axis  (optic  axis)  makes  an  angle  of  57.5667°  with  the 
R-plane.  We  performed  unpolarized  reflectance  and  transmittance  measurements 
below  350  cm~\  since  sapphire  has  a  large  number  of  infrared  active  phonon  modes 
above  380  cm~^  and  is  opaque  above  this  frequency.  The  basic  theory  of  light 
propagation  through  an  anisotropic  crystal,  based  on  a  uniaxial  transmittance 
model  at  normal  incidence  was  presented  in  chapter  2.  Here,  we  will  concentrate 
mainly  on  the  transmittance,  as  these  measurements  are  done  at  more  near  normal 
incidence,  which  was  assumed  in  the  theory. 

The  measured  reflectivity  and  transmittance  of  a  thin  Ce02  buffer  layer 
(300  A)  on  a  0.5  mm  thick  sapphire  substrate  are  shown  in  fig.  6-2.  The  fringes, 
which  have  not  been  smoothed,  are  due  to  multiple  internal  reflections  which 
occur  within  the  sapphire  substrate.  The  figure  shows  that  below  100  cm"^  the 
average  reflectivity  approaches  40%  while  the  average  transmittance  approaches 
60%.  This  implies  that  there  is  little  absorption  at  very  low  frequency  at  all 
temperatures.  A  salient  feature  of  all  the  spectra  is  the  appearance  of  beats  which 
increase  in  contrast  and  extend  to  higher  frequencies  as  the  temperature  is  lowered. 
The  beats  result  from  a  superposition  of  two  oscillatory  terms  in  the  reflectivity 
or  transmittance  associated  with  different  components  of  the  dielectric  tensor. 
The  peak  in  the  reflectivity  and  the  dip  in  the  transmittance  at  about  272  cm"^ 
corresponds  to  a  single  infrared  active  TO  phonon  mode  in  the  Ce02  buffer  layer 
[77-81]. 
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As  was  observed  in  the  case  of  MgO  the  transmittance  above  100  cm~-^ 
increases  drastically  as  the  temperature  is  lowered  from  room  temperatm-e  to  40 
K;  while  the  increase  in  the  reflectivity  is  less  drastic  over  the  same  temperature 
range.  This  drastic  suppression  of  the  absorption  above  100  cm~^  with  decreasing 
temperature,  may  be  related  to  multiphonon  difference-band  processes  in  the 
sapphire  substrate  similar  to  those  observed  in  MgO.  The  main  obstacle  to  any 
analysis  of  the  optical  data  in  this  frequency  range  is  the  lack  detailed  studies  of 
sapphire  over  this  range  [82-85].  The  majority  of  the  previous  infrared  studies  were 
focussed  on  the  phonon  modes  in  the  frequency  range  above  380  cm~^  [86-88]. 

Sapphire  belongs  to  the  D^d  space  group.  Group  theory  predicts  4  infrared 
active  phonon  modes  of  Eu  symmetry  for  the  electric  field  perpendicular  to  the 
c-axis  (o-ray),  and  2  infrared  active  modes  of       symmetry  for  the  E-field  parallel 
to  the  c-axis  (e-ray).  We  represented  the  dielectric  functions  of  the  o  and  e  waves 
by  Lorentz  models  of  the  form 


where  eo,e  is  the  dielectric  function  for  the  o  or  e  wave,  Coo  is  the  high  frequency 


is  the  damping  constant.  We  tried  to  fit  the  transmittance  of  Ce02  on  sapphire 
using  the  parameters  suggested  in  ref.  [79]  for  Ce02  and  the  parameters  in  ref. 
[88]  for  sapphire.  All  four  infrared  active  phonon  modes  as  well  as  two  electronic 
bands  at  higher  frequency  were  included  in  eo]  while  the  two  predicted  modes, 
and  two  electronic  bands  were  included  in  Cg.  In  addition  to  these,  we  added  a 
broad  oscillator  to  to  at  low  frequency  and  centered  at  209  cm"^  to  account  for  the 
strong  absorption,  particularly  at  higher  temperatures.  The  choice  of  the  center 
frequency  of  this  extra  mode  was  motivated  by  the  observation  of  a  weak  dip  in  the 
transmittance  at  approximately  the  same  frequency.  This  can  be  seen  most  easily 


(6.2) 


dielectric  constant,  Sj  is  the  oscillator  strength,  u>j  is  the  mode  frequency,  and  7j 
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in  the  300  and  200  K  spectra  in  the  lower  panel  of  fig.  6-2.  We  found  that  it  was 
not  possible  to  fit  the  substrate  transmittance  without  including  this  broad  mode 
at  all  temperatures.  . 

The  fits  to  the  transmittance  are  shown  in  fig.  6-3  along  with  the  measured 
data,  which  is  repeated  here  for  the  sake  of  comparison.  Tables  6-2-6-4  show  the 
parameters  which  were  used  in  the  fits  to  the  sapphire  substrate  and  Ce02  buffer 
layer  at  all  temperatures.  The  fits  to  the  data  replicate  the  beat  pattern  quite 
accurately,  although  the  fringe  contrast  of  the  fits  is  greater  than  the  actual  data  at 
all  temperatures.  The  difference  is  larger  at  lower  temperature  where  the  sapphire 
absorption  is  lower  and  the  dispersive  nature  of  the  substrate  dominates.  The  main 
reason  for  this  is  that  the  fits  assume  that  the  measurement  was  done  at  infinite 
resolution.  The  actual  data  was  measured  at  1  cm~^  resolution  and  so  the  fringe 
contrast  is  lower.  In  the  frequency  range  between  250  and  270  cm~\  just  below 
the  Ce02  mode  the  fits  differ  by  1-2%.  This  may  be  obscured  by  the  fringes  in  fig. 
6-3.  The  fits  fall  slightly  below  the  data  in  this  frequency  range.  The  difference  is 
less  severe  at  lower  temperatures.  The  fitted  transmittance  in  this  region  is  fairly 
sensitive  to  the  parameters  of  the  broad  mode  which  was  introduced  at  209  cm~^ 
to  account  for  the  extra  absorption  at  higher  temperatures.  This  difference  between 
the  fits  and  the  actual  data  of  the  sapphire  substrate  will  lead  to  errors  in  the 
extracted  optical  constants  of  the  film  data. 
6.2.3    YBaaCuaOe+x  on  MgO 

The  measured  reflectance  and  transmittance  of  the  500  A  film  on  MgO  are 
shown  in  the  upper  and  lower  panels  of  fig.  6-4.  The  reflectivity  falls  well  below 
unity  at  all  frequencies  in  the  normal  state.  As  the  temperature  is  lowered  the 
reflectivity  below  150  cm~^  increases  drastically,  and  approaches  unity  at  the 
two  lowest  temperatures  in  the  superconducting  state.  The  reflectivity  at  all 
temperatures  is  fairly  flat  and  featureless,  which  can  be  attributed  to  the  strong 


72 

free  carrier  contribution  at  all  temperatures.  The  free  carriers  screen  the  phonons 
which  are  normally  observed  in  the  substrate  and  in  the  film.  The  transmittance 
is,  however,  much  lower  and  never  exceeds  10%  at  any  temperature.  The  spectra, 
unlike  the  reflectivity  are  characterized  by  significant  structure,  which  is  due  to  the 
phonons  in  the  MgO  substrate.  As  the  temperature  is  lowered,  the  transmittance 
below  100  cm~^  is  reduced  and  approaches  zero  in  the  superconducting  state  for 
T<Tc. 

6.2.4    YBa2Cu306+x  on  Sapphire 

Fig.  6-5  shows  the  raw  reflectivity  and  transmittance  of  the  600  A  film  on 
sapphire  at  several  temperatures.  The  spectra  have  not  been  smoothed  and  are 
characterized  by  beats  below  270  cm~^  The  beats  are  due  to  multiple  internal 
reflections  in  the  sapphire  substrate  as  was  described  earlier.  The  peak  in  the 
reflectivity — and  dip  in  the  transmittance — at  about  272  cm~^  corresponds  to 
a  phonon  mode  in  the  underlying  Ce02  buffer  layer.  At  room  temperature  the 
fringe- averaged  low  frequency  reflectivity  is  about  55%.  The  overall  reflectivity 
increases  at  lower  temperatures,  with  the  low  frequency  reflectance  showing  an 
upturn  which  approaches  100%  at  40  K  in  the  superconducting  state.  The  onset 
of  superconductivity  below  Tc  may  also  be  inferred  from  the  changes  in  the  low 
frequency  characteristics  of  the  fringe  pattern  in  the  reflectivity.  The  fringe  contrast 
below  100  cm  ~^  is  strongly  depressed  at  the  two  lower  temperatures,  80  and  40 
K  in  the  superconducting  state  as  compared  to  the  fringe  contrast  in  the  normal 
state.  The  transmittance  at  room  temperature  increases  steadily  below  the  infrared 
edge  of  transparency  of  the  substrate,  up  to  about  17%  at  the  lowest  frequency. 
There  is  a  reduction  of  the  transmittance  below  100  cm~^  at  lower  temperatures 
with  the  transmittance  approaching  zero  below  Tc-  At  higher  frequency  the 
transmittance  is  enhanced  at  low  temperature  as  the  substrate  becomes  more 
transparent. 
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6.2.5    Two-Fluid  Model  .  • 

In  order  to  account  for  the  frequency  and  temperature  dependence  of  the 
observed  transmittance  of  both  films,  we  have  employed  a  two-component  model  in 
the  normal  state  and  a  two-fluid  model  in  the  superconducting  state.  The  models 
have  been  outlined  in  chapter  4. 

The  two-component  model  assumes  that  the  complex  dielectric  function  of  the 
film  is  composed  of  three  terms:  a  narrow  Drude  band  centered  at  zero  frequency,  a 
number  of  broad  Lorentzians  in  the  midinfrared,  and  a  high  frequency  contribution. 
The  dielectric  function  may  then  be  written  as: 

e{u)  =  €d{uj)  +  eiiuj)  +  {T  >  Tc),  (6.3) 

where  ei:i{u)  is  the  Drude  term,  eL{u)  is  a  sum  of  Lorentzians  which  are  centered  in 
the  midinfrared  and  which  account  for  the  broad  midinfrared  absorption,  and 
is  the  high  frequency  term  which  accounts  for  excitations  at  frequencies  exceeding 
10,000  cm-\ 

In  the  superconducting  state  it  is  assumed  that  a  certain  fraction  of  the  Drude 
carriers  condense  to  form  the  superconducting  condensate,  while  the  effect  on  the 
midinfrared  band  is  negligible.  We  have  used  a  two-fluid  model  dielectric  function 
of  the  form: 

e{u)  =  (1  -  fs)eD{uj)  +  fstsiuj)  +  e^Ccj)  +        (T  <  T^),  (6.4) 

where  eoi^j),  ei{u))  and  e^o  are  the  same  as  the  normal  state  values  in  eqn.  6.3, 
while  fs  is  the  fraction  of  the  free  carrier  density  which  condenses  and 


e,H  =  -^  +  ,I^iH.  (6.5) 


The  plasma  frequency  cjpo  here  is  the  same  as  the  the  Drude  plasma  frequency  in 
the  normal  state. 


The  fits  to  the  transmittance  of  the  film  on  MgO,  and  the  film  on  sapphire 
are  shown  in  the  lower  panel  of  fig.  6-4  and  in  fig.  6-6  respectively.  In  each  case 
a  single  Drude  oscillator  with  plasma  frequency  Upjj  and  damping  1/r  and  three 
Lorentzians  with  parameters  {ujpj,Uj,^j)  where  Upj,  Uj  and  7j  are  the  plasma 
frequency,  center  frequency  and  damping  of  the  jth  Lorentzian  were  used  to  fit 
the  transmittance.  We  used  the  parameters  suggested  by  reference  [3]  for  the 
Lorentzian  oscillators  and  also  fixed  eoo  =  25  as  was  done  in  ref  [9].  We  obtained 
fairly  good  fits  at  all  temperatures  for  the  film  on  MgO  by  varying  only  uipu, 
1/r  and  fg  as  well  as  the  damping  71  of  the  lowest  midinfrared  oscillator  while 
fixing  the  remaining  midinfrared  parameters.  It  was  also  found  that  Upo  and  71 
were  fairly  independent  of  temperature,  within  experimental  uncertainty.  The 
parameters  which  were  used  to  fit  the  film  data  in  the  normal  and  superconducting 
states  are  shown  in  tables  6-5  and  6-6. 

A  similar  procedure  was  applied  to  the  film  on  sapphire,  and  as  can  be  seen 
in  fig.  6-6  the  fits  are  reasonably  good  at  all  temperatures.  The  fringe  pattern  of 
the  raw  data  is  nearly  centered  on  the  fringe  pattern  of  the  fits.  In  these  fits  all 
the  midinfrared  parameters  were  fixed,  while  the  Drude  plasma  frequency  Upo, 
scattering  rate  1/r,  and  superfiuid  fraction  fg  in  the  superconducting  state  were 
varied  to  achieve  the  best  least  squares  fit  to  the  data.  It  is  nearly  impossible  to 
fit  the  fringe  pattern  very  accurately.  This  can  lead  to  substantial  errors  in  some 
of  the  oscillator  parameters  used  to  fit  the  film  data.  In  particular,  the  damping 
constant  71  of  the  lower  midinfrared  oscillator  at  center  frequency  290  cm~^  will 
accumulate  substantial  errors  if  it  is  varied  to  achieve  better  fits  in  this  region. 

Tables  6-7  and  6-8  shows  some  of  the  oscillator  parameters  which  were  used 
to  fit  the  transmittance  of  the  film  on  sapphire.  The  Drude  plasma  frequency  of 
the  film  on  sapphire  is  almost  40%  lower  than  the  plasma  frequency  of  the  film 
on  MgO.  There  are  two  main  reasons  for  this.  One  reason  for  the  difference  is 
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associated  with  the  uncertainties  in  the  film  thicknesses  which  are  used  in  the  fits. 
Such  uncertainties  produce  errors  in  the  the  plasma  frequencies  of  the  oscillators 
used  in  the  fits.  Another  possible  cause  of  the  discrepancy  may  simply  be  the 
higher  quality  of  the  film  grown  on  MgO. 
6.2.6    Free  Carrier  Relaxation  and  Condensate 

The  observed  changes  in  the  reflectivity  and/or  transmittance  of  the  films 
with  temperatiure  can  be  easily  understood  in  terms  of  the  strong  temperature 
dependence  of  the  optical  conductivity.  Fig.  6-7  shows  the  real  and  imaginary 
parts  of  the  optical  conductivity  of  the  film  on  MgO  extracted  from  the  two- 
component  fits  to  the  transmittance.  The  real  part  of  the  optical  conductivity 
ai{uj)  at  frequencies  below  200  cm~^  increases  as  the  temperature  is  lowered  in 
the  normal  state.  This  transfer  of  spectral  weight  to  lower  frequencies  is  due  to 
a  linear  decrease  in  the  free  carrier  relaxation  rate  1/r  with  temperature.  The 
transmittance  in  the  normal  state  is  strongly  infiuenced  by  ai  (tu)  which  dominates 
the  inductive  response  0-2(0;).  The  enhanced  low  frequency  conductivity  at  lower 
temperatures  leads  to  a  corresponding  decrease  in  the  low  firequency  transmittance. 
At  temperatures  below  Tc  the  onset  of  superconductivity  is  marked  by  a  strong 
enhancement  of  a2{u))  due  to  the  formation  of  the  condensate.  There  is  a  transfer 
of  spectral  weight  to  the  origin  which  forms  a  delta  function  in  ai{uj)  which 
represents  the  response  of  the  condensate.  This  redistribution  of  spectral  weight 
leads  to  a  suppression  of  ai  (w)  at  low  frequencies  as  predicted  by  the  sum  rule. 
Since  the  real  and  imaginary  parts  of  the  conductivity  are  related  by  causality  a 
delta  function  in  the  real  part  of  the  conductivity  produces  a  l/u)  dependence  in 
the  imaginary  part.  As  a;  — >  0,  a2{oj)  approaches  zero  in  the  normal  state  while  it 
diverges  as  A/u  in  the  superconducting  state  where  it  now  dominates  ai{u).  The 
quantity  A  is  related  to  the  strength  of  the  condensate  and  the  London  penetration 
depth.  The  observed  downturn  of  the  transmittance  in  the  superconducting  state, 
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shown  in  fig.  6-4  is  due  to  the  l/u  dependence  of  0-2(0;),  which  represents  the 
screening  efl'ect  of  the  condensate.  Fig.  6-8  shows  the  corresponding  conductivities 
of  the  film  on  sapphire.  The  only  difference  between  this  film  and  the  film  on  MgO 
are  the  lower  levels  of  the  conductivities  of  this  film.  This  is  associated  with  the 
lower  free  carrier  densitiy  of  this  film  which  is  reflected  in  the  lower  Drude  plasma 
frequency  cOpo- 

We  have  extracted  the  temperature  dependent  scattering  rate  of  both  films. 
The  scattering  rates  were  calculated  from  the  measured  four-probe  resistivities 
using  the  equation: 

=  ^Pdc  (6.6) 

The  results  are  shown  in  fig.  6-9  where  we  have  also  shown  the  relaxation  rates 
obtained  directly  from  the  two-component  fits.  The  figure  shows  that  the  free 
carrier  relaxation  rate  varies  linearly  with  temperature  in  the  normal  state.  This 
is  a  unique  property  of  the  carrier  dynamics  of  the  cuprates.  The  infrared  also 
probes  the  quasiparticle  relaxation  rate  below  Tc  which  is  invisible  to  dc  transport. 
The  symbols  in  the  figure  show  that  the  carrier  relaxation  rate  extracted  from  a 
two-fluid  fit  falls  faster  than  linearly  below  Tc.  Such  behavior  has  previously  been 
observed  for  YBa2Cu306+x  [61]  and  Bi2Sr2CaCu208  [63].  It  is  widely  believed 
that  the  suppression  of  the  scattering  rate  below  Tc  is  related  to  a  gap  which 
opens  up  in  the  excitation  spectrum  which  is  responsible  for  scattering  the  free 
carriers.  The  close  correspondence  in  the  energy  scales  of  characteristic  features 
which  appear  in  the  charge  and  spin  response  of  the  cuprates,  suggests  that  the 
spin  excitation  spectrum  may  be  responsible  for  this  unusual  behavior.  The  spin 
excitation  spectrum  becomes  gapped  at  a  temperature  close  to  Tc  in  the  optimally 
doped  cuprates. 

The  lower  panel  of  fig.  6-9  shows  the  temperature  dependence  of  the  super- 
fluid  fraction  fs{T).  Above  Tc  fs  =  0,  and  increases  as  the  temperature  drops  below 
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Tc.  The  symbols  in  the  figure  are  the  superfluid  fractions  extracted  from  the  two- 
fluid  fits  to  the  transmittance  data  in  the  superconducting  state  of  both  films.  The 
solid  and  dashed  curves  are  fits  to  the  data  using  an  empirical  Gorter-Casimir  two- 
fluid  model  of  the  form  /.(T)  -  /,(0)(1  -  {T/TcY),  with     =  86  K.  The  fits  show 
that  even  at  low  temperatures  fs  is  still  much  less  than  unity,  with  /^(O)  =  0.66 
for  the  500  A  film  and  0.52  for  the  600  A  film.  It  means  that  a  sizeable  fraction  of 
the  carriers  do  not  condense.  This  is  a  universal  feature  of  all  the  cuprates.  This 
residual  conductivity  which  remains  at  low  temperatures  in  the  superconducting 
state  of  the  cuprates  is  still  the  subject  of  intense  debate. 

6.3  Summary 

We  have  presented  the  results  of  optical  reflectivity  and  transmittance  mea- 
surements of  thin  YBa2Cu306+x  films  in  the  far-infrared  region  below  350  cm~^ 
The  transmittance  of  both  films  was  fitted  with  a  two-component  model,  with 
a  Drude  free-carrier  response  in  the  normal  state,  and  a  two-fluid  model  in  the 
superconducting  state.  The  Drude  scattering  rates  extracted  from  the  fits  show  a 
T-linear  variation  in  the  normal  state  and  a  more  rapid  drop  just  below  Tc,  while 
the  Drude  plasma  frequency  is  nearly  temperature  independent.  It  was  found 
that  fairly  good  fits  could  be  obtained  by  fixing  the  midinfrared  absorption  at  all 
temperatures.  The  superfiuid  fraction  extracted  from  the  two-fiuid  fits  below  Tc 
increases  with  decreasing  temperature,  but  a  sizeable  fraction  of  the  carriers  do  not 
condense  even  at  40  K. 

The  fits  to  the  film  on  sapphire  were  complicated  by  the  optical  anisotropy 
of  the  sapphire  substrate.  Although  the  fits  to  the  film  data  were  reasonable,  the 
optical  anisotropy  of  the  substrate  leads  to  errors  in  the  optical  constants  extracted 
from  the  fits.  The  uniaxial  transmittance  model  which  was  used  to  model  the 
optical  anisotropy  rephcated  the  broad  features  of  the  substrate  data. 


Although  we  have  presented  the  results  of  the  far-infrared  investigations  of 
just  two  films;  the  results  are  consistent  with  similar  measurements  which  we 
performed  on  other  films.  A  total  of  three  films  on  MgO  and  5  films  on  sapphire 
were  measured.  The  generic  features  of  the  optical  characteristics  of  all  the  films 
investigated  were  similar  to  the  ones  presented  here. 


0.6 


Q) 

O  0.4 

c 
o 

o 

^  0.2 


0.0 

1.0 

Q) 
O 

C  0.8 
O 

■+-> 

0.6 

U) 

O  0.4 
i_ 
I- 

0.2 
0.0 


r  I  I  I  I  I  I  I  I  I  I  I  I  1  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

0.3  mm  thick  MgO  substrate 


  300  K 

  200  K 

  100  K 

  60  K 

  40  K 

I  I  I  I  I  I  I  I  I  [  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

Raw  data  and       . , .  o.   300  k 

-fits  ■  -  ^   200  K 

...  A   1  00  K 

...  .0.   75  K 

.    60  K 

...  .X.   40  K 


I  J  I  i  I  I  I  I 


100  200  300 

Frequency  (cm~^) 


400 


Figure  6-1:  The  smoothed  reflectances  and  transmittances  of  MgO.  Upper  panel: 
the  measured  reflectance  of  a  MgO  substrate  at  several  temperatures  between  40 
and  300  K.  Lower  panel:  the  soUd  curves  are  the  fits  to  the  transmittance  using  a 
series  of  Lorentzian  oscillators  for  the  dielectric  function,  while  the  spectra  in  dots 
and  symbols  are  the  smoothed  transmittance  data. 
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Figure  6-2:  The  frequency-dependent  reflectance  and  transmittance  of  cerium  oxide 
on  R-cut  sapphire  (0.5  mm  thick)  at  a  few  select  temperatures. 
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Figure  6-3:  The  measured  transmittance  of  Ce02  on  sapphire,  compared  to  the  fits 
with  a  uniaxial  transmittance  model  at  several  temperatures.  The  fits  are  repre- 
sented by  solid  lines  while  the  data  is  in  dots  and  symbols. 
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Figure  6-4:  Measured  reflectance  (upper  panel)  and  transmittance  in  dots  and 
symbols  (lower  panel)  of  a  500  A  YBa2Cu306+x  film  on  a  0.3  mm  thick  MgO  sub- 
strate at  several  temperatures  above  and  below  Tc.  The  solid  curves  in  the  lower 
panel  are  fits  to  the  measured  transmittance  using  a  two-component  model. 
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Figure  6-5:  The  frequency  dependent  reflectance  and  transmittance  of  a  600  A 
YBa2Cu306+x  film  on  a  sapphire  substrate  at  a  few  select  temperatures  above  and 
below  Tc. 


84 


0.5 
0.4 
0.3 
0.2 
0.1 
0.0 

O 

c 

R  0.2 


E 

CO 

c 

D 


-1 — I — I — I — I — rg — I — I — I — I — I — 1 — I — I — I — I — I — I — r-T — i — i — i — i — i — i — i  r 

600  A  YBa2Cu30  6+x/ sapphire 


n — I — I — r 


'2^^3^6+x/ 

raw  data  and  fits 


Q   300  K 

p.   200  K 

A  100  K 


I  I  I  I  I  I  I  I  I    I  I  I  I  I  I  I  I  I 


.0.   75  K 


0.1 


0.0 


0.2 


0.1 


0.0 


I  I  I  I  I  I  I    I  I  I  I  I  I  I  I  I    I  I  I 


,*   40  K 


-I — 1 — I — I — I — I  J  I  I  I  I  I  I  I  I  I  I  I  I  1  I  I  '  '  ■  I  I  ■  I 


I  .  .  .  . 


100  200  300 

Frequency  (cnn~^) 


Figure  6-6:  Raw  transmittance  data  (dots  and  symbols),  and  the  corresponding  fits 
(solid  lines)  at  several  temperatures  above  and  below  Tc  for  a  thin  YBa2Cu306+x 
film  on  a  sapphire  substrate.  A  uniaxial  transmittance  model  was  used  to  fit  the 
substrate  data,  and  a  two-component  model  was  used  for  the  film  data. 


85 


20000 


15000 


G  10000 


5000 


15000 


G  10000 


cN  5000 

b 


T — I — I — I — 1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 1 — I — I — I — I — I — r- 


500  A  YBagCujOg+x  film 


300  K 
200  K 
100  K 
75  K 
60  K 
40  K 


I  I  I.  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


I  I  I  I 


w 

\  ■. 
\  '• 


\ 

\  ■. 
\  • 

\  . 


-I  i  I   r  I   I   I   r   I   I  I 


100  200  300 

Frequency  (cm~^) 


Figure  6-7:  The  real  part  (upper  panel)  and  the  imaginary  part  (lower  panel)  of 
the  complex  conductivity  of  the  YBa2Cu306+x  film  on  MgO  at  several  tempera- 
tures. The  conductivities  were  extracted  from  the  fits  to  the  transmittance  data. 
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Figure  6-8:  The  real  part  (upper  panel)  and  imaginary  part  (lower  panel)  of  the 
complex  conductivity  of  the  YBa2Cu306+x  film  on  sapphire  at  several  tempera- 
tures. The  conductivities  were  extracted  from  the  fits  to  the  transmittance  data. 
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Figure  6-9:  The  temperature  dependent  scattering  rate  and  superfluid  fraction 
of  the  two  films  discussed  in  this  work.  Upper  panel:  the  temperature  depen- 
dent Drude  scattering  rate,  1/r,  extracted  from  the  fits  to  the  transmittance  (in 
symbols)  and  estimated  from  the  dc  resistivity  (solid  and  dashed  fines).  Lower 
panel:  the  superfluid  fraction  fs{T)  =  ns{T)/n  as  a  function  of  temperature. 
The  solid  and  dashed  curves  are  the  best  fits  to  the  empirical  two-fluid  expression 
fsiT)  =  A(0)(1  -  (T/TcY),  where  /,(0)  =  0.66  (500  A  film)  and  0.52  (600  A  film). 
The  symbols  are  the  fs{T)  extracted  from  the  two  fluid  fits  to  the  transmittance 
data. 
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Table  6-1:  Parameters  giving  the  best  fit  to  the  transmittance  of  MgO  at  40,  75, 
100,  200  and  300  K.  The  last  oscillator  with  plasma  frequency  UpTO,  center  fre- 
quency ojto,  and  damping  constant  ^to  is  a  single  TO  phonon  mode  in  MgO.  The 
remaining  parameters  are  used  to  simulate  the  phonon  difference  band  processes 
below  400  cm~^ 


300  K 

200  K 

100  K 

75  K 

60  K 

40  K 

c<;pi(cm' 

6.4000 

3.6510 

3.6510 

3.7415 

5.0472 

7.2062 

a;i(cm" 

43.4450 

39.4450 

39.4450 

104.2421 

104.2421 

66.9503 

7i(cm" 

50.0000 

28.0000 

28.0000 

19.3637 

69.0000 

86.9352 

20.2676 

17.0927 

6.5638 

8.0721 

1.5908 

1.7637 

a;2(cm" 

102.8884 

105.1004 

104.2421 

146.5179 

146.5179 

146.5179 

72  (cm" 

14.4145 

15.6302 

15.7608 

18.8501 

16.4382 

23.6710 

a;p3(cm' 

14.8581 

8.9464 

10.5360 

7.9131 

1.5209 

0.3057 

a'3(cm" 

114.4269 

118.5110 

146.5179 

121.8676 

121.8676 

121.8676 

73  (cm" 

14.5488 

12.3915 

65.6370 

526.3826 

17.86472 

2.2729 

a;p4(cm" 

14.7325 

10.2064 

6.9068 

1.3840 

0.0024 

0.0016 

u;4(cm" 

126.2698 

128.2030 

121.8676 

183.9975 

183.9975 

183.9975 

74  (cm" 

14.6504 

13.4457 

36.3383 

13.6810 

118.8062 

171.0809 

u;p5(cm" 

14.3312 

10.2513 

10.7502 

0.0021 

0.0017 

0.0022 

u;5(cm" 

137.8251 

139.0850 

183.9975 

243.1730 

243.1730 

243.1730 

75  (cm" 

16.8833 

14.7339 

82.4147 

162.7641 

144.6945 

194.1924 

t^p6(cm" 

17.4937 

10.3808 

14.6358 

12.4206 

11.1924 

9.6699 

a;6(cm~ 

149.9949 

150.1191 

243.1730 

292.1794 

292.1794 

292.1794 

76  (cm" 

21.1419 

15.2990 

155.0362 

28.9699 

23.8954 

19.9415 

a;p7(cm' 

29.8535 

20.0011 

19.5771 

a;7(cm" 

177.3559 

173.9728 

292.1794 

77  (cm" 

45.1724 

38.4291 

50.6296 

u;p8(cm" 

48.6799 

19.2685 

u;8(cm" 

256.7527 

212.5819 

78  (cm" 

98.4499 

71.1310 

u;pro(cm 

995.0596 

1006.7790 

1016.8710 

1067.0120 

1057.5600 

1022.4100 

a;xo(cm 

409.0000 

423.0998 

425.4178 

425.3572 

419.1264 

419.1264 

7ro(cm- 

5.7792 

6.2541 

1.6389 

1.2655 

0.9469 

0.8781 

=  3.01,  thickness  d  =  0.3  mm 


Table  6-2:  Classical  oscillator  parameters  for  cerium  oxide  between  300  and  40  K. 


300  K 

200  K 

100  K 

80  K 

40  K 

UpToicm 

1200 

1080 

1100 

1100 

1020 

272 

273 

273 

273 

273 

7To(cm- 

14 

10 

10 

10 

9 

Coo  =  5.31,  thickness  d  =  300  A 
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Table  6-3:  Classical  oscillator  parameters  for  sapphire  o-ray  fits  at  all  tempera- 
tures. The  extra  mode  corresponds  to  a  mode  introduced  at  209  cm~^  to  account 
for  the  extra  absorption  at  higher  temperatures. 


mode  type 

300  K       200  K       100  K        80  K         40  K 

extra  mode 

u;pi(cm  ^) 
a;i(cm~-') 
7i(cm-i) 

56.0000      30.500      20.0000         1  — 
209          209          209          209  — 
300          300          300          300  — 

phonon 

a;p2(cm"^) 
a;2(cm"^) 
72(cm~^) 

220.7640    218.9280    216.9101    216.9101  216.9101 
384.3000    385.0949    385.8426    385.8426  385.8426 
4.1000       1.5000       0.3000       0.1000  0.0300 

phonon 

<x;p3(cm~^) 
a;3(cm"^) 
73(cm-i) 

732.8450    729.5039    725.6770    725.6770  725.6770 
438.9000    439.4219    439.8057    439.8057  439.8057 
2.6330       1.3000       0.3000       0.1000  0.0300 

phonon 

a;p4(cm"^) 
u;4(cm~^) 
74(cm-^) 

980.8650    978.6532    975.8593    975.8593  975.8593 
568.2000    568.7346    569.0360    569.0360  569.0360 
6.8180       3.5000       0.3000       0.1000  0.0500 

phonon 

a;p5(cm"^) 
a;5(cm"^) 
75(cm~i) 

241.2680    220.9885    197.1609    197.1609  197.1609 
633.6000    634.6421    635.5070    635.5070  635.5070 
6.3360       4.0000       0.3000       0.1000  0.0500 

electronic 

u;p6(cm"^) 
a;6(cm"^) 
76(cm-i) 

69149       67546       70078       70078  70078 
85723       86559       87439       87439  87439 
0.8570       0.8656       0.8656       0.8656  0.8656 

electronic 

cc;p7(cm"^) 
cJ7(cm~^) 
77(cm-^) 

164651      164651      164651      164651  164651 
137621      137621      137621      137621  137621 
1.3760       1.3762       1.3762       1.3762  1.3762 

Coo 

1.2000       1.3000       1.3000       1.3000  1.3000 

thickness  d  =  0.5  mm 
Table  6-4:  Classical  oscillator  parameters  for  sapphire  e-ray  fits  at  all  temperatures 


mode  type 

300  K 

200  K 

100  K 

80  K 

40  K 

a;pi(cm  ^) 

1040.6190 

1036.4720 

1031.8750 

1031.8750 

1031.8750 

phonon 

a;i(cm~^) 

399.5000 

401.7800 

404.1800 

404.1800 

404.1800 

7i(cm-i) 

7.9900 

6.0000 

0.300 

0.1000 

0.0500 

a;p2(cm"^) 

729.4160 

734.3642 

739.5846 

739.5846 

739.5846 

phonon 

u;2(cm~^) 

584.0000 

586.3750 

588.8750 

588.8750 

588.8750 

72(cm"^) 

20.4400 

17.0000 

0.3000 

0.2000 

0.0500 

cjp3(cm"^) 

61000 

61373 

61765 

61765 

61765 

electronic 

a;3(cm"^) 

82201 

83037 

83917 

83917 

83917 

73(cm-i) 

0.8220 

0.8304 

0.8392 

0.8392 

0.8392 

a;p4(cm-^) 

165673 

165673 

165573 

165673 

165673 

electronic 

a;4(cm~^) 

135082 

135082 

135082 

135082 

135082 

74(cm-i) 

1.3510 

1.3508 

1.3508 

1.3508 

1.3508 

,  thickness  d  —  O.b  mm 
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Table  6-5:  Parameters  used  in  the  two-component  fits  to  the  transmittance  of  the 
film  on  MgO  at  three  temperatures  in  the  normal  state. 


mode 

300  K         200  K         100  K 

Drude 

-i      /        /                    1  \ 

l/r(cm  ^) 

8626.7960    8613.1590  8614.8920 
428.7089     296.9635  160.8391 

mid-ir 

7i(cm 

3000          3000  3000 
290            290  290 
250.7474     250.0612  255.9913 

mid-ir 

a;p2(cm"^) 
a;2(cm""^) 
72(cm-i) 

9800          9800  9800 
740            740  740 
1550           1550  1550 

chain  state 

t<;p3(cm~^) 
a;3(cm"^) 
73(cm-^) 

14500         14500  14500 
3300          3300  3300 
7500          7500  7500 

Table  6-6:  The  parameters  used  in  the  fits  to  the  transmittance  of  the  film  on  MgO 
at  3  temperatures  in  the  superconducting  state.  The  higher  midinfrared  oscilla- 
tor and  the  chain  state  from  table  6-5  have  not  been  repeated  here  for  the  sake  of 
brevity. 


mode 

75  K 

60  K 

40  K 

8666.3366 

8700.3219 

8697.9634 

/. 

0.2027 

0.5029 

0.6313 

l/r(cm-i) 

109.4391 

98.9587 

99.1675 

a;pi(cm"^) 

3000 

3000 

3000 

mid-ir 

u;i(cm~^) 

290 

290 

290 

7i(cm-i) 

267.8137 

256.9662 

256.9662 

Table  6-7:  Parameters  which  were  used  in  the  two-component  fits  to  the  transmit- 
tance of  the  film  on  sapphire  at  three  temperatures  in  the  normal  state. 


mode 

300  K 

200  K 

100  K 

Drude 

ujpD^cm  ^) 

5432.5140 

5434.6850 

5486.0090 

l/r(cm-i) 

363.5505 

255.8019 

150.2427 

u;pi(cm-^) 

3000 

3000 

3000 

mid-ir 

a;i(cm~^) 

290 

290 

290 

7i(cm"i) 

250 

250 

250 

a;p2(cm"^) 

9800 

9800 

9800 

mid-ir 

a;2(cm"^) 

740 

740 

740 

72(cm-i) 

1550 

1550 

1550 

a;p3(cm-i) 

14500 

14500 

14500 

chain  state 

a;3(cm~^) 

3300 

3300 

3300 

73(cm~^) 

7500 

7500 

7500 
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Table  6-8:  Paxameters  used  in  the  two-fluid  fits  to  the  transmittance  of  the  film  on 
sapphire  at  two  temperatures  in  the  superconducting  state.  The  two  midinfrared 
oscillators  and  the  chain  state  oscillator  which  was  used  in  these  fits  are  the  same 
as  those  in  table  6-7. 


75  K 

40  K 

Updicm  ^) 

5486.0085 

5701.8970 

/. 

0.2000 

0.4948 

l/r(cm-i) 

108.4216 

96.3746 

CHAPTER  7 

FAR-INFRARED  ELECTRO-OPTIC  STUDIES  OF  YBaaCuaOg+x  FILMS 

7.1  Introduction 

In  this  chapter  we  present  the  results  of  electro-optic  transmittance  measure- 
ments on  two  YBa2Cu306+x  thin  films  in  the  frequency  range  from  40-350  cm~\ 
and  at  two  temperatures,  75  K  and  40  K,  in  the  superconducting  state.  The  details 
of  the  experiment  have  been  described  in  chapter  5.  We  will  also  present  the  results 
of  some  simple  simulations,  using  a  simple  two-fluid  model,  of  the  changes  in  trans- 
mittance which  are  to  be  expected  if  the  main  effect  of  the  current  is  to  transfer 
spectral  weight  from  the  condensate  to  the  Drude  band.  This  will  be  followed  by  a 
discussion  of  the  critical  state  model  [89,  90]  which  is  often  used  to  study  critical 
currents  in  the  cuprates.  The  model  in  many  cases  has  given  fairly  good  agreement 
with  flux  and  current  density  distributions  determined  by  a  variety  of  imaging 
techniques  [91,  92].  This  discussion  will  facilitate,  a  basic  understanding  of  the  data 
in  some  cases,  particularly  for  currents  much  smaller  than  the  critical  current. 

The  effects  of  large  transport  currents  on  the  superconducting  state  of  the 
cuprates  have  been  studied  experimentally  by  dc  transport  [93-95],  angle  re- 
solved photoemission  spectroscopy  (ARPES)  [96]  and  tunneUing  [97].  The  main 
motivation  for  this  study  is  to  investigate  whether  or  not  the  changes  in  the  su- 
perconducting state  induced  by  a  transport  current  can  be  observed  in  the  far 
infrared,  which  probes  the  low  frequency  electrodynamics  of  the  charge  carriers. 
The  results  here  will  be  discussed  qualitatively.  There  are  two  main  obstacles  to 
any  quantitative  analysis  of  the  results:  the  first  difficulty  arises  from  the  absence 
of  a  complete  theory  of  superconductivity  in  the  cuprates,  and  the  second  is  related 
to  ambiguities  in  assigning  the  observed  changes  solely  to  the  eflFects  of  the  current. 
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These  issues  will  be  discussed  in  later  sections  where  we  will  address  the  problems 
associated  with  the  inherent  inhomogeneity  of  the  current  density  distribution  in 
the  superconducting  state  and  the  possible  spurious  effects  of  contact  heating  on 
the  observed  changes  in  the  transmittance. 

.  ■  7.2    Results  and  Discussion 

7.2.1    Electro-optic  Transmittance  of  Film  on  MgO 

The  results  of  our  electrooptic  transmittance  measurements  on  a  YBa2Cu306+x 
film  on  MgO  at  75  K  is  shown  in  fig.  7-1.  The  contact  arrangement  which  was 
used  for  the  transport  measurements  is  shown  in  fig.  5-3  (b)  of  chapter  5.  The 
fringes  in  the  measured  transmittance  in  part  (A)  of  the  figure  are  due  to  the 
multiple  internal  reflections  in  the  substrate,  which  have  been  smoothed  in  part 

(B)  .  At  each  current  we  measured  the  corresponding  voltage  and  the  transmittance. 
The  results  show  that  at  all  currents  up  to  200  mA  there  is  little  or  no  change  in 
the  transmittance.  At  all  current  levels  the  voltage  was  zero  which  implies  that 
the  film  remained  superconducting  over  the  entire  current  range.  This  insensitivity 
of  the  transmittance  to  the  applied  current  at  all  current  levels  up  to  200  mA,  is 
further  highlighted  by  the  fractional  change  in  the  transmittance  plotted  in  part 

(C)  of  the  figure.  The  fractional  change  in  transmittance  is  seen  to  be  zero  at  all 
currents.  We  also  attempted  to  measure  the  transmittance  at  300  mA.  At  this 
current  level  the  voltage  was  zero  at  the  start  of  the  transmittance  measurements, 
then,  the  film  switched  to  the  resistive  state  after  a  few  seconds  and  the  voltage 
increased  rapidly.  The  substrate  cracked  and  the  film  was  destroyed.  It  is  well 
known  from  previous  attempts  to  measure  the  effects  of  large  transport  currents 
on  the  electronic  density  of  states  of  conventional  superconductors  by  tunnelling 
[98],  that  it  is  possible  to  have  resistive  switching  of  a  film  carrying  a  transport 
current  even  below  the  critical  current  /g.  This  is  attributed  to  the  superconducting 
state  being  metastable  while  carrying  a  transport  current.  If  the  current  fiows  for 
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Figure  7-1:  The  transmittance  of  a  500  A  YBa2Cu306+x  film  as  a  function  of  cur- 
rent. (A)  The  raw  transmittance  T{I)  as  a  function  of  current  7    =    0,  25,  75  and 
200  mA.  The  voltage  V  is  zero  at  all  currents.  (B)  The  smoothed  data,  with  the 
fringes  in  (A)  at  all  currents  smoothed  out.  (C)  The  fractional  change  AT/T(0)  in 
transmittance  with  current,  expressed  as  a  fraction  of  the  zero-current  transmit- 
tance r(0),  where  AT  =  T{I)  -  T{0). 


a  long  enough  time  the  film  eventually  switches  to  the  lossy  resistive  state  even 
though  the  current  is  less  than  the  critical  value.  This  is  a  possible  explanation 
for  the  switching  observed  here  at  300  mA.  The  current  density  at  300  mA  is 
J  =  5  X  10^  Acm"^,  if  we  assume  a  uniform  current  density  distribution.  The 
second  possibility  is  related  to  ohmic  losses  in  the  contacts  with  the  film.  We 
estimate  the  contact  resistance  of  each  current  contact  to  be  about  0.1  Cl.  Although 
we  did  not  perform  a  detailed  measurement  of  the  contact  resistance  we  simply 
used  a  multimeter  to  measure  the  resistance  across  the  current  leads  while  the  film 
is  in  the  superconducting  state.  This  gives  the  resistance  of  the  lead  wires  and 
the  current  contacts.  We  estimated  the  lead  wire  resistance  by  shorting  the  leads 
and  measuring  their  resistance.  The  lead  resistance  was  then  subtracted  and  the 
resulting  resistance  was  taken  to  be  the  resistance  of  the  current  contacts.  For  a 
current  of  300  mA  flowing  through  a  0.2  Q,  resistor  the  dissipated  power  is  18  mW. 
This  power  is  large  enough  to  produce  significant  heating  over  the  duration  of  a 
typical  transmittance  measurement  which  is  anywhere  between  5  to  10  minutes. 
It  is  impossible  at  this  point  to  distinguish  which  process  is  responsible  for  the 
switching,  nonetheless,  it  appears  that  the  temperature  gradient  between  the 
surface  of  the  film  and  the  backside  of  the  MgO  substrate  caused  the  substrate  to 
crack. 

7.2.2    Electro-optic  Transmittance  of  Film  on  Sapphire 

The  upper  and  lower  panels  of  fig.  7-2  show  the  results  of  our  electrooptic 
transmittance  measurements  on  the  film  on  sapphire  at  75  K.  The  contact  arrange- 
ment used  for  the  transport  measurements  is  shown  in  fig.  5-3  (c)  in  chapter  5.  The 
spectra  show  that  there  is  little  or  no  change  in  the  measured  transmittance  up  to 
10  mA.  The  measured  voltage  at  10  mA  is  zero,  which  means  that  the  film  is  still 
superonducting  at  this  current  level.  This  insensitivity  of  the  measured  transmit- 
tance to  current  at  10  mA  is  most  easily  seen  in  the  plot  of  the  fractional  change 


in  the  lower  panel.  The  fractional  change  in  transmittance  is  zero  over  the  entire 
frequency  range  at  this  current.  The  transmittance  in  the  absence  of  the  current  is 
nearly  indistinguishable  from  the  10  mA  spectrum.  There  are  tiny  differences  which 
are  within  the  experimental  uncertainty  of  the  measurement.  A  very  different 
effect  is  observed,  however,  when  the  current  is  increased  to  15  mA.  The  measured 
voltage  at  this  current  level  was  323.16  mV  which  means  that  the  sample  is  in 
the  resistive  state.  The  critical  current  therefore  lies  somewhere  between  10  and 
15  mA.  This  implies  that  the  critical  current  density  for  this  film  is  in  the  range 
3.3  X  10^  and  5  x  10^  Acm~^.  It  is  interesting  to  compare  this  current  density  with 
the  critical  current  densities  in  refs.  [94]  and  [95].  The  critical  current  densities 
in  the  latter  reference  were  between  2  x  10^  and  2.3  x  10^  Acm"^  at  77  K.  The 
transition  temperatures  of  all  their  YBa2Cu306+x  films  were  between  91  and  93 
K.  They  used  a  5  /iVcm"^  criterion  to  determine  the  critical  current  density.  The 
current  densities  J*  in  ref.  [94]  at  which  there  were  full  transitions  to  the  normal 
state  were  always  between  10^  and  10^  Acm— 2.  The  critical  current  density  cor- 
responding to  the  transition  to  the  resistive  state  at  15  mA  in  this  film  is  about  2 
or  3  orders  of  magnitude  lower  that  the  corresponding  critical  currents  in  ref.  [94]. 
One  possible  reason  for  the  discrepancy  may  be  that  the  film  in  this  experiment 
was  at  a  temperature  much  closer  to  Tc  than  the  measured  temperature  of  75  K. 
There  are  two  notable  features  of  the  transmittance  observed  at  15  mA.  The  first  is 
an  increase  in  the  low  frequency  transmittance  below  100  cm~^  when  compared  to 
the  transmittance  without  current.  The  second  is  a  depression  of  the  transmittance 
above  100  cm~^  below  the  zero-current  spectrum.  The  fractional  change  at  this 
current  shows  that  the  low  frequency  transmittance  is  enhanced  by  about  15%  at 
40  cm~^  and  is  depressed  by  a  similar  percentage  at  350  cm"^.  It  is  interesting  to 
compare  the  resistivity  calculated  from  the  fits  to  the  film  transmittance  in  chapter 
6  at  75  k,  to  that  calculated  from  the  transport  measurements  here  at  15  mA.  The 
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resistivity  is  calculated  from  the  equation 


RA 

Pdc  =  —r~i 


(7.1) 


where  R  =  V/I  where  V  is  the  measured  voltage  and  I  is  the  current,  A  is  the 
cross  sectional  area  of  the  film,  and  /  is  the  center-to-center  distance  between  the 
voltage  electrodes.  The  corresponding  dc  resistivity  at  75  K  can  be  calculated  from 
the  fits  to  the  film  transmittance,  if  we  assume  that  all  the  superfluid  is  completely 
converted  to  normal  fluid  by  the  current  without  changing  the  scattering  rate.  The 
relevant  equation  is:     ,;  - 


where  1/r  is  the  scattering  rate  and  Upu  is  the  Drude  plasma  frequency.  The 
calculations  give,  pdc  =  215.44  iJ,Q  cm  and  pfu  =  216.15  /xfi  cm.  The  agreement 
here  is  quite  good  in  spite  of  the  the  errors  inherent  in  parameters  used  to  fit  the 
film  data  at  75  K.  If  this  truly  represents  the  normal  state,  we  expect  that  the 
resistivity,  and  hence  the  fractional  change  in  transmittance  at  currents  larger  than 
15  mA  should  saturate  (i.e.  the  resistivity  and  fractional  change  should  remain 
constant  at  the  15  mA  values  at  higher  currents).  The  transmittance  at  20  mA 
exceeds  that  at  15  mA  at  low  frequency  and  falls  below  it  at  higher  frequency. 
The  fractional  change  at  20  mA  in  the  figure  approaches  25%  at  40  cm~^.  This 
exceeds  the  fractional  change  observed  at  15  mA,  which  means  that  the  fractional 
change  does  not  saturate  at  the  15  mA  value.  The  calculated  resistivity  at  20  mA 
is  P(ic(20  mA)  =  237.33  nCt  cm.  This  resistivity  falls  on  the  finear  portion  of  the 
temperature  dependent  resistivity  plotted  in  the  lower  panel  of  fig.  5-1  in  chapter 
5.  These  observations  are  a  signature  of  heating  as  neither  the  resistivity  nor 
the  fractional  change  in  transmittance  saturate  at  their  15  mA  values.  It  is  well 
known  that  currents  and  increased  temperature  act  as  pair  breakers  and  weaken 
superconductivity.  The  effects  here  may  be  a  combination  of  both.  At  this  point  it 


Pfit  = 


47r(l/r) 


(7.2) 


is  impossible  to  know  how  much  of  the  observed  changes  are  due  to  either  process 
as  both  produce  similar  effects. 

The  40  K  electro-optic  spectra  shown  in  fig.  7-3  show  quite  different  re- 
sults. The  transmittance  was  measured  at  three  currents:  50,  80  and  100  mA. 
The  voltage  reading  in  each  case  was  zero  which  means  that  the  film  remained 
superconducting  up  to  100  mA.  We  also  attempted  to  measure  the  electrooptic 
transmittance  at  150  mA  but  during  the  experiment  the  voltage  suddenly  increased 
from  zero,  and  continued  to  increase  up  to  a  few  volts.  This  was  a  sign  of  heating 
and  the  transmittance  at  low  frequency  below  100  cm"^  approached  the  normal 
state  transmittance  at  room  temperature.  This  implies  that  the  critical  current 
density,  Jc  of  this  film  at  40  K  exceeds  5  x  10'*  Acm~^.  The  critical  current  densities 
of  the  YBa2Cu306+x  films  in  ref.  [95]  exceed  10^  Acm~^  at  60  K  and  is  of  the  same 
order  of  magniture  at  4.2  K.  It  is  absolutely  hopeless  to  axihieve  such  current  den- 
sities at  40  K  in  this  experiment  without  significant  heating.  The  critical  current 
Ic  would  be  of  the  order  of  a  few  Amperes.  We  have  not  included  the  spectrum 
at  150  mA  here  as  it  adds  nothing  to  the  current  discussion.  It  can  be  seen  in  the 
lower  panel  of  the  figure  that  the  fractional  change  in  the  transmittance  at  50  mA 
increases  just  above  4%  below  50  cm~\  and  decreases  at  350  cm"^  by  about  3%. 
At  80  mA  the  reduction  of  the  fractional  transmittance  above  300  cm~^  is  stronger, 
approaching  11%.  The  lower  frequency  firactional  transmittance  change  approaches 
15%  at  this  current  at  40  cm~^  The  100  mA  transmittance  clearly  shows  an  en- 
hancement below  100  cm~^  along  with  a  further  suppression  of  the  high  frequency 
transmittance.  It  is  not  at  all  obvious  if  the  changes  here  are  associated  with  the 
current  or  heating  due  to  losses  in  the  contacts.  It  is  possible  that  the  losses  in  the 
contacts  caused  a  temperature  rise  in  the  film,  and  that  after  a  time  the  film  came 
to  a  steady-state  temperature  above  40  K  in  the  superconducting  state.  This  would 
be  consistent  with  what  is  observed  in  the  electooptic  transmittance  measurements 
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Figure  7-2:  The  electro-optic  transmittance  of  a  thin  YBa2Cu306+x  film  on  sap- 
phire at  75  K.  Upper  panel:  The  measured  transmittance  T{I)  at  currents  7    =  0, 
10,  15  and  20  mA.  The  measured  voltages  V  are  also  shown  at  each  current.  Lower 
panel:  the  fractional  change  in  transmittance  AT/T(0)  at  each  current  expressed  as 
a  fraction  of  the  zero  current  spectrum  T(0),  where  AT  =  T{I)  -  r(0). 
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along  with  the  absence  of  a  voltage  at  any  of  the  measured  currents.  In  the  next 
section  we  shall  consider  some  simple  simulations  of  the  expected  changes  in  the 
transmittance  using  a  simple  two-fluid  model  of  the  superconducting  state  with  a 
superfluid  fraction  which  varies  with  the  current. 
7.2.3    Two-Fluid  Simulations 

Currents  and  magnetic  fields  modify  the  superconducting  state  of  the  cuprates 
by  breaking  apart  the  Cooper  pairs.  This  leads  to  a  weakening  of  superconduc- 
tivity. This  concept  has  been  quantified  by  Usadel,  using  a  single  parameter,  the 
depairing  energy  F,  in  the  diffusive  limit  in  conventional  BCS  superconductors  [99]. 
The  Usadel  equations  extended  the  equivalence  of  the  effect  of  fields  and  currents 
on  the  superconducting  state  to  inhomogeneous  cases  where  the  modulus  of  the 
order  parameter  varies  in  space.  The  Usadel  equations  have  been  used  to  analyze 
the  effects  of  supercurrent  on  the  density  of  states  of  a  superconductor  studied 
by  electron  tunneUing  [98].  Despite  this,  to  the  current  authors  knowledge  there 
is  no  corresponding  theory  for  the  effects  of  a  supercurrent  on  clean  limit  d-wave 
superconductors  such  as  the  cuprates. 

We  consider  here  a  simple  model  where  the  only  effect  of  a  supercurrent  is  to 
reduce  the  superfiuid  fraction  in  the  two-fluid  fits  to  the  transmittance  in  chapter  6. 
Let  us  consider  a  simple  two  fluid  dielectric  function  of  the  form: 

e{u)  -  (1  -  fs{I))eD{io)  +  fs{I)esiu;)  +  (T  <  T,),  (7.3) 

where  all  the  quantities  here  are  the  same  as  those  defined  in  eqns.  (6.4)  and 
(6.5)  in  chapter  6,  but  with  a  superfluid  fraction  /s(/)  which  is  a  monotonic 
decreasing  function  of  the  supercurrent  /.  In  this  simple  model  the  only  effect  of 
the  supercurrent  is  to  transfer  spectral  weight  from  the  condensate  to  the  Drude 
band  without  affecting  the  Drude  scattering  rate.  The  f  sum  rule  is  automatically 
satisfied  by  this  form  of  the  dielectric  function.  Although  this  simple  model  does 
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Figure  7-3:  The  electro-optic  transmittance  of  a  thin  YBa2Cu306+x  film  on  sap- 
phire at  40  K.  Upper  panel:  The  measured  transmittance  T{I)  at  currents  7   =  0, 
50,  80,  and  100  mA.  The  measured  voltages  V  shown,  are  zero  at  all  currents. 
Lower  panel:  the  firactional  change  in  transmittance  AT/T{0)  at  each  current  ex- 
pressed as  a  fraction  of  the  zero  current  spectrum  T{0),  where  AT  =  T{I)  -  T(0). 
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not  give  a  functional  form  of  the  superfluid  function  fs{I)  it  allows  us  to  calculate 
the  transmittance  of  the  film  on  its  substrate  for  various  values  of  fs{I).  For  the 
sake  of  notational  simplicity  we  will  denote  the  zero  current  superfluid  fraction 
/s(0)  =  fso-  In  the  two-fluid  fits  to  the  film  on  MgO,  and  sapphire  fso{7^K)  ^  0.2 
at  75  K.  The  superfiuid  fraction  at  40  K  for  the  film  on  sapphire  fsoi^OK)  ^0.5. 

Simulations  for  film  on  MgO.  The  top  panel  of  fig.  7-4  shows  the  plots  of  the 
simulated  transmittance  of  the  500  A  film  on  MgO  at  75  K,  calculated  for  various 
values  of  the  superfluid  fraction  /s(/).  The  parameters  used  for  the  zero  current 
spectrum  are  the  same  as  those  which  were  used  in  the  two-fluid  fits  in  chapter  6. 
The  superfiuid  fraction  is  decreased  by  25%,  50%,  75%  and  then  100%  from  the 
zero  current  value  fso  and  the  transmittance  is  calculated  for  each  case.  Within 
this  simple  model,  the  transmittance  at  low  frequency  below  100  cm~^  increases, 
while  the  transmittance  at  higher  frequencies  decreases  below  the  zero-current 
spectrum  as  the  superfluid  fraction  decreases.  The  changes  are  most  easily  seen  in 
the  lower  panel  of  fig.  7-4  where  the  corresponding  plots  of  the  fractional  changes 
in  transmittance  axe  shown.  In  all  cases  the  fractional  change  in  transmittance 
at  low  and  high  frequencies  is  always  less  than  less  than  12%.  The  /s(/)  =  0 
curve  would  correspond  to  the  spectrum  at  or  above  the  critical  current  Ic  when 
superconductivity  is  completely  destroyed.  Above  this  current  density,  when  all  the 
spectral  weight  of  the  condensate  has  been  completely  depleted,  it  is  not  expected 
that  the  current  wiU  further  change  the  transmittance. 

Simulations  for  film  on  sapphire.  The  simulated  changes  in  the  transmittance 
for  a  film  on  sapphire  at  75  and  40  K  are  shown  in  figs.  7-5  and  7-6  respectively. 
Although  the  upper  panels  of  both  figures  do  not  show  the  changes  very  clearly,  the 
lower  panels  show  the  plots  of  the  fractional  changes,  where  the  fringes  have  been 
smoothed.  There  are  a  number  of  features  which  are  common  to  the  simulations  at 
both  temperatures.  The  first  apparent  observation  is  an  increase  in  the  fractional 
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Figure  7-4:  Simulations  of  the  changes  in  the  transmittance  of  a  film  on  MgO  at 
75  K  with  decreasing  superfluid  fraction  /s(/).  The  calculations  were  done  using 
the  two-fluid  model,  and  assuming  that  the  only  effect  of  the  current  is  to  trans- 
fer spectral  weight  from  the  condensate  to  the  Drude  band.  (A)  The  calculated 
transmittance  at  several  values  of  the  superfluid  fraction  (B)  The  fractional 

change  in  transmittance  calculated  from  the  simulations  in  part  (A),  expressed  as  a 
fraction  of  the  zero-current  transmittance. 
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change  in  transmittance  at  very  low  frequency— below  50  cm~^  at  75  K  and  below 
60  cm"^  at  40  K.  The  fractional  increase  at  75  K  at  the  lowest  frequency  40  cm"^ 
never  exceeds  5%.  This  is  in  stark  contrast  to  the  much  stronger  enhancement 
of  the  transmittance  at  40  K  which  approaches  50%  at  40  cm"^  for  the  fso  =  0 
curve.  The  second  feature  of  the  simulations  is  the  reduction  of  the  transmittance 
at  frequencies  above  60  cm"^  with  decreasing  superfiuid  fraction.  The  suppression 
of  the  transmittance  at  40  K  is  much  stronger  than  the  corresponding  effect  at  75 
K  for  similar  fractional  changes  in  the  superfiuid  fraction. 

A  comparison  of  the  experimental  electro-optic  transmittance  spectra  at  75 
K  in  fig.  7-2  with  the  simulations  in  fig.  7-5  highlight  interesting  differences.  The 
greatest  disparities  occur  at  the  two  extremes  of  the  spectral  range  at  40  and  350 
cm~^  The  enhancement  of  the  low  frequency  transmittance  in  the  simulations 
are  of  the  order  of  3  to  5%  as  compared  to  the  electro-optic  data  where  the 
enhancement  is  between  15  and  25%.  The  suppression  of  the  transmittance  at 
350  cm"^  never  exceeds  4%  in  the  simulations,  but  approaches  30%  at  the  highest 
current  (20  mA)  in  the  actual  electro-optic  data.  If  the  simulations  are  correct 
it  means  that  these  discrepancies  may  be  related  to  the  effect  of  heating  on  the 
measured  spectra,  which  causes  the  substrate  transmittance  to  increase  at  low 
frequencies  and  decrease  at  higher  frequencies. 

There  is  much  closer  agreement  between  the  fractional  changes  in  the  trans- 
mittance of  the  measured  electro-optic  data  at  40  K  in  fig.  7-3  and  the  simulations 
in  fig.  7-6.  The  fractional  change  in  transmittance  at  100  mA  approaches  40%  at 
40  cm~\  this  is  very  close  to  the  simulated  fractional  change  at  0.5 fso-  The  reduc- 
tion of  the  transmittance  at  350  cm~^  is  about  18%  for  the  100  mA  spectrum  but 
only  about  5%  for  the  0.5 fso  simulation.  One  has  to  be  wary  of  trying  to  make  any 
numerical  comparison  between  the  data  and  the  simulations.  The  parameters  which 
were  used  to  fit  the  substrate  and  film  are  uncertain  and  may  lead  to  errors  in  the 
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simulated  changes.  The  80  mA  fractional  change  in  electro-optic  transmittance 
approaches  17%  at  40  cm"\  and  the  0.75 f so  simulated  change  is  about  20%.  There 
is  also  disagreement  between  the  latter  two  spectra  at  350  cm"^  These  differences 
may  also  be  a  result  of  contact  heating.  We  cannot  unambiguously  assign  the 
changes  observed  solely  to  the  effect  of  the  current.  In  the  next  section  we  will 
discuss  critical-state  theory  which  has  often  been  used  to  predict  current  density 
distributions  in  thin  superconducting  films.  It  will  highlight  some  of  the  difficulties 
in  quantifying  these  electro-optic  results. 

7.2.4    Current  Density  Distributions  in  Superconducting  Thin  Films 

The  most  notable  feature  of  any  superconductor  is  its  ability  to  carry  an 
electric  current  where  the  losses  are  reduced  to  an  unmeasurably  small  level.  This 
persists  up  to  some  maximum  loss-free  current  density  Jc,  beyond  which  the  usual 
losses  set  in.  The  loss-free  current  only  flows  in  a  small  surface  layer  in  the  super- 
conductor at  thermal  equilibrium.  The  phase  coherence  of  the  superconducting 
state  effectively  screens  fields  and  currents  from  the  interior  of  the  superconductor. 
This  phenomenon  of  diamagnetism  is  also  responsible  for  the  high  reflectivities  and 
low  transmittance  of  superconductors  which  are  observed  in  the  far  infrared. 

The  superconductor  establishes  a  Meissner  state  when  either  a  small  external 
magnetic  field  or  a  transport  current  is  applied.  In  this  state  the  external  fields 
and  currents  are  completely  screened  from  the  bulk  of  the  superconductor.  The 
screening  occurs  because  the  superconductor  sets  up  Meissner  supercvirrents  of 
density  Jm  which  fiow  in  a  surface  layer  of  thickness  of  the  order  A,  the  London 
penetration  depth.  The  maximum  Meissner  current  density  is  the  depairing  current 
density 
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Figure  7-5:  Simulations  of  the  changes  in  the  transmittance  of  a  film  on  sapphire 
at  75  K  with  decreasing  superfluid  fraction  /s(/).  (A)  The  calculated  transmit- 
tance at  several  values  of  the  superfluid  fraction  fs{I)-  (B)  The  fractional  change  in 
transmittance  calculated  from  the  simulations  in  part  (A),  expressed  as  a  fraction 
of  the  zero-current  transmittance.  The  fringes  have  been  smoothed  to  make  the 
changes  clearer. 
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Figure  7-6:  Simulations  of  the  changes  in  the  transmittance  of  a  film  on  sapphire 
at  40  K  with  decreasing  superfluid  fraction  (A)  The  calculated  transmit- 

tance at  several  values  of  the  superfluid  fraction  (B)  The  fractional  change  in 

transmittance  calculated  from  the  simulations  in  part  (A),  expressed  as  a  fraction 
of  the  zero-current  transmittance.  The  fringes  have  been  smoothed  to  make  the 
changes  clearer. 


where  the  kinetic  energy  of  the  supercurrents  balances  the  condensation  energy 
of  the  superconducting  state.  Here  0o  =  h/2e  is  the  flux  quantum,  A  is  the 
penetration  depth  and  ^  is  the  Ginzburg- Landau  coherence  length. 

If  the  magnetic  field  H  at  the  surface  of  a  type  II  superconductor  exceeds 
the  lower  critical  field  Hd  the  superconductor  lowers  its  free  energy  by  admitting 
quantized  magnetic  flux  or  vortices.  The  entry  of  vortices  into  a  superconductor 
leads  to  a  mixed  state  where  small  quasi-normal  regions  in  the  vortex  cores  and 
superconducting  regions  exist  in  thermal  equiUbrium.  The  macroscopic  magnetic 
flux  density  B  averaged  over  a  few  intervortex  spacings  has  magnitude  B  =  n^o 
where  n  is  the  local  vortex  density.  The  flux  lines  which  enter  the  superconductor 
are  pinned  by  inhomogeneities  in  the  material.  If  the  local  current  density  J 
exceeds  the  critical  current  Jc  the  flux  lines  depin  moving  under  the  action  of 
the  Lorentz  force  F  =  J  X  B/c.  This  flux  motion  induces  an  emf  and  leads  to 
dissipation.  If  the  current  is  a  transport  cmrrent  which  is  produced  by  an  external 
source,  a  stationary  dissipative  state  is  formed.  This  Jc  represents  the  critical 
current  which  is  usually  probed  by  nonlinear  transport  measurements  which  are 
used  to  plot  the  I-V  characteristics  of  the  superconducting  state.  The  depairing 
current  density  Jo  which  is  associated  with  the  full  transition  to  the  normal 
state  forms  an  upper  bound  on  the  magnitude  of  the  depinning  critical  current 
density  Jc-  It  is  expected,  at  least  intuitively,  that  the  most  drastic  changes  in  the 
electrooptic  transmission  should  occur  at  current  densities  close  to  Jq. 

We  now  consider  the  current  density  distribution  predicted  by  the  critical-state 
model  in  a  simple  geometry  where  exact  analytic  solutions  are  available.  Fig.  7-7 
shows  an  infinitely  long  thin  film  of  width  2W  and  thickness  d  carrying  a  transport 
current  It  in  the  superconducting  state.  Here  It  flows  in  the  y  direction  and  the 
current  density  vector  has  only  one  component  in  the  y  direction  Jy{x).  It  has 
been  shown  that  when  It  is  very  weak  it  generates  a  self  field  that  curves  around 
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Figure  7-7:  Schematic  diagram  of  the  thin  film  geometry  considered  here.  The  film 
has  width  2W  and  its  thickness  d  is  smaller  than  the  London  penetration  depth  A. 
The  film  is  infinite  in  the  Y  direction  where  a  transport  cm'rent  Jy{x)  flows  in  the 
superconducting  state. 

the  sample  so  that  there  is  a  discontinuity  in  the  magnetic  induction      in  the 
X  direction  between  the  top  and  bottom  surfaces.  The  boundary  conditions  then 
mandate  that  currents  flow  across  the  entire  width  of  the  film  not  just  near  the  film 
surface  within  a  penetration  depth  A.  The  strip  is  then  in  the  Meissner  state,  and 
no  magnetic  flux  has  penetrated  so  that  B^ix)  =  0  for  -W  <  x  <W.  The  analytic 
solutions  for  the  macroscopic  current  density  and  magnetic  induction  Bz{x)  has 
been  obtained  by  conformal  mapping  techniques  and  are  given  by: 

It 


Jy{x)  = 


B,{x)  - 


T^dyJW^-x"^' 
2xIt 


-W  <x<W, 


\x\  >  W. 


(7.5) 


c|3;|\/x2  —  VK^' 

At  higher  currents  when  the  B  field  near  the  edges  of  the  film  exceeds  the  lower 
critical  field  Hd  vortices  begin  to  penetrate.  The  depth  of  penetration  depends 
on  the  critical  current  density  Jc  which  is  assumed  to  be  independent  of  B.  The 


(7.6) 
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solutions  for  the  critical  state  which  is  established  after  the  flux  penetration  is: 


,   ^arctan  (  ^/^^-^  )  ,    -a  <  x  <  a,  ^  ^ 

Jy[x)  =  <[    ^  ;  (7.7) 

Jc,  a  <  |x|  <  W, 

and  the  width  of  the  field  free  region  2a  is  given  by 


a  =  -  {IrlhY  (7.8) 

where  Ic  =  2dWJc  is  the  strips  critical  current.  The  flux  density  distribution 
becomes: 

0,  —a  <  X  <  a, 

±^/l-7W^^'  a<|x|<l^,  (7.9) 
±5/ln^,^^,  \x\>W, 
where  Bf  =  ^dJc-  The  strip  is  therefore  composed  of  two  regions:  a  vortex  free 
central  region  of  width  2a  and  a  flux  penetrated  region  on  either  side  of  the  flux 
free  region  of  width  W  —  a.  It  follows  that  as  the  current  increases  up  to  Ic  the 
width  of  the  field-free  central  region  diminishes.  The  current  density  in  the  flux 
penetrated  regions  is  Jc  but  Jy{x)  <  Jc  in  the  central  region.  As  the  transport 
current  is  increased  beyond  Ic  the  vortices  begin  to  move  and  dissipation  sets  in. 

It  is  not  very  obvious  at  this  point  how  this  inhomogeneous  current  and 
flux  density  distribution  affects  the  infrared  electrooptic  transmittance.  It  is 
not  believed  that  the  self  field  of  the  current  which  produces  vortices  in  the  flux 
penetrated  regions  will  produce  noticeable  changes  in  the  measured  electrooptic 
spectra.  The  main  reason  for  this  is  that  the  magnetic  field  in  the  fiux-penetrated 
region  is  of  the  order  of  a  few  Gauss  which  is  quite  small.  In  most  previous  studies 
of  the  effect  of  magnetic  fields  on  the  infrared  spectra  of  the  cuprates  [100-104], 
changes  in  the  spectra  are  observed  only  at  very  high  fields  of  a  few  Tesla. 


7.3  Summary 

The  results  of  our  electro-optic  transmittance  measurements  on  a  YBa2Cu306+x 
film  on  MgO  at  75  K  and  on  a  similar  film  on  sapphire  at  75  and  40  K  show  differ- 
ing results.  The  film  on  MgO  showed  no  electro-optic  effect  up  to  a  current  level 
of  200  mA.  This  was  expected  from  the  fact  that  the  measured  voltage  remained 
zero  at  all  currents  up  to  this  value  which  meant  that  the  current  had  httle  effect 
on  the  superconducting  state.  At  300  mA  the  film  displayed  switching  but  the 
underlying  MgO  substrate  cracked.  The  cracking  may  be  related  to  heating  within 
the  film  and  contacts  which  created  a  thermal  strain  on  the  substrate.  The  fairly 
high  normal-state  resistivity  of  the  cuprates  means  that  heating  will  always  be  a 
problem  when  the  system  switches  to  the  resistive  state. 

In  contrast,  the  film  on  sapphire  showed  an  electro-optic  effect  at  15  mA  and 
at  75  K.  The  measured  voltage  at  this  current  was  nonzero  which  means  the  film 
was  no  longer  superconducting.  We  believe  that  this  change  in  the  transmittance 
observed  at  15  mA  is  spurious  and  may  be  caused  by  heating.  At  20  mA  the 
fractional  changes  in  the  transmittance  did  not  saturate  and  the  voltage  increase 
exceeded  that  which  would  be  expected  for  the  normal  state  resistivity  at  that 
temperature.  This  gave  further  evidence  that  the  results  are  due  to  heating  and 
not  the  pair  breaking  effect  of  the  current.  The  40  K  electrooptic  transmittance 
spectra  showed  changes  at  currents  where  the  measured  voltage  was  zero.  These 
may  be  true  changes  associated  with  the  current,  but  they  may  also  be  an  artifact 
of  the  losses  in  the  surface  contacts.  The  larger  currents  which  were  passed  at  40 
K  as  compared  to  those  at  75  K  made  such  losses  more  likely  even  for  contacts  of 
resistance  of  the  order  of  1  Q.  The  ensuing  temperature  rise  could  possibly  lead  to 
an  effect  similar  to  that  of  the  current. 

The  simulations  using  two-fluid  fits  to  the  film  on  MgO  from  chapter  6  as 
starting  parameters,  with  a  current  dependent  superfluid  fraction  gave  valuable 


insights.  The  main  result  is  that  a  decreasing  superfluid  fraction  would  lead  to  an 
increase  in  the  low  frequency  transmittance  below  50  cm~^  and  a  depression  of  the 
transmittance  at  higher  frequencies  up  to  350  cm~^  Because  we  were  unable  to 
observe  an  electro-optic  effect  in  this  film  we  were  unable  to  make  a  comparison  of 
the  measured  data  with  the  simulations. 

We  compared  the  model  simulations  for  the  transmittance  of  the  film  on 
sapphire  at  75  and  40  K  with  the  electro-optic  data  at  the  same  temperatures. 
There  were  stark  differences  between  the  simulations  and  data  at  75  K.  The 
fractional  changes  in  the  actual  data  far  exceeded  the  fractional  changes  in  the 
simulations  at  75  K.  This  may  be  a  sign  of  heating  in  the  films.  The  simulations 
and  electro-optic  data  agree  much  better  at  40  K.  The  fractional  changes  at  40 
cm~^  and  at  80  and  100  mA,  differ  from  the  simulations  at  0.75/so  and  0.5/so 
by  just  a  few  percent.  The  difference  between  fractional  changes  between  these 
spectra  at  350  cm~^  are  much  larger.  This  may  be  an  artifact  due  to  heating  in  the 
contacts  or  due  to  errors  in  the  parameters  used  to  fit  the  film  and  substrate  in  this 
spectral  region. 

A  consideration  of  the  the  current  density  distribution  in  a  thin  superconduct- 
ing strip  using  the  critical-state  model  shows  that  the  current  density  distribution 
is  inhomogeneous  at  sub-critical  currents.  This  inhomogeneity  results  from  the 
diamagnetic  nature  of  the  superconducting  state  which  tends  to  screen  fields  and 
currents.  This  means  that  any  quantitative  analysis  of  electro-optic  transmittance 
data  at  sub-critical  currents  will  be  quite  complicated.  The  mixed  state  which  is 
formed  at  /  >  /c  will  lead  to  small  losses.  In  this  simple  experiment  we  cannot 
account  for  these  losses  which  heat  the  film  and  substrate.  The  very  long  mea- 
surement times  required  for  spectroscopic  investigations  means  that  even  a  small 
amount  of  dissipation  leads  to  a  pileup  of  heat  after  long  times.  This  heating  will 
invariably  affect  the  electro-optic  data. 


CHAPTER  8 
CONCLUSION 

We  have  studied  the  the  effect  of  a  transport  current  on  the  superconducting 
state  optical  transmittance  of  two  thin  YBa2Cu306+x  films  over  a  frequency 
range  from  40-350  cm~^  in  the  far  infrared.  The  films  were  first  characterized  by 
fitting  the  the  measured  transmittance  at  several  temperatures  in  the  normal  and 
superconducting  states,  with  a  simple  two-component  model.  The  electro-optic 
transmittance  was  then  examined  near  Tc  at  75  K  for  both  films  and  in  the  low 
temperature  regime  at  40  K  for  the  film  on  sapphire. 

The  results  of  the  two-component  fits  in  the  normal  state  show  that  the  strong 
temperature  dependence  of  the  transmittance  of  both  films  is  related  to  the  strong 
temperature  dependence  of  the  free  carrier  Drude  conductivity  centered  at  zero 
frequency.  The  mid  infrared  band  at  higher  frequency  is  fairly  inert  and  insensitive 
to  temperature.  There  is  close  agreement  between  the  linear  dc  resistivity  measured 
by  dc  transport  and  the  Drude  conductivity  cri£)(a;  =  0,  T)  extracted  from  the  fits 
and  extrapolated  to  zero  frequency.  The  fairly  narrow  Drude  peak  at  u;  =  0  has 
a  width  of  the  order  keT  which  means  the  Drude  scattering  rate  l/r  decreases 
linearly  with  decreasing  temperature.  The  Drude  plasma  frequency  Upp  was  found 
to  be  nearly  temperature  independent  for  both  films.  Within  this  simple  two- 
component  model  the  strong  temperature  dependence  of  the  optical  properties 
of  the  films  is  assigned  to  a  strong  temperature  dependence  of  the  free  carrier 
relaixation  rate.  The  origin  of  the  mid  infrared  absorption  is  still  a  subject  of 
intense  debate. 

The  supercondcting  state  response  was  analyzed  within  the  framework 
of  a  simple  two-fluid  model.  It  was  found  that  the  superconducting  state  was 

113 


114 

marked  by  a  transfer  of  spectral  weight  from  the  Drude  band  to  the  origin,  which 
formed  a  delta  function  in  the  real  part  of  the  conductivity.  This  delta  function 
response  of  the  condensate  had  a  marked  effect  on  the  measured  transmittance. 
The  slope  of  the  transmittance  changed  below  Tc  and  approached  zero  in  the 
dc  Umit.  The  strength  of  the  suppression  of  the  low  frequency  transmittance 
increased  with  decreasing  temperature  down  to  40  K.  The  most  striking  feature  of 
the  superconducting  state  is  the  observation  of  a  depression  of  the  quasiparticle 
scattering  rate  below  the  linear  normal-state  trend.  The  low  temperature  fits  also 
demonstrate  that  a  large  fraction  of  the  carriers  do  not  condense  at  40  K.  The  two 
latter  observations  are  generic  features  of  the  superconducting  state  of  all  cuprates. 

The  electrooptic  transmittance  of  both  films  appeared  to  be  compUcated 
by  the  effects  of  heating.  The  film  on  MgO  had  a  critical  current  Ic  in  excess 
of  0.3  A  at  75  K.  We  did  not  observe  any  change  in  the  measured  electro-optic 
transmittance  up  to  0.2  A.  The  film  and  substrate  cracked  at  a  current  level  of 
0.3  A,  this  effect  is  believed  to  be  a  sign  of  contact  heating.  The  critical  current 
of  the  film  on  sapphire  was  between  10  and  15  mA  at  75  K.  The  electro-optic 
transmittance  showed  no  change  up  to  10  mA,  but  showed  slight  changes  at  15 
and  20  mA  where  the  film  was  in  the  resistive  state.  There  is  ample  evidence 
that  heating  played  a  role  in  these  changes.  The  first  observation  is  the  lack  of 
saturation  of  the  changes  in  the  electro-optic  transmittance  up  to  20  mA.  The 
second  is  related  to  the  large  differences  in  the  resistivity  calculated  from  the 
two-fluid  fits  and  from  the  transport  data  at  15  and  20  mA.  Then  too,  the  low 
temperature  electrooptic  transmittance  data  at  40  K  showed  changes  at  80  and 
100  mA  even  though  we  were  well  below  the  critical  current  Ic  of  the  film.  We  are 
reluctant  to  assign  these  changes  exclusively  to  the  depairing  effect  of  the  current. 
Although  the  effects  of  contact  heating  are  not  apparent  at  this  temperature  it 
could  lead  to  similar  changes. 
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We  have  also  discussed  the  difficulties  associated  with  carrying  out  a  quan- 
titative analysis  of  electro-optic  data.  The  discussion  of  the  critical-state  model 
shows  that  the  current  and  flux  density  distributions  in  the  superconducting  state 
are  quite  complex.  It  is  quite  difficult  to  analyze  the  results  without  a  suitable 
theory  of  superconductivity  of  the  cuprates.  The  simple  simulations  of  the  changes 
in  the  transmittance  for  the  film  on  MgO  and  on  sapphire  showed  what  would  be 
expected  if  the  only  effect  of  the  current  is  to  transfer  spectral  weight  from  the 
condensate  to  the  Drude  band.  The  simulations  for  the  film  on  sapphire  at  40  K 
compared  favorably  with  the  electro-optic  data  at  the  same  temperature.  The  large 
differences  between  the  simulations  and  electro-optic  data  at  75  K  were  probably  a 
sign  of  heating. 

It  is  worth  mentioning,  that  it  has  been  proposed  that  high  temperature  su- 
perconductors carrying  a  supercurrent  axe  suitable  candidates  for  use  as  optical 
modulators  [105].  Optical  modulators  are  key  components  of  fiber  optic  communi- 
cations technology.  The  idea  is  that  by  using  pulses  of  current  the  films  can  easily 
be  switched  from  an  opaque  state  to  a  highly  transmissive  state.  For  a  supercon- 
ducting modulator  of  this  kind  one  measure  of  the  modulator's  performance  is  its 
extinction  ratio.  The  extinction  ratio  is  the  intensity  of  light  transmitted  in  the  re- 
sistive state  to  the  intensity  of  light  transmitted  in  the  superconducting  state.  Most 
current  commercial  modulators  have  an  extinction  ratio  between  6:1  and  100:1. 
The  extinction  ratios  of  the  films  in  this  study  never  exceed  1.5  if  the  simulations 
are  correct.  Prom  this  study  it  appears  that  such  devices  would  offer  no  advantages 
over  conventional  modulators.  The  two  main  disadvantages  are  the  small  changes 
in  transmittance  induced  by  the  current  and  the  high  critical  current  densities  and 
high  normal  state  resistivity  of  the  cuprates  which  lead  to  sizeable  heating. 
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Optical  modulators  are  important  elements  in  fiber  optic  communications. 
Fiber  optic  cables  have  been  known  to  be  superior  to  conventional  cables  as 
conduits  of  information.  A  modulator  becomes  more  or  less  transparent  to  light 
when  driven  by  an  external  current. 

We  have  investigated  the  feasibility  of  using  superconducting  films  as  modula- 
tors of  infrared  hght.  This  investigation  has  led  us  to  the  conclusion  that  the  use 
of  such  films  to  build  such  a  device,  is  not  a  feasible  prospect  at  this  time.  There 
are  various  practical  problems  associated  with  these  proposed  modulators  which  are 
unavoidable.  It  was  also  shown  that  optical  modulators  based  on  superconductors 
would  offer  no  advantage  over  conventional  modulators  which  are  currently  in 
operation. 


